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PREFACE 


The following pages embody, substantially as delivered, a 
course of fifteen lectures given during the Spring of 1925 at the 
Moore School of Electrical Engineering of the University of 
Pennsylvania. 

After a brief introduction to the subject of electric circuit 
theory, the first five chapters are devoted to a systematic and 
fairly complete exposition and critique of the Heaviside Oper- 
ational Calculus, a remarkably direct and powerful method 
for the solution of the differential equations of electric circuit 
theory. 

The name of Oliver Heaviside is known to engineers the 
world over: his operational calculus, however, is known to, 
and employed by, only a relatively few specialists, and this 
notwithstanding its remarkable properties and wide applicability 
not only to electric circuit theory but also to the differential 
equations of mathematical physics. In the writer’s opinion 
this neglect is due less to the intrinsic difficulties of the subject 
than to unfortunate obscurities in Heaviside’s own exposition. 
In the present work the operational calculus is made to depend 
on an integral equation from which the Heaviside Rules and 
Formulas are simply but rigorously deducible. It is the hope 
of the writer that this mode of approach and exposition will be 
of service in securing a wider use of the operational calculus by 
engineers and physicists, and a fuller and more just appreciation 
of its unique advantages. 

The second part of the present work deals with advanced 
problems of electric circuit theory, and in particular with the 
theory of the propagation of current and voltage in electrical 
transmission systems. It is hoped that this part will be of 
interest to electrical engineers generally because, while only a 
few of the results are original with the present work, most of the 
transmission theory dealt with is to be found only in scattered 
memoirs, and there accompanied by formidable mathematical 
difficulties. 
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While the method of solution employed in the second part is 
largely that of the operational calculus, I have not hesitated to 
employ developments and extensions not to be found in Heavi- 
side. For example, the formulation of the problem as a Poisson 
integral equation is an original development which has proved 
quite useful in the actual numerical solution of complicated 
problems. The same may be said of the Chapter on Variable 
Electric Circuit Theory. 

In view of its two-fold aspect this work may therefore be 
regarded either as an exposition and development of the oper- 
ational calculus with applications to electric circuit theory, or 
as a contribution to advanced electric circuit theory, depending 
on whether the reader’s viewpoint is that of the mathematician 
or the engineer. 

In the text I have not made reference to the literature of the 
subject, now fairly extensive, except where the treatment in 
the text is inadequate and I wish to refer the reader to a fuller 
and more satisfactory discussion. In the biblhography, how- 
ever, there is furnished a list of original papers and memoirs. 
While this bibliography includes all pertinent papers known 
to the author, no claim to completeness can be made and it is 
possible that important papers have been omitted. The author 
would be grateful to any reader calling his attention to such 
omissions. 

It remains a pleasurable duty to acknowledge the assistance 
I have received in the preparation of this book. To Dr. Harold 
Pender, Dean of the Moore School of Electrical Engineering of 
the University of Pennsylvania I owe the invitation to deliver 
the lectures out of which this book grew, and stimulating interest 
and advice throughout their course. To my associate Mr. R. 8. 
Hoyt I am indebted for a careful reading of the manuscript, 
the verification of most of the formulas, and valuable con- 
tributions and advice. For example, the extremely compact 
proof of formula (29) is a modification and improvement due 
to Mr. Hoyt of an earlier proof. To my assistant, Miss Dorothy 
Angell are due all numerical calculations and curves and the 
careful proofreading of the entire book. 

Finally to the Officials of the American Telephone and 
Telegraph Company and in particular to Mr. O. B. Blackwell, 
Transmission Development Engineer, I would express appre- 
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ciation of their generous attitude toward the research under- 
lying this book and my indebtedness for the opportunity and 
encouragement to put this work into shape for publication. 
JoHn R. Carson. 
New Yorks, N. Y. 
October, 1926. 
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ELECTRIC CIRCUIT THEORY 
AND THE OPERATIONAL 
CALCULUS 


CHAPTER I 
THE FUNDAMENTALS OF ELECTRIC CIRCUIT THEORY 


While a knowledge, on the reader’s part, of the elements of 
electric circuit theory will be assumed, it seems well to start 
with a brief review of the fundamental physical principles of 
circuit theory, the mode of formulating the equations, and 
some general theorems which will prove useful subsequently. 

Electric Circuit Theory is that branch of general electro- 
magnetic theory which deals with electrical oscillations in linear 
electrical networks. An electrical network is a connected set of 
circuits or meshes each of which can be regarded as made up of 
resistance, inductance, and capacitance elements. If the number 
of independent closed meshes or circuits is finite, the network 
is said to have a finite number of degrees of freedom equal to 
the number of meshes, and the problem is mathematically 
formulated by a finite number of simultaneous linear differ- 
ential equations. If the number of closed circuits or meshes 
is made infinite and the size of the elements infinitesimal the 
problem is formulated by a partial differential equation, such 
as the telegraph equation. 

There are two distinct problems involved in electric circuit 
theory. First, the derivation of the differential equations from 
the physical laws and the constants and connections of the 
system, that is, the mathematical formulation of the physical 
problem: and secondly the solution of the differential equations. 

The equations of circuit theory may be established in a number 
of different ways. For example, they may be based on Maxwell’s 
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dynamical theory. In accordance with this method, the net- 
work forms a dynamic system in which the currents play the 
role of velocities. If we therefore set up the expressions for the 
kinetic energy, potential energy and dissipation, the network 
equations are deducible from general dynamic equations. 

The simplest, and for our purposes, a quite satisfactory basis 
for the equations of circuit theory are found in Kirchhoff’s 
Laws. These laws state that 

1. The total impressed force taken around any closed loop 
or circuit in the network is equal to the potential drop due to 
(a) resistance, (b) inductive reaction and (c) capacitive reactance. 

2. The sum of the currents entering any branch point in the 
network is always zero. 

Let us now apply these laws to an elementary circuit in order 
to deduce the physical significance of the circuit elements. 

Consider an elementary circuit consisting of a resistance 
element f, an inductance element ZL and a capacity element 
C in series, and let an electromotive force E be applied to this 
circuit. If J denote the current in the circuit, the resistance 
drop is RI, the inductance drop is LdI/dt and the drop across 
the condenser is Q/C where Q is the charge on the condenser. 
It is evident that Q and J are related by the equation J =dQ/dt 
or Q= fIdt. Now apply Kirchhoff’s law relating to the drop 
around the circuit: it gives the equation 


RI+LdI/dt+Q/C=E. 
Multiply both sides by [: we get 


Gm or a We 

oe apogee 

The nght hand side is clearly the rate at which the impressed 
force is delivering energy to the circuit, while the left hand side 
is the rate at which energy is being absorbed by the circuit. 
The first term RJ? is the rate at which electrical energy is being 
converted into heat. Hence the resistance element may be 
defined as a device for converting electrical energy into heat. 

ay a : ; ‘ 

The second term 2 LI* is the rate of increase of the magnetic 
energy. Hence the inductance element is a device for storing 


Ble 


energy in the magnetic field. The third term £ Q?/2C is the 
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rate of increase of the electric energy. Hence the condenser is a 
device for storing energy in the electric field. 

In the foregoing we have isolated and idealized the circuit 
elements. Actually, of course, every circuit element dissipates 
some energy in the form of heat and stores some energy in the 
magnetic field and some in the electric field. The analysis of the 
actual circuit element, however, into three ideal components is 
quite convenient and useful, and should lead to no misconception 
if properly interpreted. 

Now consider the general form of network possessing n inde- 
pendent meshes or circuits. Let us number these from 1 to n, 
and let the corresponding mesh currents be denoted by J,, 
Iz... .I,. Let electromotive forces #,, FE. .... E, 
be applied to the » meshes or circuits respectively. Let L;;, 
R;;, C;; denote the total inductance, resistance and capacity in 
series in mesh 7 and let Lyx, Rix, Cj, denote the corresponding 
mutual elements between circuits 7 and k. Now write down 
Kirchhoff’s equation for any circuit or mesh, say mesh 1; it is 


d 1 d 1 
(Un qth g,, {att (La dt tC, rot 


d 1 
+ (Engg Rite, { dU \Ta=Bs 


Corresponding equations hold for each and every one of the n 
meshes of the network. Writing them all down, we have the 
system of equations 


d il d 1 2 
(Lngthat g. [ait ohh + (Lng thing, | a o= Bs 
RRM Rn SO 2B mr ata eras Neg oh nine nven haciahs (1) 


d 1 d Were 
(Lng ttn i ait 3 + (Lon gt Bent a | dt \In=By 


The system of simultaneous differential equations (1) 
constitutes the canonical equations of electric circuit theory. 
The interpretation and solution of these equations constitute 
the subject of Electric Circuit Theory, and it is in connection 
with their solution that we find the most direct and logical 
introduction to the Operational Calculus. 

As an example of the appropriate mode of setting up the 
circuit equations, consider the two mesh network shown in Fig. 
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1. Writing down Kirchhoff’s Law for meshes 1 and 2, respec- 
tively, we have 


d 1 Ce 
(ngtmtg, | it\L+ MG ha Bs 
d d il 
My bt (lagthatg, i iU\h=Be 


In this case the self and mutual coefficients are given by 


Dy=Ly Lo. =Le Dy. =Lei= +M 
Cu=Ci Cor =Ce Ci2= C21 =0 
Rui=Rh, Roe= Re Ry.= Ro =0. 


The conventions adopted for the positive directions of currents 
and voltages are indicated by the arrows. The sign of the 
mutual inductance M will depend on the relative mode of 
winding of the two coils. 

Now write down Kirchhoff’s Law, or the circuital equations for 
the network of Fig 2. They are 


1 | 
EDN, C ) fa I, 
-(z StRste | @)e=2o 
—(Eoiy ee a.| wh 


sli benhige gt th) +(¢ +¢,) {a} =P 


Comparison with equations (1) shows that 


Ly=L,+Ls; Ly2=L2+Ls Dy,=La= —L3 
Ru=hit+hks Roe=Rho+Rs Ry= Ro = —R3 
it Araleg el ee! 1 Onl oe i 
Cu (iC, Cu Ca" Cs Ce Cu Cs 
It should be observed that the signs of the mutual coefficients 
Ry, Liz, Cig are a matter of convention. For example if the 
conventional directions of Jz and H, are reversed, the signs of 
the mutual coefficients are reversed. 

The system of equations (1) possesses two important properties 
which are largely responsible for the relative simplicity of classical 
electric circuit theory. First, the equations are linear in both 
currents and applied electromotive forces, Secondly, the 
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coefficients Lj, Rix, Cy, are constants. Important electro- 
technical problems exist, in which these properties no longer 
obtain. The solution, however, for the restricted system of 
linear equations with constant coefficients is fundamental and 
its solution can be extended to important problems involving 
non-linear relations and variable coefficients. These extensions 
will be taken up briefly in a later chapter. 


Er} o Li Le e JE2 


Ries 1: 


Another important property is the reciprocal relation among 
the coefficients; that is Lj;,=L,;, Rj,=Rx;, and C;,.=Cy;. It is 
easily shown that these reciprocal relations mean that there are 
no concealed sources or sinks of energy. Again important cases 
exist where the reciprocal relations do not hold. Such excep- 
tions, however, while of physical interest do not affect the mathe- 


Ri Ci ui Le C2 R2 

| R3 | 

Et ) C3 (i | Ez 

| : L3 | 
Fie. 2. 


matical methods of solution, to which the reciprocal relation is 
not essential. 

Returning to equation (1) we shall now derive the equation of 
activity. Multiply the first equation by J;, the second by J., 
etc. and add: we get 


Se sD I; Li +55 6. 95 %+ 
SSReih= Say. @ 
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The right hand side is the rate at which the applied forces are 
supplying energy to the network. The first term on the left is 
the rate of increase of the magnetic energy 


Ss Lit; Tn, 


while the second term is the rate of increase of the electric energy 


1 if 

2) iC 2 
The last term, DD ale is the rate at which electromagnetic 
energy is being converted into heat in the network. Consequently 
in the electrical network, the magnetic energy is a homogeneous 
quadratic function of the currents, the electric energy is a homo- 
geneous quadratic function of the charges, and the rate of 
dissipation is a homogeneous quadratic function of the currents. 
In Maxwell’s dynamical theory of electrical networks, these 
relations were written down at the start and the circuit equations 
then derived by an application of Lagrange’s dynamic equations 
to the homogeneous quadratic functions. 

Returning to equations (1), we observe that, due to the 
presence of the integral sign, they are integro-differential equa- 
tions. They are, however, at once reducible to differential 
equations by the substitution [=dQ/dt, whence they become 


d? ane d? d 
(Ln jat Rn pt Su)Qct: 2S (Ls jet Rng tSin)Qe= Es 


Pelee ae @ dice 
mi qt nt gg ts n1 Qit ra sp Linn qpt Run gy ts r»)Qa= En 


Here, as a matter of convenience, we have written 1/Cj,.=S jx. 
It is often more convenient, at least at the outset, to deal with 
equations (3) rather than (1). 

The Exponential Solution.—In taking up the mathematical 
solution of equations (1), we shall start with the exponential 
solution. This is of fundamental importance, both theoretically 
and practically. It serves as the most direct introduction to the 
Heaviside operational calculus, and in addition furnishes the 
basis of the steady state solution, or the theory of alternating 
currents. 

To derive this solution we set H,=F,e™ and put all the other 
forces He, . . H, equal to zero. This latter restriction is a 
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mere matter of convenience, and, in virtue of the linear character 
of the equations, involves no loss of generality. 

Now, corresponding to H,=F,e™, let us assume a solution of 
the form 


p= dee" (G=1,2 .. 7) 


where J; is a constant. So far this is a pure assumption, and 
its correctness must be verified by substitution in the differential 
equations. 

Now if J;=J,e, it follows at once that 


Sth =); =I; 


and 
le waa 
| Tjdi=Tj—sJ4 0%. 


Now substitute these relations in equations (1) and cancel the 
common factor e’. We then get the system of simultaneous 
equations 


OALwtRitl/rACiy)dJit ee + (ALint+Rintl/ACin) In =f, 
ALer+Roaytl1/rACo1) Jit can + (ALDon+Ron+1/AC on) In =0, 


OLnit+Rartl/ACni) Jit as ar Lan tRantl/ACnn) dn = (0) 


We note that this is a system of simultaneous algebraic equa- 
tions from which the time factor has disappeared. It is this that 
makes the exponential solution so simple, since we can immedi- 
ately pass from differential equations to algebraic equations. 
In these algebraic equations, n in number, there are » unknown 


quantities J;, . . J». These can therefore all be uniquely 
determined. We thus see that the assumed form of solution is 
possible. 


The notation of equations (4) may be profitably simplified 
as follows: write 
AL ik +R jz + il [DC ik = 2 ik (A) =ik 
and we have 
Gad g) ete ot +Zindn=fi, 
Zoid 1+2o2Jde+ +: +2endn=0, (5) 


2niditenadot::+éandn=0. 
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The solution of this system of equations is 


Mr, Mi 
eG aa ae. 
and 
Mi nM Fi M 
ee = 6 
I; D Fye Ta (6) 
where D is the determinant of the coefficients, 
Ain Caw) iliac ag OH ooo U Zin 
(DN LO) (No BAGO oe Zon 
FAME Oe cen oprsech Gamer fi. eo 23n (7) 
Peer bit hs copa o O.0 O10 enn 


and M;, is the cofactor, or minor with proper sign, of the jth 
column and first row. 

I shall not attempt to discuss the theory of determinants on 
which this solution is based.! We may note, however, one 
important property. Since 2;,=2.;, Mj,=M;;. From this the 
Reciprocal Theorem follows immediately. This may be stated 
as follows: 

If a force Fe“ is applied in the jth mesh, or branch, of the 
network, the current in the kth mesh, or branch, is by the 
foregoing 

Mire, 


Now apply the same force in the kth mesh, or branch, then the 
current in the jth mesh is 


Comparing these expressions and remembering that M;;=M,., 
it follows that the current in the kth branch corresponding to an 
exponential impressed e.m.f. in the jth branch, is equal to the 
current in the jth branch corresponding to the same e.m.f. in 
the kth branch. This relation is of the greatest technical 
importance. 

In many important technical problems we are interested only 
in two accessible branches, such as the sending and receiving. 
In such cases, where we are not concerned with the currents in 
the other meshes or branches, it is often convenient to eliminate 

1For a remarkably concise and complete discussion of the exponential 


solution by aid of the theory of determinants, see Cisoidal Oscillations, 
Trans. A. I, E. E., 1911, by G. A. Campbell. 
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them from the equation. Thus, suppose that we have electro- 
motive forces H; and EH, in meshes 1 and 2 and are concerned 
only with the currents in these meshes. If we solve equations 
3,4, . . m,n—2in number, for J; . . I, in terms of J, and I, 
and then substitute in (1) and (2) we get 

Zul;+Zi.I2=Fi, 

Zul tZool2=LE». (8) 


The Steady State Solution.—The steady state solution, on 
which the whole theory of alternating currents depends, is 
immediately derivable from the exponential solution. Let us 
suppose that 

E,=H;= + + » =E,=0 and that Hi=F cos (wt—@). 
Now by virtue of the well known formula in the theory of the 
complex variable, cos «=3e'*+3e—**, we can write 

By=4F eo -9 i pe e—9) 
=1(cos @—1 sin 0)Fe+4(cos 0+7 sin 0)Fe~*, (9) 
=P Lape tet 
Now, by virtue of this formula, the applied electromotive force E; 
consists of two exponential forces, one varying as e’” and the 
other as e *“. Hence it is easy to see that the currents are 
made up of two components, thus 


Demet A j=1,2) 7. 2) (10) 
and we have merely to use the exponential solution given above, 
substituting for \,7w and —7w respectively. That is, 


l= ia JH 1 Ais 
la ok Gee lees a. > ee 
Un 7 Gat 0 Ta) 
‘s Fe? eee he 
= 1 ge iwt Ree — tut 
SACS Uy AT ESR yma 
The second term is the conjugate imaginary of the first, so that 
ae tere tot 
EFS ARTI 
) ie 
ay ee _pt(wt — 0) 
ARC 
| ee 
=R aa 1(wt—0— p) 
Zilia) 
F 
~ Z(Gw)| cos (wt-—6— 9). 


We thus arrive at the rule for the steady state solution: 
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If the applied e.m.f. is 7 cos (wf—6), substitute 7 for d/dé in 
the differential equations, determine the impedance function 
Z (iw) =D (iw) /M (iw) (11) 
by the solution of the algebraic equations, and write it in the 
form 


Z (iw) =|Z (iw)|e’*. (12) 
Then the required pee is 
ae Zia) cos (wt—6— 9). (13) 


This in compact form contains the whole theory of the symbolic 
solution of alternating current problems. 

The Complementary Solution.—So far in the solutions which 
we have discussed the currents are of the same type as the 
impressed forces: that is to say in physical language, the currents 
are ‘“‘forced”’ currents and vary with time in precisely the same 
manner as do the electromotive forces. Such currents are, 
however, in general only part of the total currents. In addition 
to the forced currents we have also the characteristic oscillations; 
or, In mathematical language, the complete solution must include 
both particular and complementary solutions. This may be 
shown as follows: Let J)’, . . . . I,’ be solutions of the comple- 
mentary equations, 


d i ; d 1 on 
(Engthntg. | a) Lee (Lani thing. { aera =0, 


eo de: Neate eons ual Mamie ee eee (14) 
d 1 , i 

(Lng tate. it) Ty i : se (u Gt Rae to Ca iI =0. 

Then if I; 2. « J, 18 arsolution of (1), y-iy. . dati, cis 


also a solution. 
To derive the solution of the complementary system of equa- 
tions (14), assume that a solution exists of the form 
i= lee. G=12 .% %) 
so that d/dt=\ and f di=1/). Substitute in equations (14) 
and cancel out the common factor e“. Then we have 
211(A) Jj’ + a +2in(A) Jn’ =0, 
gece iou sie ara A ae ae eee (15) 
aniA)dr + «+++ tean(A)Jn’ =0. 
This is a system of n homogeneous equations in the unknown 
quantities J;’, . . Jn’. The condition that a finite solution 
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shall exist is that, in accordance with a well known principle of 
the theory of equations, the determinant of the coefficients 
shall vanish. That is, 


Pee cee eae if (16) 


Pal (Neen 2) | 


Consequently the possible values of \ must be such that this 
equation is satisfied. In other words, \ must be a root of the 
equation D(\)=0. Let these roots be denoted by di, \2, . - Am: 
Then, assigning to \ any one of these values, we can determine 
the ratio J;’/J;’ from any (n—1) of the equations. That is to 
say, if we take 


Patel ts gue aaa +Cne”, ) 
substitution in any (n—1) of the equations determines I,’, . . I,,’. 
The m constants Ci, . . C,,aresofar, however, entirely arbitrary, 


and are at our disposal to satisfy imposed boundary conditions. 

This introduces us to the idea of boundary conditions which 
is of the greatest importance in circuit theory. In physical 
language the boundary conditions denote the state of the system 
when the electromotive force is applied or when any change in 
the circuit constants occurs. The number of independent 
boundary conditions which can, in general, be satisfied is equal 
to the number of roots of the equation D(A)=0. Evidently, 
therefore, it is physically impossible to impose more boundary 
conditions than this. On the other hand, if this number of 
boundary conditions is not specified, the complete solution is 
indeterminate: that is to say, the problem is not correctly set. 
As an example of boundary conditions, we may specify that the 
electromotive force is applied at time ¢=0, and that at this time 
all the currents in the inductances and all the charges on the 
condensers are zero. 

So far we have been following the classical theory of lnear 
differential equations. We have seen that the forced exponential 
solution and the derived steady state solution are extremely 
simple and are mere matters of elementary algebra. The 
practical difficulties in the classical method of solution begin 
with the determination of the constants Ci, . . Cn of the 
complementary solution as well as the roots \1, . . Am of the 
equation D(\)=0. It is at this point that Heaviside broke 
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with classical methods, and by considering special boundary 
conditions of great physical importance, and particular types of 
impressed forces, laid the foundations of original and powerful 
methods of solution. We shall therefore at this point follow 
Heaviside’s example and attack the problem from a different 
standpoint. In doing this we shall not at once take up an 
exposition of Heaviside’s own method of attack. We shall 
first establish some fundamental theorems which are extremely 
powerful and will serve us as a guide in interpreting and ration- 
alizing the Heaviside operational calculus. 


CHAPTER II 


THE SOLUTION WHEN AN ARBITRARY FORCE IS 
APPLIED TO THE NETWORK IN A STATE OF 
EQUILIBRIUM 


In engineering applications of electric circuit theory there are 
four outstanding problems: 

(1) The steady state distribution of currents and potentials 
when the network is energized by a sinusoidal electromotive 
force. This problem is the subject of the theory of alternating 
currents which forms the basis of our calculations of power lines 
and the more elaborate networks of communication systems. 

(2) The distribution of currents and potentials in the network 
in response to an arbitrary electromotive force applied to the 
network in a state of equilibrium, i.e., applied when the currents 
and charges in the network are identically zero. 

(3) The effect on the distribution of currents and potentials of 
suddenly changing a circuit constant or connection, such as 
opening or closing a switch, while the system is energized. 

(4) The free oscillations of the network from an arbitary initial 
configuration. This problem will not be considered in detail 
for several reasons. In the first place it is of minor technical 
importance. Furthermore, if the initial state is producible by 
electromotive forces applied at accessible terminals of the net- 
work, the problem is reducible to and can be dealt with as a 
special case of problem (3). Finally, by reason of the arbitrary 
initial boundary conditions, the formulation and solution of the 
problem is inherently incapable of reduction to the elegant sim- 
plicity of those problems directly amenable to the operational 
calculus. The formal solution of problem (4) is, however, 
straightforward and is dealt with in treatises on dynamics. 
After evaluating the roots of the equation D(A)=0 (see last 
paragraph of Chap. I), the constants of integration C,, . . Cm 
of the complementary solution must be determined to satisfy 
the initial conditions or configuration. Heaviside, by aid of 
his conjugate Energy Theorem, has given an elegant and rela- 
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tively simple expression for these constants in terms of the initial 
currents in the inductances and the initial charges on the con- 
densers of the network.+ 

We shall base our further analysis of circuit theory on the solu- 
tion of problem (2), for the following reasons: 

(A) It is essentially a generalization of the Heaviside problem 
and its solution will furnish us a key to the correct understanding 
and interpretation of operational methods and lead to an auxiliary 
formula from which the rules of the operational calculus are 
directly deducible. 

(B) The solution of problem (2) carries with it the solution of 
problem (3) and also serves as a basis for the theory of alternating 
currents. 

(C) The solution of problem (2) leads directly to an extension of 
circuit theory to the case where the network contains variable ele- 
ments: i.e., circuit elements which vary with time and in which 
non-linear relations obtain. 

Problem (2) is therefore the fundamental problem of circuit 
theory and the formula which we shall now derive may be termed 
the fundamental formula of circuit theory. 

Consider a network in any branch of which, say branch 1, a 
unit e.m.f. is inserted at time ¢ = 0, the network having been 
previously in equilibrium. By unit e.m.f. is meant an electro- 
motive force which has the value unity for all positive values of 
time ({>0). Let the resultant current in any branch, say branch 
n, be denoted by Ani(t). Ani(é) will be termed the zndicial admit- 
tance of branch n with respect to branch 1—or, more fully, the 
transfer indicial admittance. 

The indicial admittance, aside from its direct physical signifi- 
cance, plays a fundamental réle in the mathematical theory of 
electric circuits. In words, it may be defined as follows: The 
indicial admittance, A,1(t), is equal to the ratio of the current 
in branch n, expressed as a time function, to the magnitude of 
the steady e.m.f. suddenly inserted at time {=0 in branch 1. 
It is evidently a function which is zero for negative values of 
time and approaches either zero or a steady value (the d.c. admit- 


‘ The reader who wishes to study this problem further is referred to the 
following papers: 

K. W. Wagner (E. N. T., Heft 11, Band 2, 1925). 

M.S. Vallarta (A. I. E. E., Apr., 1926, ‘““Heaviside’s Proof of his Expansion 
Theorem’’). 
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tance) for all actual dissipative systems, as ¢ approaches infinity. 
It may be noted that, aside from its mathematical determination, 
which will engage our attention later, it is an experimentally 
determinable function. 

We note, in passing, an important property of the indicial 
admittance A ;(é), which is deducible from the reciprocal theo- 
rem:! this is that A ;,(¢)=A;,;(¢). That is to say, the value of the 
transfer indicial admittance is unchanged by an interchange of 
the driving point and receiving point. It is therefore immaterial 
in the expression A ;,(¢) whether the e.m.f. is inserted in branch 
j and the current measured in branch k, or vice-versa. In 
general, unless we are concerned with particular branches, the 
subscripts will be omitted and we shall simply write A(t), it 
being understood that any two branches or a single branch (for 
the case of equal subscripts) may be under consideration. 

From the linear character of the network, it is evident that if a 
steady e.m.f. H=£, is inserted at time ¢=7, the network being in 
equilibrium, the resultant current is 


E, A(t—r). 


Generalizing still further, suppose that steady e.m.fs. Ko, Fi, 
E,, . . . E, are impressed in the same branch at the respective 
times 0, 71, 72 . . . Tn; the resultant current is evidently 


n 


ByAW)+EiAb—11)+--+E,A(t—7,)= SE At—7)). (18) 


7=0 


To apply the foregoing to our problem we suppose that there is 
applied to the network, initially in a state of equilibrium, an 
e.m.f. H(é) which has the following properties. 

1. It is identically zero for ¢<0. 

2. It has the value H(0) for 0<t<At. 

3. It has the value £(0)+A,# for At<t<2At. 

4, It has the value H(0)+A,E+A2E for 2At<t<3At. 

In other words it has the increment A;# at time ¢= At. 

Evidently then the resultant current J(¢) is 


E,A(t)+A1 BA (t—At)+-++A,E.A(t—nAt). 


1 Exceptions to this relation exist where the network contains sources of 
energy such as amplifiers. These need not engage our attention here. 
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Now evidently if the interval Aé is made shorter and shorter, 
then in the limit Aft-dé and jAt = 7 and 


d 
A;E = ph ar. 


Passing to the limit in the usual manner this summation becomes 
a definite integral and we get 


I(t) =E(0)A()+ |, A(t—1)o E(o)dr. (19) 
0 


Finally by obvious transformations of the expression we arrive 
at the fundamental formula of circuit theory! 


1) =4 il ‘A(t—1)E(@)dr, (20) 
d t 
=$ |, BGS aAaE (20a) 


For completeness we write down the following equivalents? of 
(20) and (20a) 


I(t)=A(0)E(@) + | ‘Alt—d Ear, (20b) 
=A(0)E()+ i; ‘A ()Et—1dr, (20e) 
=B(0)A(t)+ || Bt—n)A(dr, (20d) 
=B(0)A()+ | 2 @s—nur. (20e) 


where the primes denote differentiation with respect to the 
argument. Thus A’(t)=d/dtA (é). 

These equations are the fundamental formulas which mathe- 
matically relate the current to the type of applied electromotive 
force and the constants and connections of the system, and 


1 This important theorem, independently derived and published by the 
author, is actually the equivalent of a much older theorem in dynamics due 
to Duhamel. Substantially the same proof of the theorem as given in the 
text was independently communicated to the author by H. W. Nichols and 
by Stuart Ballantine. 

*In case either A(t) or H(t) has discontinuities for positive values of ¢, 
we integrate through these discontinuities in the usual manner. When this 
is borne in mind, equations (20b) to (20e) are not restricted to continuous 
functions. 


SOLUTIONS WHEN AN ARBITRARY FORCE IS APPLIED 17 


constitute the first part of the solution of our problem. The 
most important immediate deductions from these formulas are 
expressed in the following theorems. 

1. The indicial admittance of an electrical network completely 
determines, within a single integration, the behavior of the net- 
work to all types of applied electromotive forces. Ags a corollary, 
a knowledge of the indicial admittance is the sole information 
necessary to completely predict the performance and character- 
istics of the system, including the steady state. 

2. The applied e.m.f. and the indicial admittance are similarly 
and coequally related to the resultant current in the network. 
As a corollary the form of the current may be modified either by 
changing the constants and connections of the network or by 
modifying the form of the applied e.m.f. 

3. Since the applied e.m.f. may be discontinuous these formulas 
determine not only the building up of the current in response to 
an applied e.m.f. but also its subsidence to equilibrium when the 
e.m.f. is removed and the network left to itself. In brief, formu- 
las (20) reduce the whole problem to a determination of the 
indicial admittance of the network. In addition, as we shall 
see, they lead directly to an integral equation which determines 
this function. 

It is of interest to show the relation between formulas (20) 
and the usual steady state equations. To do this let the e.mf., 
applied at time t=0, be EF sin (wt+6). Substitution in for- 
mula (20c) and rearrangement gives 


I(t) =A(0)E sin (wt+@) 
+E sin (ot+0) [cos wrA!(r)dr 
0 


—E cos (ot-+0) [sin wrA’(r)dr (21) 
0 
where A'(i)= u A(t). 
Now this can be resolved into two parts 


E sin (wt+4) [4 (0)+ ] cos wrA’(r)dr | 
a (22) 


—E cos (wf+8) { { sin or’ e)dr| 
0 
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which is the final steady state, and , 


—FE sin (aio) | cos wrA’(r)dr 
oe (23) 

+E cos (oto) | sin wrA’(r)dr 
t 


which is the transient distortion, which ultimately dies away for 
sufficiently large values of time. 

To correlate the foregoing expressions for the steady state with 
the usual formulas we observe that if the symbolic impedance of 
the network at frequency w/2z be denoted by Z(iw), and if we 
write 


it 
Tans 
then the steady state current is 


E [a(w) . sin (wi +0)+6(w) . cos (wi + 4)]. 


Comparison with (22) gives at once 
a(w) =A (0) +f cos wrA’(r)dr, (24) 
0 


Bio) = -{ sin wrA’(r)dr. (25) 
0 


The Integral Equation for the Indicial Admittance.—So far 
we have tacitly assumed that the indicial admittance is known. 
As a matter of fact its determination constitutes the essential 
part of our problem. It is, in fact, the Heaviside problem, and its 
investigation, to which we now proceed, will lead us directly to 
the operational calculus. 

Heaviside’s method in investigating this problem was intuitive 
and ‘‘experimental.’”’ We, however, shall establish a very gen- 
eral integral equation from which we shall directly deduce his 
methods and extensions thereof. 

Let us suppose that an e.m.f. e?’, where p is either a positive real 
quantity or complex with real part positive, is suddenly impressed 
on the network at timet=0. It follows from the foregoing theory 
that the resultant current /(t) will be made up of two parts, (1) a 
forced exponential part which varies with time as e?‘, and (2) a 
complementary part which we shall denote by y(t). The expo- 
nential or “forced”? component is simply e?'/Z(p), where Z(p) 
is functionally of the same form as the usual symbolic or complex 
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impedance Z(iw). It is gotten from the differential equations 
of the problem, as explained in a preceding chapter, by replacing 
d"/dt" by p”, cancelling out the common factor e?, and solving 
the resulting algebraic equations. The complementary or char- 
acteristic component, denoted by y(t), depends on the constants 
and connections of the network, and on the value of p. It does 
not, however, contain the factor e?' and it dies away for suffi- 
ciently large values of ¢, in all actual dissipative systems. Thus 


ert 
I) =a Ue 26 
O=z—) tv (26) 
Now return to formula (20a) and replace E(t) by e?'. We get 
t 
T(t) =“ en il A(s)e-”* dr 
dt ‘ 


which can be written as 


a | al “A (r)e—?" dr —e** [4 (r)e—?? dr \ F 
dt ‘ : | 


Carrying out the indicated differentiation this becomes 
I(é) =pert{ A(r)e-?" r— per | A(rje**dr+A(t). (27) 
0 t 


Equating the two expressions (26) and (27) for I(t) and dividing 
ee by e?! we get 


+y (te np. Aire dr—p | Weed LAe (28) 
aa ) t 


This equation is valid for all values of t. Consequently if we set 
t= and if the real part of p- is positive, only the first term on 
the right and the left hand side of the equation remain, the rest 
vanishes, and we get 


1 i < 
——_ = A (t)e?' dt. (29) 
pZ(p) 0 

This is an integral equation' valid for all positive real values of p, 
which completely determines the indicial admittance A(t). It is 


1 An integral equation is one in which the unknown function appears under 
the sign of integration. (29) is an integral equation of the Laplace type. If 
Z(p) is specified, A(t) is uniquely determined. Methods for solving the 
integral equations are considered in detail later, in connection with the 
exposition of the operational calculus. The phrase ‘‘all positive values of p” 
will be understood as meaning all values of p in the right half of the complex 
plane. 
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on this equation that we shall base our discussion of operational 
methods and from which we shall derive the rules of the oper- 
ational calculus. Equations (20) and (29) constitute a complete 
mathematical formulation of our problem, and from them the 
complete solution is obtainable without further recourse to the 
differential equations, or further consideration of boundary 
conditions. 


CHAPTER III 


THE HEAVISIDE PROBLEM AND THE OPERATIONAL 
EQUATION 


By virtue of the analysis and formulas of the preceding 
chapters, we are now in a position to take up a systematic 
exposition and development of the Heaviside operational calculus. 
This analysis may be reviewed and summarized as follows: 

The determination of the current in a network in response to 
an electromotive force H(¢), impressed on the network at reference 
time ¢=0 has been reduced to the mathematical solution of two 
equations: first the integral equation 


il io) 
sn = | AW@etat 29 
ej || Ae (29) 
and second, the definite integral 
t 
nes A (t—1)E(r)dr. (20) 
dt Jo 


In equation (29), Z(p) is the impedance function of the network 
and the characteristics, i.e., constants and connections, of the 
network enter into the problem only through Z(p), and, therefore, 
Z(p) completely specifies the network. If p is replaced by iw, 
where w/2z is the frequency, then Z(iw) is the complex steady 
state or symbolic impedance of alternating current theory. 
Thus, if Z(p) is specified, equations (20) and (29) constitute a 
complete mathematical formulation of the problem, and their 
solution automatically takes care of the equilibrium boundary 
conditions. That is to say, contrary to the classical method of 
solution, the explicit determination of integration constants is 
completely dispensed with. 

It will be observed that in deducing these equations we have 
merely postulated (1) the linear and invariable character of the 
network and (2) the existence of an exponential solution of the 
type e?‘/Z(p) for positive values of p. Consequently, while we 
have so far discussed these formulas in terms of the determination 
of the current in a finite network, they are not limited in their 
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application to this specfic problem. In this connection it may 
be well to call attention explicitly to the following points. 

The formulas and methods deduced above apply not only to 
finite networks, involving a finite system of linear equations, but 
to infinite networks and to transmission lines, involving infinite 
systems of equations, and partial differential equations: in fact 
to all electrical and dynamical systems in which the connections 
and constants are linear and invariable. 

Secondly, the variable determined by formulas (20) and (29) 
need not, of course, be the current. It may equally well be the 
charge, potential drop, or any of the variables with which we 
may happen to be concerned. ‘This fact may be explicitly recog- 
nized by writing the formulas as: 


Vis at 
pH (p) = {| h(t)e-”! dt, (30) 
nt) =% {| N= neGlee (31) 


Here E(t) is the applied e.m.f., x(¢) is the variable which we desire 
to determine (charge, current, potential drop, etc.), and 

a= E/H(p) (32) 
is the operational equation. H(p) therefore corresponds to and is 
determined in precisely the same way as the impedance Z(p), but 
it may not have the physical significance or the dimensions of an 
impedance. Similarly, in character and function, h(t) corre- 
sponds to the indicial admittance, though it may not have the 
same physical significance. It is a generalization of the indicial 
admittance and may be appropriately termed the Heaviside 
Function. Similarly H(p) may be termed the generalized 
impedance function. 

The Operational Equation.—The physical problem which 
Heaviside attacked and which led to his operational ce!eculus was 
the determination of the response of a network or electrical sys- 
tem to a “unit e.m.f.”’ (zero before, unity after time t=0) with, 
of course, the understanding that the system is in equilibrium 
when the electromotive force is applied. His problem is therefore, 
essentially that of the determination of the indicial admittance. 
In our exposition and critique of Heaviside’s method of dealing 
with this problem we shall accompany an account of his own 
method of solution with a parallel solution from the corresponding 
integral equation of the problem. 
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Heaviside’s first step in attacking this problem was to start 
with the differential equations, and replace the differential opera- 
tor d/dt by the symbol p, and the operation {dt by 1/p, thus 
reducing the equations to an algebraic form. He then wrote the 
impressed e.m.f. as 1 (unity), thus limiting the validity of the 
equations to values of f>0. The formal solution of the algebraic 
equations is straightforward and will be written as 

h=1/H(p) (33) 


) 


where h/ is the “generalized indicial admittance,’ or Heaviside 
function (denoting current, charge, potential or any variable 
with which we are concerned) and H(p) is the corresponding 
generalized impedance. Thus, if we are concerned with the 
current in any part of the network, we write 

A=1/Z(p). (34) 
The more general notation is desirable, however, as indicating 
the wider applicability of the equation. 

The equations 

h=1/H(p) 

A=1/Z(p) 
are the Heaviside operational equations. ‘They are, as yet, purely 
symbolic and we have still the problem of determining their 
explicit meaning and in particular the significance of the operator 
D. 

Comparison of the Heaviside operational equations with the 
integral equations (29) and (30) of the preceding chapter leads to 
the following fundamental theorem. 

The Heaviside operational equations 

A=1/Z(p) 

h=1/H(p) 
are merely the symbolic or short-hand equivalents of the correspond- 
ing integral equations 


le ae 
sey |, A(the-P! dt 


i a ies ey 


The integral equations, therefore, supply us with the meaning and 
significance of the operational equations, and from them the rules 
of the operational calculus are deducible. 
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By virtue of this theorem, we have the advantage, at the 
outset, of a key to the meaning of Heaviside’s operational 
equations, and a means of checking and deducing his rules of 
solution. This will serve us as a guide throughout our further 
study. 

Returning now to Heaviside’s own point of view and method of 
attack, his reasoning may be described somewhat as follows:— 
The operational equation 

h=1/H(p) 

‘is the full equivalent of the differential equations of the problem 
and must therefore contain the information necessary to the 
solution provided we can determine the significance of the sym- 
bolic operator p. The only way of doing this, when starting 
with the operational equation, is one of induction: that is, we 
must compare the operational equation with known solutions of 
specific problems and thus attempt to infer by induction general 
rules for interpreting the operational equation and converting 
it into the required explicit solution. 

The Power Series Solution 

Let us start with the simplest possible problem: the current in 
response to a “unit e.m.f.”’ in a circuit consisting of an inductance 
L in series with a resistance R. 

The differential equation of the problem is 


TOA +RA=1, #>0, 
dt 
where A is the indicial admittance. Consequently replacing 
d/dt by p, the operational equation is 


1 
es fy 


The explicit solution is easily derived: it is 
1 
=— —=$"_ p a 
A pu ao) 
where a=R/L. Note that this makes the current initially zero, 
so that the equilibrium boundary condition at t=0 is satisfied. 


Now suppose that we expand the operational equation in 
inverse powers of p: we get, formally, 


ik 1 la 1 i) la a\? a\s e\- ) 
A=— <— a a ey ee ee oe -_— ~ — aor 4! 
pLi+a/p Rpl+a/p R \p 6 a) © zs | 


by the Binomial Theorem. 
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Now expand the explicit solution as a power series in f: it is 


at [at_(@d* (ab? | 

Rill 2! 3! 
Comparing the two expansions we see at once that the operational 
expansion is converted into the explicit solution by assigning 
to the symbol 1/p” the value t”/n!. It was from this kind of 
inductive inference that Heaviside arrived at his power series 
solution. 

Now there are several important features in the foregoing which 
require comment. In the first place the operational equation is 
converted into the explicit solution only by a particular kind of 
expansion, namely an expansion in inverse powers of the operator 
p. For example, if in the operational equation 


ot 


_1_a/p_ 
~ Rit+a/p 
we replace 1/p by t/1! we get 
1 at 
om iat 


which is incorrect. Furthermore, if we expand in ascending 
instead of descending powers of p, namely 


A= 5 {1-(p/a)-+(p/a)@—- ++ 


no correlation with the explicit solution is possible and no signifi- 
cance can be attached to the expansion. We thus infer the gen- 
eral principle, and we shall find this inference to be correct, 
that the operational equation is convertible into the explicit 
solution only by the proper choice of expansion of the impedance 
function, or rather its reciprocal. 

In the second place we notice that in writing down the opera- 
tional equation and then converting it into the explicit solution 
no consideration has been given to the question of boundary 
conditions. This is one of the great advantages of the opera- 
tional method: the boundary conditions, provided they are those 
of equilibrium, are automatically taken care of. This will be 
illustrated in the next example: 

Let a “unit e.m.f.”’ be impressed on a circuit consisting of 
resistance R, inductance L, and capacity C: required the resultant 
charge on the condenser. 


26 ELECTRIC CIRCUIT THEORY 


The differential equation for the charge Q is 


jets SVEN are sp et 
ae a 


Consequently the operational formula is 
1 


OT +Rp+1/C 
eet : wh Beni 
“Lp? 1+a/p+b/p? re eT LC 


This can be expanded by the Binomial Theorem as 


1 Geno Ga \a ee DN 
O=rps I-E+) ++) ~GGtes) +} 


Performing the indicated operations and collecting in inverse 
powers of 7p, the first few terms of the expansion are:— 


OG) FOG Ge ae 
where C=a 
(=b-a? 


C3 = 2ab—a® 

c4= b?—3a?b+a! 

C3 = 3ab?—4a3b-+a5 

cs = b? — 6a*b?+5a4 b—a 
We infer therefore that in accordance with the rule of replacing 
1/p” by t”/n! the solution is:— 

La? ie i ia i 

Om pla erg ceart erg tg) |. 
Owing to the complicated character of the coefficients in the 
expansion, the series cannot be recognized and summed by inspec- 
tion. If, however, we put R=0 then a=0, and theseries becomes 


“legal sieve taiwan) - | 


Q=C{1—cos (t/+/LC)}. 
We have still to verify this solution by comparison with the 
explicit solution of the differential equation. This is of the form 


Q = C+k.e+-h.e 


whence 
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where k; and kz are constants which must be chosen to satisfy the 
boundary conditions and i, d»2 are the roots of the equation 


L\2+Rd+1/C =0. 


Now since we have two arbitrary constants we satisfy the 
equilibrium conditions by making Q and dQ/dt zero at t=0. 
whence 
C+k, +k. = 0, 
Arkitnys ko SS 0, 
and 
ky =)2C/(1—da), 
ke = Ai C/(r2 — 1). 
We have also 


__4 o\o™ 
2 


nS) 


Writing down the power series expansion of 


Q —_ C+kie'+kee, 
then 


Q=(C+ hit ke) +(e dtha ds) 4 


2 
+ (ade hehe!) + NAD 


Introducing the values of ky, kz, \1, 42 given above and comparing 
with the power series derived from the operational solution we see 
that they are identical term by term. 

This example illustrates two facts. First the power series 
expansions may be complicated, laborious to derive and of such 
form that they cannot be recognized and summed by inspection. 
In fact in arbitrary networks of a large number of meshes or 
degrees of freedom the evaluation of the coefficients of the power 
series expansion is extremely laborious. 

On the other hand, in such cases, the solution by the classical 
method presents difficulties far more formidable—in fact insuper- 
able difficulties from a practical standpoint. First there is the 
location of the roots of the function H(\), which in arbitrary 
networks is a practical impossibility without a prohibitive amount 
of labor. Secondly there is the determination of the integration 
constants to satisfy the imposed boundary conditions: a process, 
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which, while theoretically straightforward, is actually in practice 
extremely laborious and complicated. We note these points in 
passing; a more complete estimate of the value of the power 
series solution will be made later. 

To summarize the preceding: Heaviside, generalizing from 
specific examples otherwise solvable, arrived at the following 
rule :— 

Expand the right hand side of the operational equation 

h=1/H (p) 
in inverse powers of p: thus 


h—~Go-hai/ P--Ge/ Pp 4 On) p"-P 
1 : aan 
and then replace by @/n!. The operational equation is thereby 


converted into the explicit power series solution :— 
h=aotayt/1!+aet?/2!4+--+a,i*/n!4+--- (35) 
As stated above, this rule was arrived at by pure induction and 
generalization from the known solution of specific problems. It 
cannot, therefore, theoretically be regarded as satisfactorily 
established. The rule can, however, be directly deduced from 


the integral equation 
il ~*~ 
= h(t)e—”! dt. 
pH (p) | iy 


To its derivation from this equation we shall now proceed. 
First suppose we assume that h(t) admits of the power series 
expansion 
het -Pa til lt hgt?/2\er eee 
Substitute this assumed expansion in the integral, and integrate 
term by term. The right hand side of the integral equation 
becomes formally 


ho/phi/p?+he/p +: +++: 
by virtue of the formula 


aie ih 
os ee 
{| ni’ pee p>o. 


Now expand the left hand side of the integral equation asymp- 

totically in inverse process of p: it becomes 
ao/p+ai/p?+a2/pi+-+--: 

where 


Ao +ai/p+a2/p?+ - 38 
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is the asymptotic expansion of 1/H(p). Comparing the two 
expansions and making a term by term identification, we see that 
h,=G, and 

h(t) =ao+ayt/1!+agt?/2!4+---- 
which agrees with the Heaviside formula. 

This procedure, however, while giving the correct result has 
serious defects from a mathematical point of view. For example, 
the asymptotic expansion of 1/H(p) has usually only a limited 
region of convergence, and it is only in this region that term by 
term integration is legitimate. Furthermore we have assumed 
the possibility of expanding h(t) in a power series: an assumption 
to which there are serious theoretical objections, and which, 
furthermore, is not always justified. A more satisfactory deriva- 
tion, and one which establishes the condition for the existence of a 
power series expansion, proceeds as follows :— 

Let 1/H(p) be a function which admits of the formal asymp- 
totic expansion 


> n/p” 
0 | 
and let it include no component which is asymptotically repre- 
sentable by a series all of whose terms are zero, that is a function 
$(p) such that the limit, as p— 0, of p” $(p) is zero for every value 
of n. Such a function is e~?. With this restriction understood, 
start with the integral equation, and integrate by parts: we get 


1 a a —pt 1 
Ag O+ | oP A (A) dt 


where h™(t) denotes d"/di" h(t). Now let p approach infinity: 
in the limit the integral vanishes and by virtue of the asymptotic 
expansion 


1/H(p)~ >, on/p", (36) 
| 10] 
1/H(p) approaches the limit a. Consequently 
h(O) =a. 
Now integrate again by parts: we get 


p(1/H (p) —a) =h (0) + tf “ePH (dl, 
0 
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Again let p approach infinity: in the limit the left hand side of 
the equation becomes a; and we have 


h® (0) =a. 


Proceeding by successive partial integrations we thus establish 


the general relation 
h™ (0) = an. 


But by Taylor’s theorem, the power series expansion of h(t) is 
simply 
h(t) =h(0) +h (0)t/1!+h® (0)82/214------ 


whence, assuming the convergence of this expansion, we get 
h(t) = dg tayt/1!+a2t?/2!4+ --- = S\antr/n! (85) 
0 


which establishes the power series solution. It should be 
carefully noted, however, that it does not establish the con- 
vergence of the power series solution. As a matter of fact, 
however, I know of no physical problem in which H(p) satisfies 
the conditions for an asymptotic expansion, where the power 
series solution is not convergent. On the other hand many 
physical problems exist, including those relating to transmission 
lines, where a power series solution is not derivable and does not 
exist. 

The process of expanding the operational equation in such a 
form as to permit of its being converted into the explicit solution 
is what Heaviside calls “algebrizing”’ the equation. In the 
case of the power series solution the process of algebrizing 
consists in expanding the reciprocal of the impedance function 
in an asymptotic series,! thus 


1/H(p)~do-+ai/p+a2/p? +++ 


Regarded as an expansion in the variable p, instead of as a 
purely symbolic expansion, this series has usually only a limited 
region of convergence. This fact need not bother us, however, as 
the series we are really concerned with is 


ao+ayt/1 I++a t?2/21+ : 


1The symbol ~™, frequently employed in the following pages, is to be 
understood as signifying “‘is asymptotically equal to.” In other words the 
two sides of the equation are not, in general, identities but are asymptotically 
equivalent or equal, 
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It is interesting to note in passing that the latter series is what 
Borel, the French mathematician, calls the associated function 
of the former, and is extensively employed by him in his 
researches on the summability of divergent series. 

The process of ‘‘algebrizing,” as in the examples discussed 
above, may often be effected by a straightforward binomial 
expansion. In other cases the form of the generalized impedance 
function H(p) will indicate by inspection the appropriate pro- 
cedure. A general process, applicable in all cases where a power 
series exists, is as follows. Write 


1/H(p) =1 (a(*) =G(2). (36) 
Now expand G(x) as a Taylor’s series: thus formally 
G(x) =G(0) +60) 4G CO)E +. 
where 


G(0)= | £62) | (37) 


=0. 


Denote 


(n) 
a ee by an, replace x” by 1/p”, and we have 


G(x) = 1/H(p) = aot+ai/p+ac/p?+-::- 

This process of “algebrizing”’ is formally straightforward and 
always possible. As implied above, however, in many problems 
much shorter modes of expansion suggest themselves from the 
form of the function H(p). 

We note here, in passing, that the necessary and sufficient 
conditions for the existence of a power series solution is the 
possibility of the formal expansion of ((a) as a power series in a. 

At this point a brief critical estimate of the scope and value of 
the power series solution may be in order. As stated above, in a 
certain important class of problems relating to transmission 
lines, a power series does not exist, though a closely related 
series in fractional powers of ¢ may often be derived. Con- 
sequently the power series solution is of restricted applicability. 
Where, however, a power series does exist, in directness and 
simplicity of derivation it is superior to any other form of solu- 
tion. Its chief defect, and a very serious defect indeed, is that 
except where the power series can be recognized and summed, 
it is usually practically useless for computation and interpreta- 
tion except for relatively small values of the time ¢. This 
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disadvantage is inherent and attaches to all power series solu- 
tions. For this reason I think Heaviside overestimated the 
value of power series as practical or working solutions, and that 
some of his strictures against orthodox mathematicians and 
their solutions may be justly urged against the power series 
solution. He was quite right in insisting that a solution must 
be capable of either interpretation or computation and quite 
right in ridiculing those formal solutions which actually conceal 
rather than reveal the significance of the original differential 
equations of the problem. On the other hand, the following 
remark of his indicates to me that Heaviside had a quite exagger- 
ated idea of the value and fundamental character of power 
series in general: “I regret that the result should be so compli- 
cated. But the only alternatives are other equivalent infinite 
series, or else a definite integral which is of no use until it is 
evaluated, when the result must be the series (135), or an equiv- 
alent one.” As a matter of fact the properties of most of the 
important functions of mathematical physics have been investi- 
gated and their values computed by methods other than series 
expansions. I may add that in technical work the power series 
solution has proved to be of restricted utility, while definite 
integrals, which Heaviside! particularly despised, have proved 
quite useful. 


The Expansion Theorem Solution? 

We pass now to the consideration of another extremely 
important form of solution. Heaviside gives this solution with- 
out proof: we shall therefore merely state the solution and then 
derive it from the integral equation. 

Given the operational equation 

h=1/H(p) 
which has the significance discussed above: i.e., the response of 
the network to a “unit e.m.f.” The explicit solution may be 
written as 


n 
1 eri! 


Ce H(0) a ed Dy, H" (x) 
1 


(38) 


' Vide a remark of his to the effect that some mathematicians took refuge 
in a definite integral and called that a solution. 

* This terminology is due to Heaviside. A more appropriate and physi- 
cally significant expression would be ‘‘The Solution in terms of normal or 
characteristic vibrations,” 
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where pi, p2 . . . pn» are the n roots of the equation 
H(p)=0 
and 
H"(p.) =| ple | : (39) 
dp D = Dk 


As remarked above, this solution, referred to by him as The 
Expansion Theorem, was stated by Heaviside without proof; 
how he arrived at it will probably always remain a matter of 
conjecture.! Its derivation from the integral equation is, 
however, a relatively simple matter, though in special cases 
troublesome questions arise. 

The derivation of the expansion solution from the integral 


equation 
if io0) 
== | A(ierdt 
pH (p) |, ee 


follows immediately from the partial fraction expansion 


n 


1 i i 
= 40 
pH(p) pH(0) a TEED a 
= 
where pi, po . . . Pn are the roots of the equation H(p)=0, 
and 
’ d 
We)=| FH@)| (41) 
dp P=Di 


Partial fraction expansions of this type are fully discussed in 
treatises on algebra and the calculus and the conditions for their 
existence established. Before discussing the restrictions imposed 
on H(p) by this expansion, we shall first, assuming its existence, 
derive the expansion theorem solution. 

By virtue of (40) the integral equation is 


1 - 1 pt, 
pHO)* 2@=peH@)™ ii herd. (42) 


The expansion on the left hand side suggests a corresponding 
expansion on the right hand side: that is, we suppose that 


h(t) =ho(t) +hilt) +he(t) ++ + +ha(t) (43) 


1Since this was written it has been independently shown by K. W. 
Wagner and M. 8, Vallarta (see footnote to p. 14) that a proof based on his 
“conjugate theorem” can be reconstructed from scattered papers of 
Heaviside. 
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and specify that these component functions shall satisfy the 
equations 
1 o 
a tice (44) 
1 eS : 
on h; (t)e-?* dt = oes (45) 
Cae i@ d 
It follows at once from (43) and direct addition of equations (44) 


and (45) that (42) is satisfied and hence is solved provided 
ho, . . hn can be evaluated from (44) and (45). 


Now since 
“oer dt= + (46) 
ee =—_ 
[pete 


provided the real part of \ is not positive (a condition satisfied 
in all network problems), we see at once that equations (42) 
and (48) are satisfied by taking 


1 
ho(t) =ho= H(0) (47) 
Fee = 02 
()=— =o Wee 
BES SE 
Consequently from (43) and (47) it follows that 
1 - epit 
h(t) == Ba ee an 48 
= HO) apt’) Sr 


which establishes the expansion theorem solution. 

As implied above, the partial fraction expansion (40), on 
which the expansion theorem solution depends, imposes certain 
restrictions on the impedance function H(p). Among these are 
that H(p) must have no zero root, no repeated roots, and 1/H(p) 
must be a proper fraction. In all finite networks these conditions 
are satisfied, or by a slight modification, the operational equation 
can be reduced to the required form. The case of repeated 
roots, which may occur where the network involves a unilateral 
source of energy such as an amplifier, can be dealt with by 
assuming unequal roots and then letting the roots approach 
equality as a limit. Without entering upon these questions in 
detail, however, we can very simply and directly establish the 
proposition that the expansion theorem gives the solution 
whenever a solution in terms of normal or characteristic vibrations 
exists. The proof of this proposition proceeds as follows. 
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It is known from the elementary theory of linear differential 
equations that the general solution of the set of differential 
equations, of which the operational equation is h=1/H(p), is of 
the form 


Opereres SS C,er# 
1 


where p; is the jth root of H(p)=0, and Cy, Ci . . . Cn are 
constants of integration which must be so chosen as to satisfy 
the system of differential equations and the imposed boundary 
conditions. The summation is extended over all the roots of 
H(p) which is supposed not to have a zero root or repeated roots. 

Now substitute this known form of solution in the integral 
equation of the problem and carry out the integration term by 
term. We get 


=C eae 49 
Hi) +7 Dp Di co 
Setting p=0, we have at once 

Co=1/H(0). (50) 


To determine C; let p=p;+gq where q is a small quantity ulti- 
mately to be set equal to zero, and write the equation as 


CoH (p + SEa, =f (51) 


If now p=p;+q and gq approaches zero, this becomes in the 
limit 


pi H'(p;)C;=1 (52) 
or 
1 
oa Hey Oke 
whence h(t) = ath = a — a (54) 


which is the Expansion Theorem Solution. 

We shall not attempt to discuss here cases where the expansion 
solution breaks down though such cases exist. In every such 
case, however, the breakdown is due to the failure of the imped- 
ance function H(p) to satisfy the conditions necessary for the 
partial fraction expansion (40), and correlatively the non- 
existence of a solution in normal vibrations. Furthermore, it is 
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usually possible by simple modifications to deduce a modified 
expansion solution. It may be added here, that while the 
proof given above is also limited implicitly to finite networks, 
the expansion solution is valid in most transmission line problems. 

Let us now illustrate how the expansion solution works by 
applying it to a few simple examples. Take first the case 
considered in the preceding chapter in connection with the 
power series solution. Required the charge Q on a condenser 
C in series with an inductance L and resistance R in rseponse 
to a “unit e.m.f.’”? The operational equation is 


o il 
~ Lp?+Rp+1/C 
or app een! 
Lp?+2ap+e? 


where a = R/2L and w? = 1/LC. 
The roots of the equation H(p) =0 are the roots of the equation 
p*+2ap+tw?= 
whence 
p= —atVa?—w= —at8, 
Do= —a—-Va?—w?= —a—B. 
Also H'(p) =2L(p+a), so that 
H'(p1) =26L 
H"(p2) = —26L 


and 
1/0) = 1G. 
Inserting these expressions in the Expansion Theorem Solution 
(38), we get 
—at Bt eft 
ay jas ( SS rctatie 
Be oar oot mare 
It is now easy to verify the fact that this solution satisfies the 
differential equations and the boundary condition 
Q=0 and dQ/dt=0 at time t=0. 
If w>a, 8 is a pure imaginary 


B=twV/1—(a/w)? = tw! 


XN 


and 
einen w’cos wt-+a sin w’t 
w’L ata!” 
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In connection with this problem we note two advantages of the 
expansion solution, as compared with the power series solution: (1) 
it is much simpler to derive from the operational equation, and (2) 
its numerical computation is enormously easier. A table of expo- 
nential and trigometric functions enables us to evaluate Q for any 
value of t almost at once whereas in the case of the power series 
solution the labor of computation for large values of tis very great. 
A third and very important advantage of the expansion solution 
in this particular problem is that without detailed computation 
we can deduce by mere inspection the general character of the 
function and the effect of the circuit parameters on its form: 
an advantage which never attaches to the power series solution. 

This last property of the particular solution above is extremely 
important. The ideal form of solution, particularly in technical 
problems, is one which permits us to infer the general character 
and properties of the function and the effect of the circuit con- 
stants on its form, without detailed computation. A solution 
which possesses these properties, even if its exact computation is 
not possible without prohibitive labor, is far superior to a solution 
which, while completely computable, tells us nothing without 
detailed computation. It is for this reason that some of the 
derived forms of solution, discussed later, are of such importance. 
In fact a solution which requires detailed computation before it 
yields the information implied in it is merely equivalent to an 
experimentally determined solution. 

Unfortunately the advantages attaching to the expansion 
solution of the specific problem just discussed, do not, in general, 
characterize the expansion solution. The following disadvan- 
tages should be noted. First, the location of the roots of the 
impedance function H(p) is practically impossible in the case of 
arbitrary networks of more than a few degrees of freedom. In 
the second place, when the number of degrees of freedom is large 
it is not only impossible to deduce the significance of the solution 
by inspection, but the computation becomes extremely laborious. 
In such cases, the practical value of the expansion solution 
depends, just as in the power series solution, on the possibility of 
recognizing and summing the expansion. This will be clear 
in the case of transmission lines, where the roots of H(p) are 
infinite in number and the direct computation of the expansion 
solution (except in the case of the non-inductive cable) is quite 


impossible. 


CHAPTER IV 


SOME GENERAL FORMULAS AND THEOREMS FOR THE 
SOLUTION OF OPERATIONAL EQUATIONS 


We have seen that the operational equation 
h=1/H(p) 
is the symbolic or short-hand equivalent of the integral equation 


ff oo a 


and from the latter we have deduced two very important forms of 
the Heaviside solution. In recognizing the equivalence of these 
two equations we have a very great advantage and are able, in 
fact, to base the operational calculus on deductive instead of 
inductive reasoning. In this chapter we shall employ this equiy- 
alence to establish certain general formulas and theorems for 
the solution of operational equations. That is to say, we shall 
make use of the principles that (1) any method applicable to 
the solution of the integral equation supplies us with a corre- 
sponding method for the solution of the operational equation, 
and (2) a solution of any specific integral equation gives at once 
the solution of the corresponding operational equation. We 
turn therefore to a brief discussion of the appropriate methods 
for solving the integral equation. 

It has been noted above that any solution of the integral equa- 
tion supplies a solution of the corresponding operational equation. 
This principle enables us to take advantage of the fact that a very 
large number of infinite integrals of the type 


[ soem 
0 


have been evaluated. The evaluation of every infinite integral of 
this type supplies us, therefore, with the solution of an operational 
equation. 

Of course, not all the operational equations so solvable have 
physical significance. Many however, do. Below is a list of 
infinite integrals with their known solutions, accompanied by 
the corresponding operational equation and its explicit solution. 
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All of these solutions are directly applicable to important tech- 
nical problems. It may be remarked in passing that the infinite 
integrals have for the most part been evaluated by advanced 
mathematical methods which need not concern us here. 


Table of Infinite Integrals, the Corresponding Operational Equa- 
tions, and Their Explicit Solutions 


(a) fre oS dba ty! 
—p-—nX 
ack ae 
(b) [ota 
= =i hn) 
(c) ye i eas 
: 0 Vat Vp 
f= A/D =) at. 
"a (20 dt 
@ J © 185: @n—1) Vat vrs 
tees (2t) 1 


p” °1.3.5---(n—1)/ qt 


2 00 fn 1 
=p aoe jp 
(e) {| ee dt erence 
Set eee Ue ey 
Ce ee 


“ past - 
[re Pot ens 
0 T ty/ t 


oO eke ag 2p 


a ota) ends 
e) 0 rt Vp 


es ON 
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Pp 
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: ae eae 
(i) || e-”* cos At dt Petr? 


h= aBe =cos At. 

(}) ee cos At dt See 
h= ant ane =e! cos MM. 

(ke) [omew sin Malm egy 
h= Canoe sg # aim XE, 


e-?*Io(dt)dt = = 
[em doQnd= 7 


Pp 
== =Jo(Mt). 
242 of ) 


Wi coe ea 

m cake V Pe 2 C= SS 

m) [endo oe 
h=—P__ ¢->Vv*F1=0 for t<r 


Vp+l oe 
=J)(~/#2—d?) for t>n. 


(n) || em Sq(ntydt=1(">P) rP=p+nr, 
0 ry 


Sais 
nah ; ) Jn(dt). 


(p) [ete o(ndat= 
0 
1 
h=—- 
Ay eee) /p 
In formulas (1), (m), (n), J,(v) denotes the Bessel function of 
order n and argument x. In formula (p), /o(~) denotes the Bessel 
function Jo(ix) where i=~+/ —1. 

This list might be greatly extended. As it is, we are in posses- 
sion of a set of solutions of operational equations which occur in 
important technical problems and which will be employed later. 

The foregoing emphasizes the practical and theoretical impor- 
tance of recognizing the equivalence of the integral and operational 
equations. With this equivalence in mind, the solution of an 


il 
V p?+2nrp 
=e T(t). 
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operational equation is often reduced to a mere reference toa table 
of infinite integrals. Heaviside did not recognize this equivalence. 
As a consequence many of his solutions of transmission line 
problems are extremely laborious and involved and in the end 
unsatisfactory because expressed in involved power series. 

Not all the infinite integrals corresponding to the operational 
equations of physical problems have been evaluated or can be 
recognized without transformation. This statement corresponds 
exactly with the fact that a table of integrals is not always suffi- 
cient but must be supplemented by general methods of integra- 
tion. We turn, therefore, to stating and discussing some general 
theorems applicable to the solution of operational equations. 

In the derivation of the operational theorems, which constitute 
the general rules of the operational calculus, the following proposi- 
tion, due to Borel and known as Borel’s theorem, will be fre- 
quently employed.! 

If the functions f(t), filt), and fo(t) are defined by the integral 
equations 


F(p) =| S(e->*at 
Fane {f fuer dt 
raw) = | fo(t)e”* dt 


and if the functions F, F, and F, satisfy the relation 
F(p) =F \(p).F2(p) 
then 


f(t) = {i (r)folt —1)dr 
=a ik (r)fi (t —r)dr. 
0 


The operational theorems will now be stated and briefly proved 
from the integral equation identity. 


1 For a proof of this important theorem the reader is referred to Borel, 
“Tecons sur les Séries Divergentes” (1901), p. 104; to Bromwich, “‘ Theory 
of Infinite Series,” pp. 280-281; or to Ford, “Studies on Divergent Series and 
Summability,” pp. 93-94 (being Vol. II of the Michigan University Science 
Series, published by Macmillan). The proof depends on Jacobi’s trans- 
formation of a double integral: see Edward’s “Integral Calculus,” 1922, 


Vol. II, pp. 14-15. 
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\ 


If in the operational equation 
h=1/H(p) 
the generalized impedance function H(p) can be expanded in a 
sum of terms, thus . 


1 1 1 if 
H(p) Hyp) Help) * Hap) 
and if the auxiliary operational equations 
1 

| Ne ges ecia 

Hp) 
h a 

* Hap) 


can be solved, then 
h=hythot-:-:>thny. 

This theorem is too obvious to require detailed proof: in fact 
it is self-evident. The power series and expansion theorem 
solutions are examples of its application. In general, however, 
the appropriate form of expansion of 1/H(p) will depend on the 
particular problem in hand. The theorem, as it stands is a 
formal statement of the fact that solutions can often be obtained 
by an appropriate expansion whereas the equation cannot be 
solved as it stands. 


THEoREM II 
If h=h(t) and g=g(t) are defined by the operational equations 
n= Li (op) 
g=1/pH(p) 


gt) = [rover 
0 


To prove this theorem we start with the integral equations 


1 co 
a hiije-?* dt, 
pil (p) J, 
il re co 
ee te-P* dt. 
pH (p) i) av 
The second of these is in form for an immediate application of 
Borel’s theorem since 


then 


1 C—¥ 1 . a <1) 
pip) ip pH (p) 
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The functions f; and f, of Borel’s theorem then satisfy the 


equations 
1 oe 
= te” di, 
Pp if flo 


sae) Om 


It follows at once that 
fi) =1 
fo(t) =h(t) 


whence by Borel’s theorem 


t 
o= | h(r)dr. 
0 
THEOREM III 


If h=h(t) and g=q(t) are defined by the operational equations 
h=1/H(p) 
g=p/H(p) 
then 
() =F) 
Uae gt 


provided h(0) =0. 
The integral equations of the problem are 


1 i>.) 
ee —pt d 
pH (p) J, ue te 


il Le) 
a = He? dt. 
H(p) J, — 
Integrating the first of these by parts we have, 


| 1 NY pee 
——~ =—h(0)+—]| A'(He?' dt 
pH (p) pr ) all @ 
where h’(t) =d/dt h(t). 

If h(0) =0, we have at once 


me (Ao 


Comparison with the integral equation for g(é) shows at once 
that g(t)=h'(t), since the integral equation determines the 


function uniquely. 
Theorems II and III establish the characteristic Heaviside 


t 
operations of replacing 1/p by | dt and p by d/dt. 
0 
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THrorem IV 
If in the operational equation 
h=1/H(p) 
the generalized impedance function can be factored in the form 
H (p) = Hi(p) - H2(p) 
and if the auxiliary operational equations 
hi=1/H,(p) 
h2=1/H2(p) 
define the auxiliary variables hy and he, then 


t 
h(t) Al hi(r)he(t—7)dr 


= a i T)hy(t—r)dr 
2 t U . 
al a : e 


This theorem is immediately deducible from Borel’s theorem and 
theorems II and III, as follows. 
The integral equations are 


ae eee “n(tenPt dt 
pH(p) pHi(p) pHa(p) pe 
1 “2 
— = | halt)e-P*dt 
pH (p) i UL 


1 92) 
Se | fee at, 
pH»(p) | ates 


Now define an auxiliary function g(t) by the operational 
equation 
it 


I= DH) 
1 l ie 
+ = t)e-P* dt 
pinto) pHa) J, 9% 
and by Borel’s theorem 


g(0= | na(e)natt—nar 


Then 


a 


t 
= | he(r)hi(t—r)dr. 
0 


From this equation it follows that g(0)=0, and hence comparing 


the operational equations for h and g, we have by aid of Theorem 
Il 


a) =" g(t 
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and hence 
t 
hee iI eed ede 
i), 


t 
=< , he(r)ha(t—r)dr. 

This theorem is extremely important, although not stated or 
employed by Heaviside himself. We shall make use of it in estab- 
lishing two important general theorems and shall have frequent 
occasion to employ it in specific problems occurring in connection 
with the subsequent discussion of transmission theory. 


THEOREM V 
If h=h(t) and g=g(t) are defined by the operational equations 


il 
h= 
H(p) 
ae 1 
9* H(p+n) 


where ts a positive real parameter, then 
t 
ga -tn i; dt)e-™ h(t). 
0 


To prove this theorem we start with the integral equations 


il Lo /2) 
i) =|] eee 


1 a 
pH(p+r)_ | BOE aets 
0 


In the first of these equations replace the symbol p by q+): 
we get 


1 1 e 
Se ee oe —tp—at 
th HGH J, Momma 
and then to preserve our original notation replace the symbol q 
by p, whence 


I = t et 
CERT i h(t)e™ e-P dt. (a) 


The integral equation in g(t) can be written as 


r 1 Coa Bs 
(+5 earner j, one a 
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Comparing equations (a) and (b) it follows at once from 
theorems I and II that 


g(t) = (1 +h | aU) 


From the foregoing, the following auxiliary theorem is immedi- 
ately deducible. 


THEOREM Va 
If h=h(t) and g=g(t) are defined by the operational equations 
jae 
H(p) 
g = eo 2 ns 
(p+r)H(p+n) 
gt) =h®e—™. 
The proof of this theorem will be left as an exercise to the reader. 


then 


THEeoREM VI 


If h=h(t) and g=g(t) are defined by the operational equations 


k= A@®) 
g=1/H (Xp) 
where \ 1s a positive real parameter, then 


g(t) =h(t/A). 
We start with the pai equations 


nl =[- h(t)e-P* dt 


Tan =(. g(t)e-P! dt 


and in the first of these equations we replace p by \q and t by 
7/X, whence it becomes 


1 oT 


Now replacing the symbols q and + by p and ¢ respectively, we 


have 
1 Go 
——— ee NX p—pt 
pH(p) i! h(t/r)e-”! dt 


whence by comparison with the integral equation in g(t) it 
follows at once that 


g(t) =h(t/r). 
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This theorem is often useful in making a convenient change in 
the time scale and eliminating superfluous constants. 
TurorEM VII 


If h=h(t) and g=g(t) are defined by the operational equations 


eel 
h=T(p) 
Brees 
I~ H(p) 


where d is a positive real quantity, then 
Ae) = Om tor ban 
=h(t—r) for t>n. 

This is a very important theorem in connection with trans- 
mission line problems where retardation, due to finite velocity 
of propagation, occurs. Its proof proceeds as follows: 

If the auxiliary function k=k(t) is defined by the operational 
equation 

k=e—” 


then by Theorem IV, 
d t 
g(t) = all k(r)h(t—7)dr. (a) 


Now, corresponding to the operational equation k=e~*? we 
have the integral equation 


ePe . 
—= {| k(t)e-”* dt. 
jy 0 
The solution of this integral equation, which is easily verified 
by direct substitution in the infinite integral, is 

k@)=0 for#<)d 

=1 fort>d. 
Hence equation (a) becomes 
g(t)=0 fort<r 


d t 
=h(t—d) for t>. 


Theorem IV, employed in the preceding proof, as stated above, 
is extremely important and we shall have frequent occasion to 
employ it in specific problems. We shall now apply it to deduce 
an important theorem which extends the operational calculus 
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to arbitrary impressed forces, whereas heretofore the operational 
equation h=1/H(p) applied only to the case of a ‘unit e.m.f.” 
impressed on the system. 

It will be recalled from a previous chapter that if x(¢) denotes 
the response of a network to an arbitrary force f(¢), impressed 
at time t=0, and if h(t) denotes the corresponding response to a 
“unit e.m.f.,” then 


art 
“=a | h(r)f(t—r)dr (31) 


and 


aie = | h(t)e-?* dt. (30) 


Now f(t) may be of such form that the infinite integral 


[so dt 
0 


can be evaluated and has the value F(p)/p: thus 
is 1 
{| f(He?' dt =— F(p). (55) 
0 P 


This is possible, of course, for many important types of applied 
forces, including the sinusoidal. 

It follows at once from Theorem IV that x(t) satisfies and is 
determined by the integral equation 


om | x(t)e* dt. (56) 
We have thus succeeded, by virtue of Theorem IV in expressing 
the response of a network to an arbitrary e.m.f. impressed at 
time ¢=0, by an integral equation of the same form as that 
expressing the response to a “unit emf.” That is to say 
we have, at least formally, extended the operational calculus 
explicitly to the case of arbitrary impressed forces. 

We now translate the foregoing into the corresponding 
operational theorem. 


Turorem VIII 
If the operational equation 


h=1/H(p) 
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expresses the response of a network to a “unit emf.” and if an 


arbitrary e.m.f. E impressed at time t=0, is expressible by the 
operational equation 


E=V(p) 
or the infinite integral 
i E(te-» at = 
0 P 


then the response x of the network to the arbitrary force ts given by 
he operational equation 
AVIS 
H(p) 
and x(t) is determined by the integral equation 


TeV 
D HeAe i ae(te™ dt. 


x 


THEeoremM 1X 
If the operational equation 
h=1/H(p) 
is reducible to the form 
Beas 
1+)K(p) 


where is a real parameter, and if the auxiliary functions f=f(t) 
and k=k(t) are defined by the auxiliary operational equations 


f=F(p) 
k=K(p) 


then h(t) ts determined by the Poisson integral equation 


d (*t 
h(t) =f(t) Nil h(r)k(t—r)dr. 


This theorem is of considerable practical importance in 
connection with the approximate and numerical solutions of 
operational equations when the operational equation and the 
equivalent Laplace integral equation prove refractory. In such 
cases, as will be shown later, the numerical solution of the Poisson 
integral equations can often be rapidly and accurately effected, 
and in many cases the qualitative properties of h(¢) can be 
deduced from it without detailed numerical solution. 
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The proof of this theorem proceeds as follows: 
By virtue of the relation h=1/H(p) the operational equation 


haa Oe 
1+)K(p) 
can be written as 


K(p) _ 
Ue ae =F (p) 


K(p) 
h=F(p)-—\=*: 
(p) H(p) 
A direct application of Borel’s theorem or Theorem IV gives at 
once the explicit equivalent 


a(t) =f) 15, i REG Ree 
J0 ‘ 


The preceding theorems, together with the power series and 
expansion theorem solutions formulate the most important rules 
of the operational calculus, and are constantly employed in 
the solution of electrotechnical problems. On the other 
hand, the table of infinite integrals furnishes the solution of a set 
of operational equations, which are of the greatest usefulness in 
the systematic study of propagation phenomena in transmission 
systems which will engage our attention. Before taking up 
this study, however, we shall first solve a few specific prob- 
lems which will serve as an introduction to asymptotic and 
divergent solutions involving Heaviside’s so-called “fractional 
differentiation.” 

Problem A: Current Entering the Non-Inductive Cable.— 
The non-inductive cable is a smooth line with distributed resist- 
ance FR and capacity C per unit length; for the present we neglect 
inductance and leakage. A consideration of cable problems 
leads to some of the most interesting questions relating to 
operational methods, particularly to questions regarding diver- 
gent expansions. It would seem best to allow specific problems 
to serve as an introduction to these general questions, 

The differential equations of the cable are 


ss 
igl= cre 
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where x is the distance, measured along the cable from any fixed 
point, J is the current at point x, and V the corresponding 
potential. 


Replacing d/dt by the operator p, we have 


0 
RI= Seria 


(58) 
a 
pCV = —— I. 


Eliminating, successively, V and J from these equations, we 
get 


Rey 
lea 
and 
a2 


These equations have the general solution 


V =Vye-%7+ Ve (59) 

r= [PC rie Vae (60) 
where hase 

y=V RC. (61) 


The term in e€”™” represents the direct wave and the term in 


e” the reflected wave. V, and V2 are constants which must be 
so chosen as to satisfy the imposed boundary conditions at the 
terminals of the cable. 

For the present we shall assume that the line is infinitely long 
so that the reflected wave is absent. We shall also assume that 
a voltage H is impressed directly on the cable at x=0: we have 
then, 

V =Ee-tV0CR = Re-V ew (62) 


1a PC Bev 008 = PC pv (63) 


where a denotes «?RC. 

To convert these to operational equations let us suppose that 
E is a “‘unit e.m.f.” (zero before, unity after time t=0). We 
have then, in operational notation 


V =e Van (64) 
[= 4 eoe-vae (65) 
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Now suppose that +=0 so that a=0, in other words consider 
a point at the cable terminals. Then 


I = (26. (66) 


The first of these equations means that V is simply the impressed 
voltage, zero before, unity after time t=0, as of course, it should 
be from physical considerations. 

Corresponding to the operational equation 


r=). (66) 


we have the integral equation 


VEve if “T(de-*dt. (67) 


The solution of this is known (see formula (c) of the preceding 
table of integrals): it is 


i: Vacs Ay hee a (68) 


Heaviside arrived at this solution from considering the known 
solution of the same problem in the theory of heat flow. He 
therefore inferred that the operational equation 


I=V>p 
T=1/Vri. 


This is correct; we, however, have derived it directly from the 
integral equation of the problem and the known integral 


Taal le: (69) 
Vp 0 Vat 

We thus see from the foregoing that, if a ‘unit e.m.f.’’ is 
impressed on the cable terminals, the current entering the cable 
is initially infinite and dies away in accordance with the formula 
VC/Rt. The case is, of course, idealized and the infinite 
initial value of the current results from our ignoring the dis- 
tributed inductance of the cable, which, no matter how small, 
keeps the initial current finite, as we shall see later. 

Now let us go a step farther; suppose that in addition to 
distributed resistance R and capacity C, the cable also has 


has the explicit solution 
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distributed leakage G per unit length. The differential equations 
are now 


(70) 
0 


Consequently it follows that in the operational equation for the 
current entering the cable we need only replace Cp by Cp+G. 
Therefore, when leakage is included, equation (66) is to be 
replaced by 


Se ECAG © ve a 
T=4| = = RV Pts 71) 
where \=G/C. 
The corresponding integral equation is, of course, 
CV p+n i “ 
= I(tye-?' dt. 72 
ara Re (72) 


We shall give two solutions of this problem; first the solution of 
the integral equation, and second the typical Heaviside solution 
directly from the operational equation. 

Equation (72) may be written as 


CU) a ae 
I(t)e-” dt. (73) 
VE Vj Vp+r 0 
Now suppose that J(¢) is the solution of the equation 
i 
J (the! dt 74 
V ptr “Il oe, 


it follows at once from Theorems (I) and (II) of the preceding 


chapter that 7 
I(t) -/S( 354 j au). (75) 


Also from formula (c) of the table of integrals and Theorem 
(Va) the solution of (74) is 


I= - (76) 


10-51, | NS, wae (77) 


whence 
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The integral appearing in (77) can not be evaluated in finite 
terms; it is easily expressible as a ‘series, however, by repeated 


integration by parts. Thus 
—rt t ~ fe t = 
reo 2 [lonavin2vi ewan | em v/t dt. 
oVi 0 0 
Proceeding in this way by repeated partial integration we get 
for the integral term of (77) 


= 2 
avior (148 ONT... (78) 


The straightforward Heaviside solution is obtained by expand- 
ing the operational equation as follows: 


- (Sa) v 
= pl1t5 3 mG y tase) ve 


Identifying ~/p with 1/+/zt (from known solutions of allied 
problems) and substituting for 1/p" multiple integrations of the 
nth order we get 


pal OF (2M) _(2M)*, 1.3(2M)*_ i 
‘ Va { +2 934934567 | Ce 


It can be verified that this solution is convergent and equivalent 
Ose doe 

This problem, while simple and of minor technical interest, 
will serve to introduce us to the very important and interesting 
question of asymptotic series solutions. 

An asymptotic series, for our purposes, may be defined as a 
series expansion of a function, which, while divergent, may be 
used for numerical computation, and which exhibits the behavior 
of the function for sufficiently large values of the argument. 

Let us return to equation (77). We observe that the series 
solution (78) of the definite integral becomes increasingly labor- 
ious to compute as the value of ¢ increases. This remark applies 
with even greater force to the Heaviside solution (79) on account 
of the alternating character of the series. Right here we have 
an excellent example of what I regard as Heaviside’s exaggerated 
sense of the importance of series solutions as compared with 
definite integrals. Consider the solution in the form of (77) as 
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compared with Heaviside’s series solution (79). The former is 
incomparably easier to interpret and to compute, either by 
numerical integration or by means of an integraph or planimeter. 
In fact the series (79) is practically unmanageable except for 
small values of f. 


Returning to the question of an asymptotic expansion of the 
solution (77), we observe that the definite integral appearing in 
that equation can be written as, 


oon jaa { ois ie. dt (80) 


provided is eet as it is in this case. Now the value of the 
infinite integral is known; it is ~/z/\. Consequently 


ee ata alee ame ee dt: (81) 


e™ iG feet bee pan? ee OY Bete ae 
et = | 46 = =Oh 
svt MP Oe vt all t/t 
Integrating again by parts we get 
ice Lee ola Gee 
oo dt 
AV/t 2QW2tr/t a erele ; ae 
Continuing this process, we get 
fee, Ga E i 3 135 


furthermore, 


Leia Bees: 


2rt (2nt)? (2nr1)8 
1.3.5. ..(2n—1) 
sae (— L) “(2nt)" 


el reper etl) | ee lag 
d “22X)" , Pte 
Now this series is divergent, that is, if we continue out far 
enough in the series the terms begin to increase in value without 
limit. On the other hand, if we stop with the nth term the error 
is represented by the integral term in (82) and this is less than 


(a) 1.3.5 (Qn), =, 
Wt Quy se, 


Consequently the error committed in stopping with any term in the 
series is less than the value of that term. Therefore if we stop with 
the smallest term in the series, the error is less than the smallest 
term and decreases with increasing values of t¢. 


(82) 
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We can therefore write the ra (77) as 


A Ke é Ape 3.5 

\'s +e fe ~ ant)? * ney? CO 
The first term, since \=G/C, is simply »/G/R, the d.c. admit- 
tance of the leaky cable. The divergent series shows how the 
current approaches this final steady value. 

In this particular problem no asymptotic solution is derivable 
directly from the operational equation, at least by the straight- 
forward Heaviside processes. Asymptotic solutions, however, 
constitute a large and important part of Heaviside’s transmission 
line solutions. We shall therefore discuss next a problem for 
which Heaviside obtained both convergent and divergent series 
expansions. 

Problem B: Terminal Voltage on Cable with “Unit E.M.F.”’ 
Impressed on Cable through Condenser.—We now take up a 
problem for which Heaviside obtained a divergent solution, and 
which will introduce us to the theory of his divergent solutions 
and so-called ‘‘fractional differentiation.”” We suppose a ‘‘unit 
e.m.f.’’ impressed on an infinitely long cable of distributed resist- 
ance F& and capacity C per unit length through a condenser of 
capacity Cy: required the voltage V at the cable terminals. 
The operational equation of the problem is derived as follows:— 

We know from the problem just discussed that the current 
entering the cable whose terminal voltage is V, is, in operational 


notation 
\ Gas 
But the current flowing into the condenser is 
Cop(l—V) 
since the voltage across the condenser is 1—V. Equating these 


two expressions we get 


= peo 
pias eae —— 85) 
796; ite /pC/R \ 
which is the operational equation of the problem. 
This may be written as 
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where 
Ps lets = 
Va=o VC/R. 
Now expanding this by the binomial theorem 
2 
Ve=l— ee a eae +(2) i, 
P Pp p P 
=149+(2) +. 
—a4e+(2) +. wl 


Sinope eae (86) 
7 oat (at) (2at)* = 
( in nO) 


by the usual Heaviside rules of ‘‘algebrizing.”’ 
It is worth while verifying this from the integral equation of the 


problem. We have 
iL 1 a 
= = V (t)e?* dt. 87 
P1+Va/p \| ee 


The left hand side can be written as 


det ae 
De es ND) 


and by the formulas and theorems given in a preceding section the 
solution can be recognized at once as:— 


V(t) =e" — (88) 


a eut 
Tv 


This can also be written as 


vinmaftet [rar (89) 


If the definite integral of (88) is evaluated by successive partial 
integrations it will be found in agreement with the Heaviside 
solution (86). 

Now the solution (86) is in powers of ¢ and while absolutely 
convergent becomes progressively more difficult to interpret 
and compute as the value of é increases. From (89), however, 
we can derive a divergent or asymptotic solution applicable both 
for interpretation and computation, when the value of ¢ is suffi- 
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ciently large. Asin the example discussed before, the asymptotic 
expansion results from repeated partial integrations; thus 


ie ie eal ane! 
wee 4 t V7 
nienens vl ce 


pelt = see 


OG 


Ton —=de-* 
"avi 1 TV T 
Cam at e —at il. 3 eau 
eel eae <0 
avV/t 2a2tv/t 2a” al Ta/7 
and finally 
Ca i 1G oe es) ] 
pate Sy ee a ae Oe ee ee 0 
av/t | Te oat ane (Qat)3 | C0) 


The series (90) is divergent just as is (82) of a preceding problem 
and the error committed by stopping with the smallest term, is of 
the same character and subject to the same discussion. With this 
understanding we write the solution (89) as 

1 1 | es ae Be es | 
Va 1 att (Qat)? (at)? t [ ee 
For large values of ¢ (a#>5) this series is accurately and rapidly 
computable. Furthermore it shows by mere inspection the 
behavior of V(t) for large values of ¢t, and that it ultimately 
approaches zero as 1/+/rat. 

Let us now see how Heaviside attacked this problem and how 
he arrived at a divergent solution from the operational formula. 
Returning to the operational equation (85), it can be written as 


ya VV p/o 
1+V/p/ /a 


Now expand the denominator by the binomial theorem: we get 


formally 
: ; 2 ) 
oe ne D4 B_P [P(2) i te 
a a a a a a 
eae ata ee 
(1+2+(?) +... )a/? (3 


VO 


(92) 
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Heaviside’s procedure at this point was as remarkable as it was 
successful. He first discarded the second series in integral powers 
of p as meaningless. He then indentified ~/ p with 1/ / at 
and replaced p” by d”/dé” in the first series, getting 


il@ . il ¢ il 
V=({1+-— Sal Se 9 eee 4 
( ue di’ a dt? cc Nae ae 


or, carrying out the indicated differentiation, 


ye (-s tana) 
/ rat 20» (Zan) (2a). 
which agrees with (91). 

This is a typical example of a Heaviside divergent solution for 
which he offered no explanation and no proof other than its 
practical success. His procedure in this respect is quite unsatis- 
factory and in particular his discarding an entire series without 
explanation is intellectually repugnant. We shall leave these 
questions for the present, however; later we shall make a syste- 
matic study of his divergent solutions and rationalize them in a 
satisfactory manner. First, however, we shall take up a specific 
problem for which Heaviside obtains a divergent solution without 
discarding any terms. 

Problem C: Current Entering a Line of Distributed L, R and C. 
Consider a transmission line of distributed inductance L, 
resistance FR, and capacity C per unit length. The differential 
equations of current and voltage are 


d Ouee 
Erp ees Fa 
(95) 
a eee 
dt One 
Replacing d/dt by p, we get 
(pL+R)I=—7-V 
. (96) 
Cp-V=—-<1. 


Equations (96) correspond exactly with (58) for the non- 
inductive cable, except that we must replace R by pL+R. For 
the infinitely long line, therefore, the operational formula for 
the current entering the line is 


Tae eV, (97) 
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where V» is the voltage at the linesterminals. If this is a “unit 
e.m.f.”’ we have, as our operational equation, 


jeg ee (98) 


which can be written as 


C 1 
L/1+2n/p 
where \=F/2L. 
The corresponding integral equation is 
Or i. 
eee e-Pt I t dt. (100) 
L/p+2rp Jo i 


From either equation (99) or (100) and formula (p) of the table of 
integrals, we see at once that the solution is 


T=a[Se-™ 100) (101) 


where Jo(Aé) is the Bessel function Jo(iAt), where i=+~/—1. 
(The function is, however, a pure real.) 

Heaviside’s procedure, in the absence of any correlation 
between the operational equation and the infinite integral, was 
quite different. Remarking, with reference to equation (99), 
that “the suggestion to employ the binomial theorem is obvious,”’ 
he expands it in the form 


oe agi bey Oo ge beyay eae 
fee ji ) ~135(%) ie (102) 


and replaces 1/p" by é"/n! in accordance with the rule discussed 
in preceding sections. The explicit solution is then 


1.3 1.3.5 

r=/6 L=N+ (Qa)? (ga Qe)? 3 (103) 
a convergent solution in rising powers of ¢. As yet, however, he 
does not recognize this series as the power series expansion of 
(101), which it is. He does, however, recognize the practical 
impossibility of using it for computing for large values of ¢, 
and remarks “ But the binomial theorem furnishes another way 
of expanding the operator (operational equation), viz. in rising 
powers of p.”” Thus, returning to (99), it can be written as, 


p=, |C VP/2r_ 


L/itp/2x 0s) 
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Now expand the denominator by the binomial theorem: we get 


C | 13 (pt) ae ls 5/0 || 
rayfe| mat: +2) - Bl NaN pe Peete) 
He now eae V/p/2 with 1/+/2rM and replaces p” in the 
series by es dt", thus getting finally 


1 1 (1.3)? , (1.8.5)2 
r / Int | IT gt aicgnty?? 3(SM)* oe uy, 
This series solution is divergent: Heaviside recognizes it, how- 


ever, as the asymptotic expansion of the function e~™ J)(At), and 
thus arrives at the solution 


[= af Se Tolas) (101) 


which we have obtained from our table of integrals. 

Now the divergent expansion (106) is the well known asymp- 
totic expansion of the function e~ “ J)(t), which is usually derived 
by difficult and intricate processes. The directness and simplic- 
ity with which Heaviside derives it is extraordinary. 

We note in this example that no integral powers of p appear 
in the divergent expansion: consequently no terms are discarded. 
Otherwise Heaviside’s process is as startling and remarkable as 
in the example discussed in the preceding section. 

We shall later encounter many problems in which asymptotic 
solutions are derivable as in the preceding example. We have 
sufficient data, however, in these two typical examples to take 
up a systematic discussion of the theory of Heaviside’s divergent 
solution of the operational equation.+ 


1 Quite recently an entirely different treatment of this whole problem, 
based on the Fourier Integral identity, has been developed by Dr. Norbert 
Wiener. A brief report on this work was made before the American Mathe- 
matical Society in May, 1925; it has since been published in Math. Annalen, 
Bd. 95, Heft 4, 1926. 


CHAPTER V 


THE THEORY OF THE ASYMPTOTIC SOLUTION OF 
OPERATIONAL EQUATIONS 


A study of Heaviside’s methods, as exemplified in the preceding 
examples and in many problems dealt with in his Electromagnetic 
Theory, Vol II, shows that they may be divided into two classes: 
(I) those of which the operational equation is of the form 


h=F(p)V/p (I) 
and (II) those of which the operational equation is of the form 
h= (p'v/p) (II) 


where & is an integer. 

Heaviside himself does not distinguish between the two classes, 
but employs the following rule! for obtaining asymptotic expan- 
sion solutions: 

If the operational equation 


h=1/H@) 
can be expanded in the form 


R= ao Gy p02 p?--* "Ga po = * 
(bo+bip+bap?+:+bap"+---)Vp, (107) 


a solution, usually divergent, is obtained by discarding the first 
expansion entirely, except for the leading constant terms ao, replacing 
Vp by 1/+/xt and p™ by dv/dt" in the second expansion, whence 
an explicit series solution peg 


1 
pee sia i dt? een . Me (108) 
=ay+— 


1.3.5 
Salo athe g a Os one Tt ) oe 


The present chapter is devoted to the mathematical theory of 
this type of expansion solution of the operational equation, and 
to an investigation of its practical and theoretical significance. 

1Tt does not appear that Heaviside makes a formal statement of this rule. 
It is a generalization, however, of his procedure in a large number of specific 
problems. 

62 
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Heaviside himself gives no information which would serve us as 
a guide in informing us when the rule is applicable and when it is 
not. Consequently it becomes a matter of practical importance, 
not only to investigate the underlying mathematical philosophy 
of the rule and to establish it on the basis of orthodox math- 
ematics, but also to develop if possible a criterion of its appli- 
cability. In this investigation we shall have recourse to the 
integral equation of the problem. 

We shall take up first the type of problem (Class I) in which the 
operational equation is 


h= ie ) =F (p)Vp (110) 
and assume that F'(p) admits of the formal power series expansion 
F(p)=bo+bip+bep?+b3p3+ - -- (111) 

The corresponding integral equation is 
|e erat (112) 


Vp 0 


We now assume the existence of an auxiliary function k(é), 
defined and determined by the auxiliary integral equation 


F(p)= il Ope (113) 


Now since 
1 = dt 
asl ——— (114) 
Vp 0 rt 


it follows from (112), (113), and (114) and Borel’s Theorem, or 


Theorem IV, that 

ey eG) 

h(t) = —— dr. ~ (ks 
Oa | ae (115) 
Now if we differentiate (113) repeatedly with respect to p and put 
p=0, it follows from the scaraenwen that 


b,=(—1)" Hb) (116) 


This equation presupposes, it Sat be noted, the convergence of 
the infinite integrals for all values of n, and therefore imposes 
severe restrictions on k(¢) and hence on F(p). We shall suppose 
that these restrictions are satisfied and discuss them later. 

Now (115) can be written as:— 


es a, 7)(1—7/t)-”2 
n= | der@)C [ty (117) 
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It can be shown that, if /(¢) satisfies the restrictions underlying 
(116), the integral (117) has an asymptotic solution obtained 
as follows:—Expand the factor (1—7/t)’? by the binomial 
theorem, replace the upper limit of integration by # , andintegrate 
term by term: thus 


11 1?t 
nO~F | i katy, [ eCOat 


oe fee 
earl ER()dt-+- |. (118) 
Finally from (116) we get 
1 1 1.3 1.3.5 \ 
ever {bob toy tbe (22°? (apa or | (119) 


which agrees exactly with the Heaviside rule for this case. 

The foregoing says nothing regarding the asymptotic character 
of the solution. It is easy to see qualitatively, however, that 
(118) and therefore (119) does represent the behavior of the 
definite integral (117) for large values of ¢, provided k(é) con- 
verges with sufficient rapidity. 

The foregoing analysis may now be summarized in the following 
proposition; which will be called Rutz I. 

If the operational equation h=1/H(p) ts reducible to the form 


h=F(p)Vp 
and if F(p) admits of power series expansion in p: thus 
F(p) =botbi pbs p?+ Nee +bnp?+ oieale 
so that formally, 
h=(botbiptbsp2+ - + +bap"t+ - - )Wp 


an explicit series solution, usually asymptotic, is obtained by replac- 
ing Vp by 1/+/xt and p(n Sa by d"/dt", whence 


il 
h(t)™{ bo+b 22 \——— 
() — +b ot ) a 
1 13.0 
~ amet ag tbe apes aye > 


provided the function k=k(t), defined by the operational equation 
k=F(p), and the infinite integrals 


[exoat (n=1,2,....) 
0 


exist. 
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We shall now apply the foregoing theory to a physical problem 
discussed in the last section: namely, the current entering an 
infinitely long line of inductance L, resistance R and capacity C 
per unit length. It will be recalled (see equation (100)) that the 
integral equation of this problem is 


@ 1 i . 
eee rei(Out 
TAs te oe 


where \=#/2UL, and that the solution is 
C 
= per Too), 


We can derive the solution in another form appropriate for our 
purposes by writing 


Oa aay {| eI (A) dt 
0 


L~/pVp+2r 
Now since 
1 a . dt 
/p 0 ! Vat 
and 
1 oo ena 


sd i 
Sree Nhe ee 


it follows from Borel’s er that 


e72Ar 
————_ dr. 
Voz, Viv t—7 
Now subject this definite integral (omitting the factor »/C/L) 
to the same ae joer to (117): we get 


L gage 4 2nt 
ee e7 2M di 
Ay, /t atts ie 
Te iia/e 
pan Se Ay 


The infinite integrals are known and have been evaluated. Sub- 
stituting their values this series becomes :— 
ie Hl 17557 AN ae 
Tea Tet aie? sens 

which is in fact the well known asymptotic expansion of the 
function e~™ I(nt). 

A second example may be worth while. Consider the case 
of an e.m.f. e~- impressed at time t=0 on a cable of distributed 
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resistance R and capacity C: required the current entering the 
cable. The required formula is! ~~ 
Crd eo 2 
Hea ina hy ey as 


= oe ee af, att] om 


by obvious transformations. 
Asymptotic expansion of the definite integral as in the preced- 
ing example gives the asymptotic formula 
eC. | 1 ee ee 
Be ee Greene Tonpet | 
The operational formula of the problem is 


a “e caer 


_ jC pyr 
= 6 Powe 


EO) =)'--|ve 


Applying the Heaviside Rule, we get the asymptotic expansion 
[ca Es ee ee | 

rRt|2nt' (@nt)?* nat | 

which agrees with the preceding formula, derived from the 
definite integral. 

We shall now discuss a specific problem in which the Heaviside 
Rule breaks down. For example let us take the preceding prob- 
lem, and replace the applied e.m.f. e~” by sin wt. The formula 
corresponding to (120) is now 


ee ee a2) 


If we now attempt to expand the definite integral of (121) in the 

same way as that of (120), we find that the process breaks down 

because each component of the infinite integral is now itself 

infinite. In fact no asymptotic solution of this problem exists. 
Let us, however, start with the operational formula: since 


“s A) 
e-”' sin wtdt =——— 
0 p?+w? 


1 The derivation of the formulas in this problem is left as an exercise for the 
reader. 


l= 


ja 
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Fe eo aV D- 


Now expand this in accordance with the Heaviside Rule: we get, 
operationally, 


se ORT Eo Vein @ 
rGl()-€)+0)- Jus 
and explicitly 
- Cot te 135 | 
us: Nee ret 
which is quite incorrect.! The incorrectness of the result will 
be evident when we remember that the final value of the current 
is the steady state current in response to sin wt, or 
mz 
“) oF (cos wt-+sin wt). (122) 


This result can be derived directly from (121) by writing it as 
Tas “eos ot {att sin Af ental 


If the time is made indefinitely great the upper limits of the 
integrals may be replaced by infinity. The infinite integrals 
are known: substitution of their known values gives (122). 

This example illustrates the care which must be used in apply- 
ing Heaviside’s rules for obtaining divergent solutions and the 
importance of having a method of checking the correctness of 
his processes and results. 

We now take up the discussion of the asymptotic expansion 
solution of operational equations of the type 


h= 6(p*v/p) (k integral). (123) 
The Heaviside Rule for this type of operational equation may 
be formulated as follows (Ruuz II): 
If the operational equation h=1/H(p) ts reducible to the form 
h= ¢(p'v/p) 
and if ¢ admits of power series expansion in the argument, thus 
h = dota, p* VJ ptasp**1+a3 ptt VJ p+ ce 


1 While this series is incorrect as an asymptotic expansion of the current it 
has important significance, as we shall see, in connection with the building up 
of alternating currents. 
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a series solution, usually divergent and asymptotic, 1s obtained by dis- 
carding integral powers of p, and writing 

h=Go4- (1p +039" "4-0, po Vp. 
The explicit series solution then resulis from replacing Vp by 
1/+/rt, and p by d"/di, whence 


q3k+1 (ok+2 ii 

ke Botan aaa qpeee tS VA 

ee (—1)* iO ean(ai— 1) 1.3...(6k+1) 
dot al (a® (2t)* a3 (2t)3e+1 ae ) 


The theory of this series solution will be based on the follow- 

© .—pt ag 

ing proposition, deducible from the identity i Sit /V p. 
0 Vat 


If the function F(p) of the integral equation 
Ft)= [serra 
0 


approaches 1/+/p as p approaches zero, then f(t) ultimately behaves 
as 1/+/zt: that is, if F(p)1/+/p as p— 0, then f(t)~1/+/zt as 
t— oo, provided that f(t) converges to zero, and contains no term or 
factor which is ultimately oscillatory. 
To illustrate what this condition means suppose that 
fine a ia COS wt 


/ xt See at 


then 
i) f(Hetdt >a/+/p as p— 0, 
0 


and the oscillatory term in f(t) converges to a higher order. The 
presence of such oscillatory terms vitiate, therefore, the Heaviside 
Rule: in the following discussion we shall assume that they are 
absent. 

We are now prepared to discuss the operational equation 

h=$(p'v/p) 
and for convenience shall assume that k =0 so that the operational 
equation becomes 
= $(v/p) 


of which the corresponding or equivalent integral equation is 


nov p= { h(t)e-”* dt. (123b) 


THE THEORY OF THE ASYMPTOTIC SOLUTION 69 


We assume that $(+/ p) admits of formal power series expansion 
in the argument: thus 
o(Vp) =do+a1 V/ptarptas pV p+as | Ot ak 
without, however, implying anything regarding the convergence 
of this expansion. 
We now introduce the series of auxiliary functions, g,91,g2, 


gs, . . . defined by the following scheme 
g(t) =h(t) —ao 
ay 
No S= 
gi(t) =g(2) va 
Gs 1 as 
galt) =t.ilt) +5 = 
= LEMS, 
a(t) =t90 ~ ge (123¢) 
1.355) G7 


ga(t) =t.gs(t) + 


2° A/ at 


Successive substitutions in the integral equation (123b) and 
repeated differentiations with respect to p, lead to the set of 
formulas, 


a) ay 
g(t)e-?' dt~—= as p— 0 
J, Vp 


oo a3 
t.gi(t)e"' dé~— ——~ as p> 0 (123d) 
J, 2V/p 
os Hee as 
t.go(t)e—”' dt~—, —j=as p—0 
J 2° /p 


rd 1.3.5 az 
t.ga(t)e7?! di — ———— —Las p90 
J DENT 


Sel els: te) os, w) so) lo Sa) 0, @) ele) @) el mle! (a) po: 8) 6) 6: 1m (6) 2) ope Ae 


Now assuming that h(t) satisfies the restrictions stated in the 
preceding proposition, it follows from that proposition, that 


g(t)~a1/W rt as t—> 0% 


ROE 7a 88 t—>"00 (123e) 
Tv 
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From the set of equations (123d) and (123e) it follows by 
successive aay that ‘ 


tee ale 1.3.5 
MO) ~a0t (a 3 5,145 ene" pat’ : ) 


which agrees with the series gotten by applying the Heaviside 
Rule. 

The defect of this derivation, which, however, appears to be 
inherent, is that it requires us to know or assume at the outset 
that h(t) satisfies the required restrictions. Consequently an 
automatic application of the Heaviside Rule may or may not 
give correct results. On the other hand if we know that an 
expansion solution in inverse fractional powers of ¢ exists, the 
Heaviside Rule gives the series with extraordinary directness 
and simplicity. | 

The type of expansion solution just discussed will now be 
illustrated by some specific problems. The first problem is that 
of the propagated voltage in the non-inductive cable in response 
toa “unit e.m.f.” It will be recalled that in a preceding chapter 
we derived the operational formula 

V=e-Ver (124) 
where a=2?RC, for the voltage at distance x from the terminal 
of a non-inductive cable of distributed resistance R and capacity 
C, in response to a “unit e.m.f.’’ impressed at pointw=0. Heavi- 
side’s solution of this operational equation proceeds as follows: 

Expansion of the exponential function in the usual power 


series gives 
V ap ap _oPV ap (ap)? _ 
1! 2| 3! A! 
which may be rearranged as 


y= 1—(14 8p oP Van 
aps (ap) 7 (ep): ee 
+(P4 Ale ee 
Heaviside now discards the series in integral powers of p entirely, 
replaces +/p by 1/+/zt and p" by d"/dé in the first series, and 
thus gets 


a? q? a 
a (14 epee ) rt 


oy eee a fe LS (a \e Seopa | 
ss 20 He) Hats) a @ eta (126) 


V=i—- 
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Forme i etl a\ Gel ae \ommlayo\* 
re val ke) rela) +-): (127) 


This solution is correct, as will be shown subsequently. 

A rather remarkable feature of this solution—a point on 
which Heaviside makes no comment—is that it is absolutely 
convergent. In other words, a process of expansion which in 
other problems leads to a divergent or asymptotic solution, here 
results in a convergent series expansion. 

To verify this solution we start with the corresponding integral 
equation of the problem 


1 Van = if V (t)e-?* dt. (128) 
Pp 0 
It follows from this formula and theorem (II) that 


vo | aoa 


where ¢(¢) is determined by the integral equation 
ven | o(te”* dt. 
0 


Now from formula (f) of the table of integrals 


or 


eval |e a ae ee 
2NiaJo t 
whence 7 
Sl Ke e7e/st 
$0 =r 7 
and finally 


V= ade, where t’ =4t/a. (129) 


Alec 


To convert this to the form of (127) we write 


de (ees Lf ter 
Sasa |) eee 130 
V@=| oe Fl, ae oe 


The value of the infinite integral is known to be ~/z so that 
1 


a ae 
V=1- isi 
alle ee ey 


Now in the integral term of (131) expand e~" in the usual 


exponential power series and then integrate term by term; the 
series solution (127) results. This series, while absolutely con- 
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vergent, is difficult to compute for small values of t; an asymptotic 
expansion, which can be employed for computation for small 
values of ¢t is gotten as follows: 

Write (129) as 


t! ; 1 eee 
=~ /—e-1/? — — | ——dr. 
Ne 2V/rJo0 Wr 


Repeated partial integrations of this type lead to the series 


Vaa)bewe | 1-(E)41a(6)'—--}. 32) 


It is interesting to note, in passing, that an asymptotic solution 
of this type does not appear to be directly deducible from the 
operational equation. We observe also that, in this problem, 
the series in inverse powers -of ¢ is convergent while the series 
in ascending powers of / is divergent: the converse is the case in 
the problems discussed previously. 

A second specific problem may be stated as follows: 

Let a ‘unit e.m.f.’’ be impressed on an infinitely long non- 
inductive cable of distributed resistance R and capacity C per 
unit length through a terminal resistance Ro: required the voltage 
V on the cable terminals. The formulation of the operational 
equation of this problem is very simple. It will be recalled 
that the operational formula for the current entering the cable 
with terminal voltage V is V\/Cp/R. But the current is clearly 
also equal to (1—V)/o: equating these expressions we get 


+o = VV pG/R 

e 0 

whence 
aoa! Re 
~Vp/ntl ids 


where 1/+/A=Ro\V/C/R. This is the required operational 
formula. 

To derive the Heaviside divergent expansion, expand (133) by 
the binomial theorem: thus 


VSL p/N(B/\) =D. N) ere 
=1—(1+p/d+(p/d)2+--)Vp/r (134) 
+ (p/ A+ (p/ rd)? + (p/ A) 


THE THEORY OF THE ASYMPTOTIC SOLUTION 73 


Discard the second series in integral powers of p; replace /p 
by 1/+/zt and p" by d"/dt" in the first series, thus getting 


x 1 
V=1- (hae Aa (135) 


7 ihe 
3 Talis nT ON ) 8e) 


which is the asymptotic solution of the problem. 
To verify this solution we shall consider the more general 
operational equation 


| 
A a 
a form of equation to which a number of fairly important 
problems is reducible. (The parameter of equation (133) can 
be eliminated from explicit consideration by means of theorem 
VI.) 
Multiplying numerator and denominator of equation (137) 
by p"+/p—1, it becomes 


(n integral) (137) 


ae A asd Bee fe 


One a ptt — 
p ife 


(138) 


pett—] 1 


and by direct partial fraction expansion, this is equivalent to 


= 27 
ie Vp ae 


an+1l< ae 75 ~ Qn+1 P—Dm ee) 
0 


where 
. 2mr 


Dm=e2tl (m=0,1,2...2n). 


Write, for convenience, 


and consider the operational equation 


1 Ds. = Pm ) 
hn= eA): 140 
sens” eu? P—Pm ca 
By the rules of the operational calculus, fully discussed in 


preceding chapters, the solution of this is 


au topm(t—r) 


Nps 
hal ari i We 


dr +1 ~er') (141) 
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We have now to distinguish two cases: (1) when the real part 
of Dm is positive, and (2) when the,real part is negative. 
Taking up case (1) first, the preceding can be written 


il ion t oPmr 
hm (t) =——~| 1+? m'\ — dr—1 142 
SO al Cen val aa 
n+1 OO. 

{p é Pmt \ 
Ltermts fm eireat 
Pen. : a Aer 0 Vr 2 J 
yi 2n + 1 Dee K * =P mT 


dr 
Jr t Wr 
=— ——™_ EP m = F 144 
rane ens 4 \/r i Une 


Repeated integration by parts of the definite integral leads to an 
asymptotic series, identical with that obtained by applying the 
Heaviside Rule to the operational equation (140). 

If, on the other hand, the real part of pm is negative, we write 


(141) as 
7 ern 
hm) =s- Dae if ePmr } (145) 
al ee : Ves 

The term e””’ ultimately dies away, and the definite integral 
can be expanded asymptotically in accordance with the theory 
discussed under Rule 1, again leading to an asymptotic series 
identical with that given by direct application of the Heaviside 
Rule to the operational equation. 

Consequently since the operational equation in h» can be 
asymptotically expanded by means of the Heaviside Rule, the 
operational equation in h= Din is similarly asymptotically 
expandible, and the Heaviside Rule is verified for equation (137). 

We have now covered, more or less completely, the theoretical 
rules and principles of the operational calculus in so far as they 
can be formulated in general terms. We shall now apply these 
principles and rules to the solution of important technical 
problems relating to the propagation of current and voltage 
along lines. In doing so, while we shall take advantage of our 
table of integrals with the corresponding solutions of the oper- 
ational equation, we shall also sketch Heaviside’s own methods 
of solution. 

We shall close this discussion of divergent and asymptotic 
expansions with a general expansion solution of considerable 


(143) 


t 
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theoretical and practical importance in the problem of the 
building-up of alternating currents. It will be recalled from 
Theorem VIII that the response of a network of generalized 
operational impedance H(p) to an e.m.f. H(t) impressed at time 
é=0 is given by the operational formula 


a) 
*~ H(p) 


where H=V(p) is the operational equation of the applied e.m.f.: 
that is, analytically 


1 tet (aes ee 
‘vey= | E(t)e-”' dt. 


Now suppose that the impressed e.m.f. is sin wf: then by formula 
(h) of the table of integrals 


ja ga awed 
V@)=s3 403 (146) 
and denoting x by wz 
wp 1 
Fens en oe 147 
* p +o? H(p) a 


If, on the other hand, the impressed e.m.f. is cos wt, then by 
formula (i) 


V as 148 

and 
dee oe 149 
eS Fe tot HD) ae 


Now let us consider the operational expansion suggested by the 
Heaviside processes: 


eee OVE (2 Aa} fe 1 5 
Goal +(2) A herecanTiesy te 


and 


i Oma ee | | wey (151) 
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Now let us identify 1/H(p) with h(t) and replace p* by d"/dé": 
we get \ 
Mid siden ied 


Jie armenia ri ae 
and 
[id 14): 14 | 
= 1 de ota sae {7 a 


We have now to inquire into the significance of equations (152) 
and (153), derived from the operational equations for the response 
of the system to an e.m.f. sin wf and cos wt respectively, impressed 
at time ¢=0. From the mode of derivation of these expansions 
from the operational equations it might be inferred that they are 
the divergent or asymptotic expansions of the operational 
equations (147) and (149). This would certainly not be an 
unreasonable inference in the light of the Heaviside expansions 
we have just been considering. This inference is, however, not 
correct: on the other hand, the series (152) and (153) have a 
definite physical significance, as we shall now show from the 
explicit equations of the problem. 

By equation (31), the explicit equation for x,, given operation- 
ally by (147), is 


t t 
a -2/ sin wr.h(t—7)dr -{ sin w(t—7)h'(r)dr+h(0) sin wt (154) 
0 0 


where h’(t)=d/dt h(t). By a well known trigonometric formula, 
this is 


t t 
2,°= 61 ot | cos wt.h’(t)dt —cos ot | sin wt.h’(t)dé-+h(0)sin wt. 
0 0 


Writing 
t (or ro) 
{| dt= il a | dt 
0 0 t 


ie) ive) 
%,=sin ot { cos wt.h’(t)dt—ecos ot { sin wt.h’(t)dé 
0 0 


this becomes 


+h(0)sin ot— | sin w(t—7)h’(r)dr. (155) 
t 


The first three terms are simply the steady state response to 
the impressed e.m.f. sin wé: that is, they represent the ultimate 
steady state value of x, when the transient oscillations have 
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died away. The last term, which we shall denote by 7’, repre- 
sents the transient oscillations which are set up when the e.m.f. 
is applied. Thus 


iva) 


T.= ={ sin w(t—r)h'(r)dr. (156) 


Now from (156) 
= -1f h'(r).d.cos w(7 —#) 
and integrating by parts 
1d ie d? 

si = at +4 cos o(r—t) 7 h(r)dr. (157) 
Repeating the process of partial integration, we get: 

if, Ge {| Se oer (158) 

*  wdt ow Ji dr* , 
Repeating me ee eee 
T.=— 5 ot) <aahl)+= ale sin w(7 — ee Nee 


This process can be repeated indefinitely, and we get 


1 1 3 il 5 n—1 J2n—-1 
T.=( ee a 


wdt wdt?  w*dt® aye} 


—1)[". 
ue | sin w(r—)a ea h(r)dr. (159) 


The series expansion (159), except for the remainder term, is 
identical with the series expansion (152) derived directly from 
the operational equation. This series may be either convergent 
or divergent, depending on the frequency w/2m and the character 
of the indicial admittance function h(t). In the important 
problems of the building-up of alternating currents in cables 
and lines we shall see that, even when divergent, the series is 
of an asymptotic character and can be employed for computation. 

We thus arrive at the following theorem: 

If an e.m.f. sin wf is impressed at time 4=0 on a network or 
system of generalized indicial admittance h(é), and if the transient 
distortion, T,, is defined as the instantaneous difference between 
the actual response of the system and the steady state response, 
then 7, can be expressed as the series, 

eee Lae (—1)""'d 
(4 od tw 5dt5 a wen oo (160) 
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with a remainder term 

aby fe .. qamti 

ein } s1n ott) anti h(r)dr. 


If the impressed e.m.f. is cos wt, the corresponding series for the 
transient distortion, 7, is 


id? ld aled (—1)" a) 
es SN Ei ee 161 
& dt? w*dt4 ia dt® a2” dt2” h(t) ( 6 ) 
with a remainder term 
— 1)" Me 2nt+1 
"eo cos o(t 1) angi U(r) dr. 


The second part of this theorem, relating to the transient 
distortion, 7, in response to an e.m.f. cos wt, is derived from 
formula (31) by processes precisely analagous to those employed 
above in deriving the series expansion for 7. The derivation 
will be left to the reader. 

To summarize the preceding discussion of the divergent 
solution of operational equations, it may be said that the theory 
is as yet rather unsatisfactory. To the physicist it is unsatis- 
factory because he requires an automatic rule giving a correct 
asymptotic expansion by purely algebraic operations without 
investigations of remainder terms or auxiliary functions. Fur- 
thermore, the precise sense in which the expansion asymptotically 
represents the solution cannot be stated in general, but requires 
an independent investigation in the case of each individual 
problem. 

On the other hand, when an asymptotic expansion is known 
to exist, the Heaviside Rule finds this expansion with incom- 
parable directness and simplicity, the problem of justifying the 
expansion being a purely mathematical one, which usually need 
not trouble the physicist. Furthermore, on the purely mathe- 
matical side, the Heaviside Rule is of large interest and should 
lead to interesting developments in the theory of asymptotic 
expansions. 

APPENDIX TO CHAPTER V 

It will have been observed that the asymptotic expansions 
derivable from the operational equation by the direct applica- 
tion of the Heaviside Rules, as discussed in the preceding chapter, 
are all of the form: 


il 
BRB ar Sales Mags ais Sau) a> a) 
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In other words, they are all simple series in inverse fractional 
powers of ¢. In addition to this limitation, the direct applica- 
tion of the Heaviside Rules, as we have seen, does not always 
lead to correct results. In this appendix we shall briefly discuss 
some extensions and modifications of the direct Heaviside method 
whereby the limitations, inherent in the direct application of the 
Heaviside Rules, may often be successfully avoided. 

The asymptotic solution of operational equations is of great 
importance from the standpoint of both pure and applied mathe- 
matics. Its importance in practical applications to technical 
and physical problems does not require emphasis, while from the 
mathematical standpoint it not only is of value in the solution 


ie) 
of infinite integrals of the type j f(t)e-¥' dt, but promises to be a 
Jo 


useful tool in the general theory of asymptotic expansions. 
I. We start with the operational equation 


F(p) 
i= ry " il 
G(p) (1) 
and assume that F(p) admits of the formal series expansion 
F(p) =bo+b; ptbe p?+b; p?=- este (2) 


It is evident that equation (1) is a generalization of (110) of the 
preceding chapter, into which it degenerates when 1/G(p) is 
put equal to ~/p. A consideration of (1) will therefore lead to a 
formal generalization of Rule I. 

Let the function g =g(t) be defined by the operational equation 


g=1/G(p) (3) 
and an auxiliary function k=k(p) by the operational equation 
k=p°F(p). (4) 


Assuming that these functions exist, we have by Borel’s theorem 
t 

h@)= { k(r)g(¢—71)dr (5) 
J0 


Now returning to (1) let us write, by virtue of the expansion 


(2), 
h=(aotaiptarp?+:: ls (6) 
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Now, guided by the usual Heaviside procedure, let us replace 


1/G(p) by g(é), in accordance with (3), and replace p” by d*/dt" 
in the series factor. The right hand side of (6) then becomes 


d d? - 
(dota 7 tae Get -++)g(@) (7) 
We have now to inquire into the significance of this expansion. 
Expand the function g({—7) as a power series in 7, and substi- 
tute in (5): we get 
t t 
h(t)=g(t) | k(r)dr—g'(t) | Gk(r)dr 
F ot! 
ton 
$e (=1)9gQ) (i 7" k(a)dr (8) 
0 . 


t 
+f k(r)R, (¢—7)dr 
where F,,(¢—7) is the remainder in the expansion of g(¢—r). 


We now suppose that k(¢) converges in such a way that the 
integrals 


|, “k(r)dr — (n=1,2,3. .1) 
o 2! 
exist. On this assumption it follows from the integral equation 
i) k(t)e-?' di =F (p) 
0 


and the series expansion (2), that 


|| 7 k(@)dr = (—1)"bp (9) 


Ihelhll 


equation (8) becomes by virtue of (9), 


Now writing 


id= (dot n§, Eh fa dt? ai “+0, asa 


+ | R\Ra(t—a)dr (10) 


=| dr Bo) tot -+(- 1 Fae (9 
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The first series of (10) is identical with (7), while the rest of 
the right hand side of (10) represent remainder terms. The 
question as to whether the series (7) represents the function 
h(t), asymptotically or otherwise, can only be answered by an 
independent investigation of the remainder terms in each specific 
problem. The expansion (10) therefore formulates a general- 
ization of the problem underlying Rule I, into which it degener- 
ates when g(t)=1/+/zt. The fact that no general statement 
can be made as to the way in which the series (7) represents 
h(t) indicates the care which must be taken in the interpretation 
of the operators. Often, however, the expansion (7) is valid 
and quite useful: in addition it possesses considerable theoretical 
interest. 

II. If the variable z(t) is defined by the definite integral 


x(t) =< if ‘OO n(n)dr (2.1) 


it is known that, provided the real part of \ is not negative, it 
may be written as 


2) = FR) tv (2.2) 


Furthermore, in a number of important cases the component 
y(t) admits of asymptotic expansion. On the other hand, the 
operational equation 
D 1 
~ p—-hH(p) oe 
corresponding to (2.1) does not admit of direct asymptotic 
expansion and the direct application of the Heaviside Rule 
is invalid. We can, however, set up an operational equation 
for y(t) which leads to the required asymptotic expansion, where 
such an expansion exists. The procedure is as follows: 
From (2.2) and (2.3) we have the usual integral equation: 


I. dte-? | vO a On’ ls HH) 


Integrating the second term in the integral and rearranging 


gives 
io) a 1 1 1 
{ die” y(t) = p— 1 I H(p) H(A) | 


and the corresponding operational equation is 


ae : i 7 1 
¥= 5 | We) ae 2) 
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As an example, let \=dw and h(t)=1/+/zt. The problem is 
then that discussed on page 66 and formulated by equation 
(121), concerning which it was remarked that the straight- 
forward procedure gave incorrect results. In this case 

1/H(p)=Vp; — 1/HOA)=Viw 
and the operational equation is 


y=—?_|Vp- Viel 


Pp 


eee ae 
V/V p+ V iw V iw i fe 
or, expanding in powers of ~/p, 
ee eos Dig ee t) 
Vie tw Ven (iw) (iw) V iw 
Dropping the integral powers of p, according to Rule II (p07) 
and replacing ~/p by 1/+/zt and p” by d"/dt" gives 


.j1d AL oh ha 
Yu adt ot dB [a 
1d “lead 
eee ope Ce 


which is correct. 

Note that, if in the last equation we generalize 1/+/zt into 
h(t), the series becomes the equivalent of (160) and (161) for the 
transient distortion in the case of sinusoidal impressed forces. 

III. In a number of important problems, the function h, 
defined by the operational equation h=1/H(p), has an asymp- 
totic expansion of the form 


AM 
h~ (tort ar/ttas/P+- ; -) (3.1) 


In such cases the direct application of the Heaviside rules 
breaks down, due to the presence of the exponential factor. If, 
however, the form of the asymptotic expansion is known, the 
following procedure enables us to derive the asymptotic expansion 
by operational methods. 
Let us write 

h(t) =e™ g(t) (3.2) 

and substitute in the integral equation: we get 


See hee 
ie dt e~?' eM g(t) = pH) 
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and, by Theorem Va, 


—pt — ted 
{ Oe ENE (ph oe 


whence, operationally, 


a (es : 
9 DENA +) oe) 


If the assumed form of expansion exists, then g must admit of 
expansion in inverse fractional powers of ¢, and this expansion 
must be derivable from the operational equation by the direct 
application of the Heaviside rules. 
Example 

As an interesting example of the preceding we shall derive 
the usual asymptotic expansion of the Bessel function Jo(dt) 
from the operational equation (n) and the corresponding integral 


equation 
uf 


= pt = A 

ii dt e~?* J (Nt) ory 
It is known that the function J9(A¢) can be represented as 
u(t) cos t+ v(t) sin At (3.6) 
and that uw and v admit of asymptotic expansions in inverse 


fractional powers of t. 
Write (3.6) as 


x Ono” sem u(t) +20(2)] (3.7) 


(3.5) 


and substitute in (3.5): we get 
: i dt e-Pte™ [u(t) —év(t)] (3.8) 
0 


ee pie | 
7 ib POE SAG ION a rice 


whence, replacing p—7A by q and then replacing qg by p, we get: 
; i di e-?*[u(t) —iv(t)] (3.9) 
0 


+ An dt e~Pt ee [u(t) +iv(t)]= Waiae is 


Now the function 

e~™ Tu(t) +iv(6)] 
is oscillatory and does not admit of asymptotic expansion in 
inverse fractional powers of ¢. Consequently if the right hand 


84 ELECTRIC CIRCUIT THEORY 


side of (3.9) is expanded in fractional powers of p, we shall get 
therefrom by the procedure discussed in Chapter V the asymp- 
totic expansion of the function 


1 ; 
p lu(t) —in(2)] 
That is to say, in an asymptotic sense only, we can set up the 
operational equation 
; 2p 
1.1 = 210 
Vp? +i2rp ( ) 


2 2 iD. —}4 
= ox | Vp 


The factor e~'"/* suggests that u —iv be replaced by e—'"/*(P —7Q), 
whence 


p~ig= P| 1-12 |v (3.11) 


Expanding the right hand side by the Binomial Theorem, replac- 
ing p” by d”/dé” and /p by 1/+/xt, we get: 


Pi hi eee 
me N ay 8vt_ 2!(8dt)? 
12.32.52 12.32.52.72 


TT Ee ABE | Cie 
whence 
2 12.32 
Pry) pp | 
ner (3.13) 
on lee Lael Be 
U~a/7| Bat 31(8d8) | 
and finally 
Jo(dt)~P cos (w-7)+¢ sin( =) (3.14) 


which is the usual asymptotic expansion. 

The foregoing discussion is tentative and far from complete. 
It does, however, suggest that important developments of the 
theory of asymptotic expansion may result from a study of the 
operational equation or rather its integral equation equivalent. 


CHAPTER VI 


PROPAGATION OF CURRENT AND VOLTAGE IN THE 
NON-INDUCTIVE CABLE 


The principal practical applications of the operational calculus 
in electrotechnics are to the theory of the propagation of current 
and voltage along transmission systems. Of such transmission 
systems the simplest is the non-inductive cable. The theory of 
the non-inductive cable is not only of great historic interest, 
relating as it does to Kelvin’s early work on the possibility of 
transatlantic telegraphy, but is also of very considerable practical 
importance today, and serves as a basis for the theory of submarine 
telegraphy over long distances. We shall therefore consider the 
propagation phenomena in the non-inductive cable in some detail. 

The propagation phenomena in any type of transmission 
system are isolated and exhibited in the clearest possible manner 
when we confine attention to the infinitely long line, with voltage 
applied directly to the line terminals. Furthermore, as we 
shall see later, the solution for the infinitely long line is funda- 
mental and can be extended to the more practical case of the 
finite line with terminal impedances. We therefore, in this 
chapter, shall confine our attention to the case of the infinitely 
long cable with voltage applied directly to the cable terminals. 

Consider a cable of distributed resistance R and capacity C 
per unit length, extending from x=0 along the positive x axis. 
From a previous chapter (see equations (64) and (65), we are in 
possession of the operational equations of voltage and current; 
they are, for the infinitely long line, 


T= Jap e-VeV =| Pave Vo (163) 
Vides : R , 


where a=2?RC, and Vo is the terminal cable voltage at x=0. 
Let us now assume that the terminal voltage Vo is a ‘unit 
e.m.f.;” then = 
VeeV2, (164) 
aes yp e-~V ap 
Hi = BV we : (165) 
85 
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The solution of (164) for V was considered in some detail in the 
preceding chapter; it is, by (129) 


Vets | ae 166 
avd ev on 


where 7 =4t/a=4t/x?RC. Series expansions of this solution were 
also given. Another equivalent form is, by (131) 


2 (1/7 
va1- Z| eT dr. (167) 
0 


This last form, recognizable also from inspection of the series 
expansion (132), is useful because the integral term is what is 
called the error function and has been completely computed and 
tabulated. 

Before discussing these formulas and the hight they throw on 
propagation phenomena in the non-inductive cable, we shall 
derive the solution for the current. A very simple way of doing 
this is to make use of the differential equation (57) 


ole 8 
Rox 
Now from (166) and the relation 


oo _drd 
ax dxdr 


[= Ke 


we get : 
Oe ee 
Ox Va rv/7 da vRC 


whence 


I = —— 2 Daa G =, 
tRV a2 Vr By otk ae elo), 
It is worthwhile verifying the formula by direct solution from 


the operational equation (165). From formula (g) of the table 
of integrals, we have 


h=e2V* veo 


pane (6: 
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Comparison with the operational equation shows that they are 
identical, within a constant factor provided we put \=a/4. 
Consequently the solution of (165) is 
ce coh 
= Ee Bay eS Be 5/7 
: rRi” wR 

which agrees with (168). This, it may be remarked, is an excel- 
lent example of the utility of the table of integrals in solving 
operational equations. 

This formula is easily calculated for large values of t by 
expanding the exponential function; it is 


Edep-C40)--] 


The propagation phenomena of the non-inductive cable are 
therefore determined by the pair of equations 


tp—1/r 1/vV7 
ee ie dra1— | e-"'dr (169) 
Vi Jo 


ne 
and 
2 ea l/r vis : 
pe eS el a ee a 170 
i VreR Vr eRe uo) 


At 
where T= 4t/a= apa 


Now an important feature of these formulas is that the voltage 
: 4 : She 
at point x is a function only of 2RC t; that is, of 4¢ divided by the 


total resistance and capacity of the cable from x=0 tox=az. The 
same statement holds for the form of the current wave: its mag- 
nitude, however, is inversely proportional to xR, or the total 
resistance of the cable up to point x. Consequently a single 
curve, with proper time scale serves to give the voltage wave at 
any point on the cable. Similarly a single curve, with proper 
time and amplitude scales, serves to depict the current wave at 
any distance from the cable terminals. These curves are given 
in Figs. 3 and 4. 

Referring to the curve depicting the current wave, we observe 
that it is finite for all values of t>0; consequently, in the ideal 
cable, the velocity of propagation is infinite. This is a con- 
sequence, of course, of the fact that the distributed inductance 
of the cable is neglected. Actually, of course, the velocity of 
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propagation cannot exceed the velocity of light. The error, 
however, in neglecting the inductance in the case of long cables 
is appreciable only near the head of the wave provided we confine 
attention to d.c. or low frequency voltages. This point will be 
discussed and explained more fully in connection with the 
transmission line. 


fics aa 


| 
a 


o 


T=4tAFRC 
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Seen ee 
T= 4t/PRC aes 
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Fie. 4.—Voltage in non-inductive cable (G =0) unit e.m.f. applied. 


The current, while finite, is negligibly small until r reaches the 
value 0.2. In the neighborhood of this point it begins to build 
up rapidly; reaches at 7=2 its maximum value 

2 ¢-05 2 
= = —=—— (0.429 
nein eA eras | 
and then begins to decrease, ultimately dying away in accordance 
with the formula 


Senator} 
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Its subsidence to its final zero value is very slow; for example, 
when 7 = 100 its value is still 


Turning to the voltage curve, Fig. 4, we see that it is negligibly 
small until 7 reaches the value 0.25, at which point it begins to 
build up. Its maximum rate of building up occurs when 7 = 2/3, 


Fre. 5.—Power transmitted in non-inductive cable (4 =0). 


after which it builds up more and more slowly. Its approach to 
its final steady value is in accordance with the formula 


2 Lee lar! 
Sy foe ea ee ee 
Even, therefore, when 7 is as great as 100, V differs sensibly 
from its ultimate value, unity, its value being 0.8876. 


: : eee IAG” S. 
Since the actual time is Saas it follows that the speed of 


building up is inversely proportional to the square of the length 
of the cable. 

The power curve V-I is given in Fig. 5. V-I is the rate at 
which energy is being transmitted past the point x of the cable. 

The fact that the form of the current and voltage waves 
depends only on 4t/xz? RC is at the basis of Kelvin’s famous “KR” 
law, long applied to cable telegraphy and sometimes incorrectly 
applied to telephony. When the first transatlantic telegraph 
cable was under consideration, Kelvin attacked the problem of 
propagation along the non-inductive cable and arrived at for- 
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mulas equivalent to (169) and (170). From these formulas he 
announced the law that the “‘speed”’ of the cable, i.e., the number 
of signals transmissible per unit time, is inversely proportional 
to the product of the total capacity and total resistance of the 
cable (KR in the English notation). To see just what this 
means requires a little digression into the elementary theory of 
telegraph transmission. 

Telegraph signals are transmitted in code by means of ‘‘dots”’ 
and ‘‘dashes.”’ The ‘‘dot’’ is the signal which results when a 
battery is impressed on the cable for a definite interval of time, 
after which the cable is short circuited. A ‘‘dash”’ is the same 
except that the time interval during which the battery is con- 
nected to the cable is increased. The “dots” and “dashes”’ are 
separated by intervals, called ‘‘spaces,’”’ during which the cable 
is short circuited. Now when the cable is short-circuited we may 
imagine a negative battery impressed on the cable in series with 
the original battery. Consequently the current in the cable, 
corresponding to a signal composed of a series of dots, dashes and 
spaces, will be represented by a series of the form 

I(t) —I(t—t)) +] (t-te) —T(t—ts) +1 (t-t4) —- - - (171) 
where, in the cable under consideration, /(t) is given by (168). 
t; is the duration of the first impulse, f:—t, of the first space, 
t3—t, of the second impulse, etc. 

Now by (168) 


Bec 
tO 


tR~/r Ve sani soi): 


rT is, of course, 4t/x?CR=4t/KR (in the English notation). Now 
suppose that 


4h 
3 ea 
T2= oop’ etc. 
Then the Bae can be written as 
ae {$(r)— o(r—71) + o(7—72) — + + +}. (172) 
Now if the relative time intervals 71, 72 . . . are kept constant 
(as the length of the cable is varied), the actual time intervals 
ti, lg . . . are proportional to z*CR or to KR, and the wave 


form of the total signal is independent of KR, when referred to 
the relative time scale 7. Hence, if 7 is the total time of the 
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signal, T is proportional to z?CR (or to KR). That is to say, if 
the duration of the component dots, dashes, spaces of the signal 
are proportional to the “KR” of the cable, the wave form of the 
received signal, referred to the 7 time scale, is invariable, and the 
total time required to transmit the signal is proportional to 
the “KR” of the cable. Now the maximum theoretical speed of 
transmission on the cable is limited by the requirement that the 
received signal shall bear a recognizable likeness to the original 
system of dots and dashes: in other words there is a maximum 
allowable departure in wave form between received and trans- 
mitted signals. If, therefore, the actual speeds of two cables 


Fia. 6.—Elementary telegraph signals in non-inductive cable. 


are inversely proportional to their ‘‘KRs,”’ the wave form will be 
the same. ‘This establishes Kelvin’s “KR” law. Asa corollary, 
if the length of the cable is doubled the speed of signaling is 
reduced to one-quarter, assuming the same definition of signals. 

The foregoing will be somewhat clearer, perhaps, if we refer to 
curves, 4, 5, 6, 7 of Fig. 6 which illustrate the distortion suffered 
by elementary dot signals in cable transmission. Curve 4 shows 
the dot signal produced by a unit battery applied to the cable 


2 
terminals for a time interval {=2 oe while curves 5, 6 and 7 
are the corresponding dot signals when the battery is applied 


; : CRC Lene 
for the time intervals Ta A and i ae Any further 


decrease in the duration of the impressed dot, beyond that 
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shown in curve 7, does not affect the shape of the transmitted dot, 
which means that the cable speed has reached its theoretical 
maximum. These curves, it should be observed, can be inter- 
preted in two ways. First, we can regard the length x of the 
cable as fixed and the duration of the impressed dot as varied. 
On the other hand, we can regard the actual duration of the 
impressed dot as constant and the length of the cable as varied. 
From the latter standpoint the curves illustrate the progressive 
distortion of the signal as it is transmitted along the cable. 
The dot signal of relative duration 7 can be written as 
D=I(r), ee li 
=I(r)—I(r—T), tor 

and the second expression can be expanded in a Taylor’s series, 
giving 


d fe hg 
D=T IL) — a gal@t: a 
If T is sufficiently short this becomes 
D=T'I'(r). (173) 


Hence when the dot signal is of sufficiently short relative duration 
T, the wave shape of the received signal is constant, J’(r), and 
its amplitude is proportional to the relative duration of the dot. 

This can be generalized for any type of transmission system: 
Let the dot signal be produced by an e.m.f.f(t) of actual duration 
T. Then the received dot signal, by formula (31), is 


a re 
Dai | otter, Ea 


ff 
all f(r) L[(t—7)dr, tr 


For t>T this becomes 


The TF 
=r | feeddr—2" |. th(r)dr+--- 


and for sufficiently short duration 7’, we have approximately, 
T 

p=. S@r)ar. (174) 
0 


Hence for a sufficiently short duration of the impressed e.m.f. 
the received dot signal is of constant wave form, independent of 
the shape of the impressed e.m.f., and its amplitude is propor- 
tional to the time integral of the impressed e.m.f. These princi- 
ples are of considerable practical importance in telegraphy. 
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The leaky cable, that is, a cable with distributed leakage 
conductance G in addition to resistance R and capacity C, is of 
some interest. The differential equations of the problem are 
given in equations (70); the operational formulas for the case of 
voltage directly impressed on the terminals of the infinitely long 
line are 

Vm e~*VORDFRG V, 
pO G o-*V ORp ERG 
7 +p° VCRp+RG VY, 
Writing CRxz?=a and RGx?=8, G/C=, and assuming a “unit 
e.m.f.’’ impressed on the cable, this becomes 


V =e- Vents, (175) 
Taa[ovpthe-Vorrs (176) 


These equations are readily solved by means of the table and 
formulas given in a preceding chapter. 

But first let us attempt to solve the operational equation (175) 
for the voltage by Heaviside methods, guided by the solution of 
the operational equation 


jf 


V =e- Ven (124) 


of the preceding chapter. Expand the exponential function in 
(175) in the usual power series; it is 


1—Vop p+ CEth)_ eee ee 2 Kebep 


Now discard the integral terms and write 


fas i | Ven 8. (178) 


5! 


We have now to interpret the expression \/ap+f8. We have by 
ordinary algebra 


ne Bt ok B 1/2 orat nN 1/2 a. 3 
vapta=(1+5) Vap=(1+>) V ap 
LIS Lone = 
=|dteae aie as (2) aa |vao. 


Now identify ~/p with 1/+/zt in accordance with the Heaviside 
rule, and 1/p with dt. We get 


Vap+b= NE f14+¥—° a +540) s- a} (180) 


(179) 
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Now in the terms of the expansion (178) identify p” with d"/dt" 
and substitute (180); we get 
try 
(181) 


y/o d 
Yat | tts (agte) tg Grr ae 
Ae MO of 


This series is hopelessly complicated to either interpret or com- 
pute. It is, in fact, an excellent illustration of the grave dis- 
advantages under which many of Heaviside’s series solutions 
labor. We shall therefore attack the solution by aid of the 
theorems and formulas of Chapter IV. The simplicity of the 
solution which results is remarkable. 
The operational formula for the voltage is 
V =e- Var +8, (175) 
Now the operational formula for the voltage in the non-leaky 
cable is (see equation (164) ) 
V =e- Vee, 


In order to distinguish between the two cases, let us denote the 
voltage in the latter case by V®; thus 
V9 =e- Van, (182) 
Now by theorem (Va) and equation (182) we have 
ogp—-M=_P_ 4—-Valp tn), 
V-e E ay é P 


(183) 
PP o-Vapts. 
ae 
Now write (175) as 
= PER, Peo Vents, 
: re (184) 
[list a, rae e~ Vapté, 


It follows at once by comparison with (183) and the rule that 
1/p is to be replaced by f dé, that 


t 
ae +n iv) Pope. (185) 
0 


By a precisely similar procedure with the operational formula 
(176) for the current, we get 


ey 
= (14a au)Ie™ (186) 
0 
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where I° is the current in the non-leaky cable. Now by formulas 


(169) and (170) 
Vee 1 acme, (169) 
Ve i Va 


wie 
v= _~ _p-a/4t 
I pee : (170) 


which completes the formal solution of the problem. 

Formulas (185) and (186) are extremely interesting, first as 
showing the superiority of the definite integral to the series 
expansion—compare (185) with the series expansions (181)— 
and secondly as exhibiting clearly the effect of leakage on the 
propagated waves of current and voltage. We see that in both 
the current and voltage the effect of leakage is two-fold; first it 
attenuates the wave by the factor e~-™“, (\=G@/C), and secondly it 
adds a component consisting of the progressive integral of the 
attenuated wave. This, it may be remarked, is the general 
effect of leakage in all types of transmission systems. Its effect 
is, therefore, easily computed and interpreted. 

Formulas (185) and (186) are very easy to compute with the 
aid of a planimeter or integraph; or, failing these devices, by 
numerical integration. However, for large values of t, the 
character of the waves is more clearly exhibited if we make use 


of the identity 
t oo (e) 
fia={ a dt 
0 0 t 
v=(inf iv) v9e%—y f Ve “dt (187) 
0 3 
r=(4f i)Ire™ xf Ie “ak. (188) 
0 t 


The first two terms of these formulas are clearly the ultimate 
steady state values of the voltage and current waves, and can 
be determined by evaluating the infinite integrals. A far simpler 
and more direct way, however, is to make use of the fact that the 
ultimate steady values of V and J are gotten from the operational 
formulas by setting p=0. That this statement is true is easily 
seen if we reflect that the steady d.c. voltage and current are 
gotten from the original differential equations of the problem by 
assuming a steady state and setting d/dt=0. 


whence 


and 
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From the operational formulas we get, therefore 
(i4af dt)VPe-Mae~ Vim e-*V, (189) 
0 
z Cn ne _ WRG 
On=i— =V/ Ge. | oa RG 
(14 au)ire Re Re ; (190) 


Introducing these expressions into (187) and (188) respectively, 
we get 


Vaertvno—y | V°e-dt, (191) 
t 
Ga /RG Pitons 
HN est ae PT ati (192) 
t 


The definite integrals can be expanded by partial integration; 


thus 
=u) Vee-atm | Vode-™ 
t t 


ea 0 : ee 70 
V°e I; e dt dt. 


Continuing this process we get 


= p—tW RG px , 0 
V=e e (tart yrs (193) 
G 


T= [Ge-2V 0 em(it4 aE ean re (194) 


Using the values of V? and J°, as given by (169) and (170), it is 
extremely easy to compute V and J, for large values of t, from 
(193) and (194). 

So far we have considered the current and voltage waves in 
response to a “unit e.m.f.,’”’ impressed on the cable at x=0. 
It is of interest and importance to examine the waves due to 
sinusoidal e.m.fs., suddenly impressed on the cable, particularly 
in view of proposals to employ alternating currents in cable 
telegraphy. 

We start with the fundamental formula 


z(t) = aC —rt)h(r)dr 


t 
= ‘i f(t—r)h (r) dr, 
0 
provided h(0) =0, which is the case in the cable. 
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If f(£) =sin wt, we write 


t 
aia ot | cos wt.h’ (t)dt 
0 


(194-a) 
—cos wt { sin wt.h’(t)dt. 
0 
Similarly, if the impressed e.m.f. is cos wt, 
t 
X(t) =cos ot | cos wt.h’(t)dt 
: (194-b) 


t 
+sin ot sin wt.h’(t)de. 
0 


The investigation of the building-up of alternating currents and 
voltages, therefore, depends on the progressive integrals 


t 
C =|| cos wt.h’(t)dt, 
0 
: (194-c) 
s=| sin wt.h'(t)dt. 
0 


For the case of the voltage wave on the non-inductive, non-leaky 
cable these integrals, by aid of equations (169), become, if we 


write w’ =aw/4, 
1 Te-1/ cos w'7 
C dr, 
0 


Te rW/1 


Ge 1 { e-it sin o'r (194-d) 
Vaso INT 
where, as before, 7 =4t/a. 
For the current wave we have, by (170), 
2 T/1  1\ e-*/ cos o'r 
al Ge ) Cera 
GE 2 
V4 xR 0 tr (194-e) 


pe (Lal) eo ain or 
Vrtk Jo \t 2 tr 


For small values of 7 and w’ these integrals can be numerically 
evaluated without great labor. Mechanical devices, such as the 
Coradi Harmonic Analyzer, are here of great assistance. In fact 
the Coradi Analyzer gives these progressive integrals auto- 
matically. It may be said, therefore, that a complete mathe- 
matical investigation of the building-up of alternating current 
and voltage waves in the non-inductive cable presents no serious 
difficulties, although the labor of computation is necessarily 
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considerable. One fact makes the complete investigation 
much less laborious than might be supposed. This is, if the 
foregoing integrals are calculated for a given value of w’, the 
results apply to all lengths of cable and all actual frequencies 
w/2r, such that aw is a constant. Thus if we double the length 
of the cable and quarter the frequency, the integrals are 
unaffected. 

The solid curve of Fig. 7 shows the building-up of the cable 
voltage in response to an e.m.f. cos wt, impressed at time ¢=0. 
The frequency w/2z7 is so chosen that w’ = aw/4 = 27, and the curve 
is calculated from equations (194-b) and (194-e). The dotted 


0 t=4t/x2?RC 


=035 
Fia. 7.—Non-inductive cable (4 =0), building-up of alternating current. 


Applied e.m.f. cos wt; w = 27 ane 

curve shows the corresponding steady state voltage on the cable; 
that is, the voltage which would exist if the e.m.f. cos wt had 
been applied at a long time preceding t=0. We observe that, 
for this frequency, the building-up is effectually accomplished 
in about one cycle, and that the transient distortion is only 
appreciable during the first half-cycle. 

The case is very much different when a higher frequency is 
applied. Fig. 8 shows the building-up of the alternating current 
in the cable when an e.m.f. sin wt is applied at time #=0. The 
frequency is so chosen that w’=aw/4=107. The outstanding 
features of this curve are that the initial current surge is very 
large compared with the final steady state, and that the transient 


CURRENT AND VOLTAGE IN NON-INDUCTIVE CABLE 99 


distortion is relatively very large. It is evident that the fre- 
quency here shown could not be employed for signaling purposes. 
This curve has been computed from the steady state formulas, 
and equations (160) and (161) for the transient distortion. 

If the applied frequency w/2z is very high, the steady state 
becomes negligibly small, and the complete current is obtained 
to a good approximation by taking the leading terms of (160) 


een 
oe ae 
ome 
Ea 
NAA AAD AD pp | 


I 2.0 3.0 40 
Fig. 8.—Non-inductive cable (G=0). Building-up of alternating current. 


. ; ley es oe. 
Applied e.m.f. sin wt; w = 107 2RC 


and (161). Thusif the applied e.m_f. is sin wt, and w is sufficiently 
large, the cable current is 
Dein GaeC 


Vaahe! dr Jz 
by (160) and (170) while, if the impressed e.m.f. is cos wé, it is 
2 1\2d? e-1/ 
by (161) and (170). Here w’=aw/4 and r=4t/a. 


CHAPTER VII 


THE PROPAGATION OF CURRENT AND VOLTAGE IN 
THE TRANSMISSION LINE 


We now take up the more important and difficult problem of 
investigating the propagation phenomena in the transmission line. 
The transmission line has distributed series resistance R and 
inductance L, and distributed shunt capacity C and leakage con- 
ductance G. It is the addition of the series inductance L which 
makes our problem more difficult and at the same time introduces 
the phenomena of true propagation with finite velocity, as 
distinguished from the diffusion phenomena of the cable problem. 
The cable theory serves very well for the problems of trans- 
oceanic telegraphy! but is quite inadequate in the problems of 
telephonic transmission. 

If J denotes the current and V the voltage at point x on the line, 
the well known differential equations of the problem are: 


d ous 
(Lgtk)I=— ae (195) 


d Fee 
Replacing d/dt by p, these become 


(Lp+R) I= -2Y, 
(196) 
0 


From the second of these equations 


ON ae eo 
Ox Co+G dx? 


1 With the installation of the new submarine cable, continuously loaded 
with permalloy, this statement must be modified. In this cable, the induc- 
tance plays a very important part, and is responsible for the greatly increased 
speed of signaling obtainable. (See paper by Buckley, A. I. E. E., Aug., 
1925, ‘The Loaded Submarine Telegraph Cable.’’) 
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and substitution in the first gives 


(Lp +R)(Cpt+@)= 241. (197) 
Similarly if we eliminate J, we get 
2 
(ip eR Cpe OV = Ss V. (198) 


If we assume a solution of the form 
Vis Ace Bek 


where A and B are arbitrary constants, substitution shows that 
the solution satisfies the differential equation for V provided 


y?=(Lp+R)(Cp+@). (199) 
From equation (196) it then follows that 
[= iE ,(Ae77® — Bev") (200) 
SoU 


“(Ae771* — Ber), 
VY 


Now restricting attention to the infinitely long line extending 
along the positive x axis, with voltage V» impressed directly on 
the line at x=0, the reflected wave vanishes and we get 

V=Vo Fauee 


Ves, (201) 


y?=(Lp+R)(Cp+G). 
Now let us write 


1 
v= [(p+p)*—o7] (202) 
where . 
v=1/V/LC, 
Ree 
P=ort 90” 
pe ns 
SOF. OE! 


Then setting Vo=1, the operational equations of the problem 
become 


Vas 1 V (we)? 0?, (203) 


CV go Vereme 
1=1( 0+ ) SSS is (204) 


P 
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Now consider the operational equation, defining a new variable F: 


~ FV OF 


(205) 


It follows at once from our operational rules, and (203) and (204), 


that 
t 
f= i(c +6) au)P, (206) 


va-of GM (207) 


Our problem is thus reduced to evaluating the function F, from 
the operational equation (205). This equation can be solved by 
aid of the operational rules and formulas already given. The 
process is rather complicated, and there is less chance of error if 
we deal instead with the integral equation of the problem 


a Sea ={ “P (t)ePt dt. (208) 
ae +p)? 2g? 0 
Now let us search through our table of definite integrals. We 
do not find this integral equation as it stands, but we do observe 
that formula (m) resembles it, and this resemblance suggests 
that formula (m) can be suitably transformed to give the solution 
of (208). We therefore start with the formula 


ee -{ T(E MDa (m) 


This, regarded as an integral equation, defines a function which is 

zero for é<) and has the value Jo(+/t?—X?) for £>X, Jo being the 

Bessel function of order zero. We now transform (m) as follows: 
(1) Let \p=q and ¢/A=¢,. Substituting in (m) we get 


eV ate : 
Nea +)? =[e at Jy (OV t? —1)dt,. 


Now, in order to keep our original notation in p and ¢, replace 
q to p and ¢; by t; we get 
e~V p+ 


at AV/t2—1)dt. m.1 
a | ( ) (m.1) 


(2) In (m.1) make the substitution p=q-+, and then in the 
final expression replace g by p; we get 


e- V (pF u)?+™ 


[ere OVE Dat 


aD 
Viceun ravine 
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(8) In (m.2) make the substitution p=—q and f==t, and 


ultimately replace q by p and ¢, by t; we get 


x 


- a) given 
{ EP! emt To Ve = at = ———— (m.3) 
x / v V (p+ u?+re 


v ; 
where A; =e Nand py =" wu. (They are, of course, as yet, arbitrary 


parameters, except that they are restricted to positive values.) 
(4) Now if we compare (m. 3) with the integral equation (208) 
for F, we see that they are identical provided we set 
MiP, ee 
\1=te =o —1. 


Introducing these relations, we have 


i. V Par ACE 
e?t ert I (c t?— 2? /y?2) dt = G (m.4) 
2/2 V (p+p)?—3" 
Here J, denotes the Bessel function of imaginary argument; thus 
J (iz) = J(2). 
It follows from (m.4) and the integral equation (208) that 
FG) =0 tor t<2/», a (209) 
=e" Io(ov/t?—2?/v?) for t>2/v. 
Having now solved for F =F (é), the current and voltage are gotten 
from equations (206) and (207). Thus 
T= ort <2/9, (210) 


a ‘ 
= JSro +f F(t)dt fort2=2/v. 


The corresponding voltage formula is 
V =0 fort<z/v, 


gong a8 (DOB) 
OS fy V7? 27/v? 

Here I,(0\/r?—.2?/v2) is the Bessel function of order 1: thus 
—iJ, (iz)=I,(z). The function is entirely real. The derivation 
of formula (211) is a little troublesome, owing to the discontinuous 
character of the function Ff’; the detailed steps are given in an 


appendix to this chapter. 
The preceding solution depends for its outstanding directness 
and simplicity on the recognition of the infinite integral identity 


drfort>zx/v. (211) 
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(m), into which the integral equation of the problem can be trans- 
formed. Whensuch identities are known their value in connection 
with the solution of operational equations requires no emphasis. 
On the other hand, we cannot always expect to find such an iden- 
tity in the case of every operational equation; and, particularly 
in such an important case as the transmission equation it would 
be unfortunate to have no alternative mode of solution. For- 
tunately a quite direct series expansion solution is obtainable from 
the operational equation, and this will now be derived. As a 
matter of convenience we shall restrict the derivation to the 
voltage formula 


Vea oe (203) 
As a further matter of mere convenience we shall assume that 
G= 0, so that =p and (203) becomes 

V =e7V2?+20p (203-a) 
where 7=2/v. 


The method holds equally well for the current equation (204) and 
for the general case o¥p. 
Write (203-a) as 
eS ent? (1+2p/p)1/2 


and expand the exponential factor (1+2p/p)!/? by the binomial 


theorem; thus 
3 
(1+2p/p)!2=1+! prea 8 a +an(#) oe 


so that 


a2Tp” _asTp® aaTp* 
Pp p? p° 


V =e?.e-. exp fe eS eg 


Now the operational equation 


Fees . 
p p° Pp 
can be expanded in inverse powers of p; thus 
TOR aeeegre 
the power series solution of ae is 
v=14ey peat +e eee 
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It follows at once from the preceding and Theorem VII that 
V(t) =0 fort <7 
=i +6 


sah eee 1 : for t>7. 


If the coefficients 81, 8. . . . are evaluated, a simple matter of 
elementary algebra, the foregoing expansion in the retarded time 
t—7 will be found to agree with the solution (211) when a is put 
equal to p. 

We shall now discuss the outstanding features of the propaga- 
tion phenomena in the light of equations (210) and (211) for the 
current and voltage. We observe, first, that we have a true 
finite velocity of propagation v=1/~/LC. No matter what the 
form of impressed e.m.f. at the beginning of the line (x=0), 
its effect does not reach the point x of the line until a time t=2/v 
has elapsed. Consequently v=«/t is the velocity with which the 
wave is propagated. This is a strict consequence of the dis- 
tributed inductance and capacity of the line and depends only on 
them, since v=1/+/LC. It will be recalled that in the case of the 
cable, where the inductance is ignored, no finite velocity of propa- 
gation exists. 

The question of velocity of propagation of the wave has been 
the subject of considerable confusion and misinterpretation when 
dealing with the steady state phenomena. It seems worth while 
to briefly touch on this in passing. 

As has been pointed out in preceding chapters, the symbolic or 
complex steady state formula is gotten from the operational 
equation by replacing the symbol p by zw where 7=~+/—1 and 
w/2mr is the frequency. If this is done in the operational equation 
(203) for the voltage, the symbolic formula is 

Vine vV Gate —o piot 
If the expression ~/(iw+p)?—o? is separated into its real and 
imaginary parts we get an expression of the form 


5 zx 
V =e gio(t a) . 
where 


ents aR Cir ine 
NS, se . 2w? 
and 


a= p/Bo. 
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Now if we keep the expression ¢ — 8 constant, that is, if we move 


along the line with velocity dx/dt =v/8, the phase of the wave will 
remain constant. This is interpreted often as meaning that the 
velocity of propagation of the wave is v/8. Now since # is 
greater than unity and only approaches unity as the frequency 
becomes indefinitely great, the inference is frequently made that 
the velocity of propagation depends upon and increases to a 
limiting value v, with the frequency. This velocity, however, is 
not the true velocity of propagation, which is always v, but is the 
velocity of phase propagation in the steady state. This distinction 
is quite important and failure to bear it in mind has led to serious 
mistakes. 

Returning to equation (210) and (211) we see that after a time 
interval t=2/v has elapsed since the unit e.m.f. was impressed on 
the cable, the voltage at point x suddenly jumps from zero to the 
value e~**’” while the current correspondingly jumps to the value 


C 
ne ef The exponential factor px/v is 


e(E+@)vie=o(Sy9)2) =e 


which will be recognized as the steady state attenuation factor for 
high frequencies. Similarly \/C/Lis the steady state admittance 
of the line for high frequencies. The sudden jumps in the current 
and voltage at time t=2/v are called the heads of the current and 
voltage waves. If, instead of a unit e.m.f., a voltage f(t) is 
impressed on the line at time ¢=0, the corresponding heads of the 
waves are f(O)e~** and +/C/L f(O)e~** for voltage and current 
respectively. These expressions follow at once from the integral 
formula 


= | 1—2)ar ar 
0 
=fO)h() + | “pt—a)a(r)dr, 


The tails of the waves, that is, the parts of the waves subse- 
quent to the time t=/v, are more complicated and will depend on 
the distance x along the line and on the line parameters p and o. 
The two simplest cases are the non-dissipative line, and the 
distortionless line. 
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The ideal non-dissipative line, quite unrealizable in practice, is 
one in which both RF and G are zero. In this case p=o =0, and 
formulas (210) and (211) become 

I=0 for t<2/v, 


nC 
= LOR 070, 


V =0 for t<z/v, 
=1 for t>2z/v. 


Both current and voltage jump, at time t=2/v, to their steady 
values. If an e.m.f. f(¢) is impressed on the line at time ¢=0, the 
corresponding current and voltage waves are 

I=0 fort<2a/v, 


= Sr —z/v) for t>2/2, 


V =0 for t<z/2, 
=f(t—2/v) for t>2/v. 

Consequently the ideal non-dissipative line transmits the waves 
with finite velocity v, without attenuation or distortion. Such a 
line is, of course, the ideal transmission system. 

The non-dissipative line is, of course, purely theoretical and 
unrealizable in practice; the distortionless line is, however, approx- 
imately realizable, and as the name implies, transmits without 
distortion of wave form. The distortionless line is one in which 
the line constants are so related that 


Layee 
ONE POC 
If this condition is satisfied, formulas (210) and (211) become 
l=0fort=<z2/, 
C —ax 7 }o) 
= {2c 1Or t= 0/0: 
V=O for 1 =< 2/0, 


=¢@-° for t>2/v. 
Furthermore, if the impressed e.m.f. is f(t), the corresponding 
current and voltage waves are: 
I=0 for t<z/v, 


= [Sens —£/v) for t>2/v, 


V=0 for t<2/7, 
=¢ 47 f(t—2/v) fori>2/2. 
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The distortionless line, therefore, transmits the waves without dis- 
tortion of wave form, but attenuates the waves by the factor e~°”. 
Such a line is an ideal transmission system as regards preservation 
of wave form, but introduces serious attenuation losses. For 
example, if a line has normally negligible leakage, and leakage is 
introduced to secure the condition R/L=G/C, the line is thereby 
rendered distortionless but the attenuation is doubled. 

One of Heaviside’s most important contributions to wire 
transmission theory was to point out the properties of the dis- 
tortionless line, its approximately realizable character, and to 
base on it a correct theory of telephonic transmission. 

The character of the wave propagation when the parameters 
p and o are not restricted to special values, can only be roughly 
inferred from inspection of the formulas, and then only when 
the properties of the Bessel function J) and J; have been studied. 
Fortunately these functions have been computed and tabulated 
for small values of the argument, and have simple asymptotic 
expansions for large values. It is therefore a simple matter to 
compute and graph a representative set of curves which show 
the current and voltage waves for various values of p, o and z. 
For this purpose it is convenient to introduce a change 
of variables and write: 


T=0t 
G=p/v 
b=c/¥ 


whence the formulas for current and voltage become: 


ie [See 1oVF=a) 
(210-a) 
aC {| ‘et (dV 2dr, 


V= ena + hy CRONE ESE 


a a/72— 42 

Figs. (9) to (18) give a representative set of curves illus- 

trating the form of the propagated current and voltage waves 

for different lengths of line, and different values of the line 
parameters a and b, or p and o. 

The curves of Figs. (9) and (10) show the current entering the 

line in response to a unit e.m.f. applied at time ¢=0. The line 


(211-a) 
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is assumed to be non-leaky (b=0) and is computed for two 
different values of the parameter a. We see that the current 
instantly jumps to the value ~/C/L and then begins to die 


0 =) Values jof wb=t 


0 20 40 60 80 100 120 140 {60 {80 200 


Fie. 9.—Current entering line; EAC =H) (AMEE =O) 


Multiply ordinates by ~/C/L. 


away, the rate at which it dies away depending on and increasing 


_R Ae 
with the parameter a ONT, 


0 Values of vt= 
0 20 40 60 60 {00 120 {40 160 {80 200 


Fig. 10.—Current entering line; # ¢ =a =(0.2645; G—0. 


Multiply ordinates by VC/L. 


If we now consider a point « out on the line, the current is 
zero until r=, at which time it jumps to the value ~W/C/L e~™. 
It then begins to die away provided 2 and a are such that ax<2. 
If, however, we are considering a point at which ar>2, the 
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current begins to rise instead of fall after the initial jump, and 
may attain a maximum value very large compared with the head 


Values of vb=t 
200 220 240 260 280 300 


320 
Fig. 11.—Propagated current in line; « =200; BAe = (—(.OloonG —OE 


Multiply ordinates by a/C/L.e.72:84 


0 Values of vb=t 
180 200 220 240 260 280 300 320 340 360 380 


Fig. 12.—Propagated current in line; x =200; BRIE =0.02645; G =0. 
Multiply ordinates by /C/L.e75-29, 
before it starts to die away. This is shown in the curves of 


Figs. (11), (12) and (13), also computed for the non-leaky line 
(b=0). From these curves we see that, as the length of the line 
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and the parameter a increase, the relative magnitude of the tail, 
as compared with the head of the wave, increases. Finally 
when the line becomes very long, the head of the wave becomes 


ta4t/xPRC 
20 


Fig. 13.—Propagated current in line; aN co=—103G —0. 
Multiply ordinates by 2/Rz. 
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t=4t/x? RC 
1.6 20 


Fia. 14.—Propagated current in cable. 
Multiply ordinates by 2/Rz. 


negligibly small, and the wave, except in the neighborhood of 
its head, becomes very close to that of the corresponding non- 
inductive cable. This is shown in curves (13) and (14), for the 
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line and the corresponding cable, which are plotted to the same 


time scale and ordinate scale to facilitate comparison. 


(15) shows the effect of leakage in eliminating the tail. 


line is not quite distortionless but nearly so. 


180 200 220 240 260 280 300 320 340 360 
Fie. 15.—Propagated current in line; z =200. 


Ee CanG lee 
a= BAC +E, IE 0.08883 


= B/C_ o/h. : 
b= FN Es Vg 70001765 


Multiply ordinates by /C/L.e77-8, 


0 2 4 6 8 {0 {2 14 16 18 


Fira. 16.—Propagated voltage in line; aN z=axr=3;G=0. 


Curve 
This 


An interesting feature of both current and voltage waves is 
that when a sufficient time has elapsed after the arrival of 
the head of the wave, the waves become closer and closer to the 
waves of the corresponding non-inductive cable; that is, to the 
cable having the same Rk, CandG. Consequently the inductance 
plays no part in the subsidence of the waves to their final values, 
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Curves (16), (17) and (18) illustrate the voltage wave for 
several conditions. After the arrival of the head, the wave 
slowly builds up to its final value. Curve (18) represents the 


ERR Sees 
RES eee ass 
Jae aeeese 
SeveGeasoen. 
27 es 
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Fig. 17.—Propagated voltage in line; a w=azr =6; G=0. 


case where the line is very nearly distortionless, showing how 
completely the distorting tail of the wave is eliminated. 

So far we have confined attention to the current and voltage 
waves in response to a unit e.m.f. applied at time ¢=0 to the 


line terminals. Of much greater technical importance is the 
question of the waves in response to a sinusoidal e.m.f. suddenly 
applied to the line terminals. In order to investigate this 
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important problem it is convenient to divide the expressions for 
the current and voltage waves as given by equations (210-a) 
and (211-a) into two components.s We write for r>z2, 


[Ne e-9?-- J (i), (210-b) 

V=e"**+ W(t), (211-b) 
where, by definition, J(¢) and W(t) are the differences between 
the total waves and their heads. The advantage of analyzing 
the waves into these components is that the distortion of the 
waves is due to J(¢) and W(t) respectively, while the first compo- 
nents of (210-b) and (211-b) introduce merely a delay. Thus, 
if the e.m.f. impressed at time é=0 is f(¢), the corresponding 
waves for />2/v or r>2, are 


t= 4[e-eft-2/2) 


fe ON 212) 
a/v 
V=e*f(t—2/0) 
+" sa—ww'andty, (213) 
x/v 


ane) =53( and Wi) =f Wit). 


The integrals of (212) and (213) can be computed and analyzed 
in precisely the same way as discussed in connection with the 
non-inductive cable problem, and are of very much the same 
character as the corresponding waves of the cable. In the 
total waves, however, as given by (212) and (218), a very essen- 
tial difference is introduced by the presence of the first terms, 
which represent undistorted waves propagated with velocity v. 
Thus, if the impressed e.m.f. is sin wt, (212) and (213) become 


Or ees aes 
T=a/Se sin w(t—2/v) 


i 
+{ sin w(t—t,) J’ (1) dt, for t> x/v (214) 
a/v 
V =e-** sin w(t—2/v) 
be 
+] sin w(t—t) W(t) dt, for t> x/v, (215) 
v/v 


Now the first terms of (214) and (215) are simply the usual 
steady state expressions for the current and voltage waves when 
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the frequency is sufficiently high to make the steady state atten- 
uation constant equal to a and the phase velocity equal to v. 
Furthermore the integral terms become smaller and smaller 
as the applied frequency w/ Qa is increased. It follows, there- 
fore, that for high frequencies the waves assume substantially 
their final steady value at time t=2/v, and that the tails of the 
waves, or the transient distortion, becomes negligible. This is 
a consequence entirely of the presence of inductance in the line, 
and shows its extreme importance in the propagation of alter- 
nating waves and the reduction of transient distortion. 

It should be pointed out, however, that if the line is very long 
and the attenuation is very high, the integral terms of (214) and 
(215) are not negligible unless the applied frequency is corre- 
spondingly very high. For example, on a long submarine cable, 
the a.c. attenuation is so large that the first terms of (214) and 
(215) are very small, and J(é) is very large compared with 
1/C/L e~**. Consequently here there is very serious transient 
distortion and alternating currents are therefore not adapted 
for submarine telegraph signalling. 

This discussion may possibly be made a little clearer, without 
detailed analysis, if we recall the discussion of alternating current 
propagation in the non-inductive cable of the preceding chapter. 
From that analysis it follows that, when the applied frequency 
w/2r is sufficiently high, the integral term of (214) becomes 
approximately 


F'(0 
(69) 
and the complete current wave is 
Gee sin N= ead (216) 
L w 
and similarly the voltage wave is 
EArt 2 W(t). (217) 


Now if the total attenuation az is large the last terms of (216) 
and (217), before they ultimately die away, may become very 
large compared with the first terms, which represent the ultimate 
steady state. 

APPENDIX TO CHAPTER VII 


Derivation of Formula (211).—The only troublesome question 
involved in deriving (211) from (207) and (209) is that we have to 
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differentiate with respect to x, in accordance with (207), the 
discontinuous function F(t). To accomplish this we write (209) 
in the form 


F(t) =$(t—2/v)e* Ino Poa?) (209-a) 


where ¢(t) is defined as a function which is zero for ¢<a2/v and 
unity for?>a/v. Clearly this is equivalent to (209) and permits 
us to deal with F(t) as a continuous function. Now, in accord- 
ance with (207), perform the operation of differentiation upon 
(209-a): we get 

OF _ 


0 pee ee 
05 ay $(t—a/v) eI (o-V/ 2? — 2?/v?) 
—09(t— 2/0) 56 1o(oV P= 24/02). 
The first expression follows from the fact that 
0 ted 
ag Pt — 2/2) ec a? (t — 2/2). 
0 
We observe also that a? E—2/e) =0 except at t=2x/v, when it is 
infinite. We also observe that, for ¢>2/v, 
ta 
wet 2/v)at=1 


and that the whole contribution to the integral occurs at t=2/v. 
With these points clearly in mind, the expression 


t 
V=-0 {| Le 
0 Ox 


reduces to (211) without difficulty. 


CHAPTER VIII 


PROPAGATION OF CURRENT AND VOLTAGE IN 
ARTIFICIAL LINES AND WAVE FILTERS 


The artificial line here considered is a periodic structure, 
composed of a series of sections connected in tandem, each section 
consisting of a lumped impedance 2z, in series with the line, and 
a lumped impedance z, in shunt across the line. In the artificial 
line which we shall consider it will be assumed that the voltage 
is applied at the middle of the initial or zeroth section, as shown 
in Fig. 19. This termination is chosen because of its practical 
importance, and because also of the fact that the mathematical 
analysis is simplified thereby. Furthermore any other termi- 


22, 


Fie. 19. 


nation can be regarded and dealt with as an additional terminal 
impedance, so there is no essential loss of generality involved. 

A study of the properties of the artificial line is of practical 
importance for several reasons: 

1. The artificial line is often used as a model of an actual 
transmission line and it is therefore of importance to determine 
theoretically the degree of correspondence between the two. 

2. The solution for the corresponding transmission line with 
continuously distributed constants is derivable from the solution 
for the artificial line by keeping the total inductance, resistance, 
capacity and leakage constant or finite, and letting the number 
of sections approach infinity. 

3. The artificial line is very closely related, in its properties 
and performance, to the periodically loaded line, and its solution 
is, to a first approximation, a working solution for the loaded 


line. 
Aly/ 
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4, The structure is of great importance in its own right, and 
when the impedance elements are properly chosen, constitutes 
a ‘wave filter.” 

We shall now derive the operational and symbolic equations 
which formulate the propagation phenomena in the artificial 
line. Let J, denote the mesh current in the n‘ section of the 
line; Z,_; the mesh current in the (n—1) section etc. Now 
write down the expression for the voltage drop in the n** section; 
in accordance with Kirchhoft’s law we get: 


(2; +220)In a mre ice) =) (218) 


where, of course, the impedances have the usual significance. 

Now this is a difference equation, as distinguished from a 
differential equation, but the method of solution is essentially 
the same. We assume a solution of the form 


I,=Aeé™ +Be™ (219) 


where A, B and I are independent of n, and substitute in (218). 
After some simple rearrangements we get 


{ (21 +222) —2z2. cosh T}-{ Ae” +Be"™} =0 (220) 


Equation (218) is clearly satisfied by the assumed form of solu- 
tion, and furthermore leaves the constants A and B arbitrary 
and at our disposal to satisfy any boundary conditions, provided 
I’ is so chosen that 


cosh T= ee 
220 


=1+2p 


(221) 


where p = 21/422. 
Now by reference to equation (219) it is easily seen that I is the 
propagation constant of the artificial line, precisely analogous to 
the propagation constant y of the smooth line. In terms of the 
impedances z; and z2, the propagation constant of the artificial 
line is determined by (221). This equation may either be 
regarded as an operational equation or a symbolic equation, 
depending on whether the impedances are expressed in terms of 
the operator p or in terms of 7w, where w is 27 times the frequency. 
Now suppose in (221) we write e"=2; the equation becomes 


a+1/z=2(1+2p) 
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and solving for x we get 
x=e' =(1+2p)++/(1+2p)?—1 


=(Vitp+vo)?=(Vitp—Vp) 


which is an explicit formula for I. 

Now return to equation (219) and let us assume that the line 
is either infinitely long, or, what amounts to the same thing, 
that it is closed by an impedance which suppresses the reflected 
wave. We assume also that a voltage Vo is impressed at mid- 
series position of the zero” section (n=0). Equation (219) 
becomes 


(222) 


In=Ae™ 
and the currents in the zero’” and 1** sections are 
Io=A, I;=Aer. 


Now, by direct application of Kirchhoff’s law to the zero‘ 


section, we have 
Vo= (S21 +22)L0 —2e!i, 


whence 
A{kz,+2.(1—e")} =Vo. (223) 
But 
il 
ljg=Al= KY 
is Vo —nT 
= K° . 


where K is the characteristic impedance of the artificial line (at 
mid-series position). Hence by (223) and (222) 
hie = all 1 


Kaz ql =O 
cot iat es ee Le 
222~/ptp? Vee V1+p 
By aid of the preceding the direct current wave can be written as 
— Vo EAT eV ie (225) 
V 21 Cp) V1 +>\p 


This formula is not so physically suggestive as its equivalent 


(224) 


iL 
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but is useful when we come to the solution of the operational 
equation. 

Before proceeding with the operational equation, and the 
investigation of transient phenomena in artificial lines, it will be 
of interest to deduce from the foregoing the unique and remark- 
able properties of wave filters in the steady state. For this 
purpose we return to equation (221) 


cosh [ =1+2p. 


Now suppose that the series impedance 2; is an inductance L and 
the shunt impedance 22 a capacity C, so that, symbolically, 
: w? LC 
2:=twL, of IY gi (aay \ 
and 
cosh T=1—$w?LC. (226) 


Now let us write [=76, where 1=~+/ —1; the preceding equation 
becomes 

cos 6=1—4w?LC (227) 
and the ratio of currents in adjacent sections ise~’. Consequently 
if 62s a real quantity the ratio of the absolute values of the currents 
in adjacent sections 1s unity, and the current 1s propagated without 
attenuation. 

Inspection of equation (227) shows that 6 is real provided the 
right hand side lies between +1 and —1; or that w lies between 
0 and 2/+1/LC. Consequently this type of artificial line trans- 
mits, in the steady state, sinusoidal currents of all frequencies 
from zero to 1/r+/LC without attenuation. It is known as the 
low pass filter. 

If we invert the structure, that is, make the series impedance 2; 
a capacity C and the shunt impedance 2 an inductance L, so 
that 


Pee ae nis pete 
06 2 ean ae, 
we get, corresponding to (226) and (227), 
1 
1 
cos 6=1-a oro’ (228-a) 


This type of artificial line transmits without attenuation currents 
of all frequencies for which the right hand side of (228-a) lies 
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between +1 and —1; that is, all frequencies from infinity to a 
lower limiting frequency 1/4r\/LC, while it attenuates all 
frequencies below this range. It is known, on this account, as the 
high pass filter. 

It is possible by using more complicated impedances to design 
filters which transmit a series of bands of frequencies. We 
cannot, however, go into the complicated theory of wave-filters 
here, which has been covered in a series of important papers.! 
One point should be noted, however: transmission without 
attenuation implies that the impedance elements are non- 
dissipative. Actually, of course, all the elements introduce 
some loss, so that in practice the filter attenuates all frequencies. 
Careful design, however, keeps the attenuation very low in the 
transmission bands. 

We shall now derive the indicial admittance formulas for 
some representative types of artificial lines and wave-filters 
from the operational formula 


ii a . 
Ape eA) —2n, 229 
Jae ee +p+vol (229) 
This equation follows directly from (225) on putting V)=1. 
We start with the so-called low pass filter on account of its 
simplicity and also its great importance in technical applications. 
This type of filter consists of series inductance L and shunt 
capacity C. The general case which includes series resistance 
R and shunt leakage G has been worked out (see Transient 
Oscillations, Trans. A.I.E.E., 1919). The solution is, however, 
extremely complicated and will not be dealt with here. We shall, 
instead, consider the important and illuminating case where the 
series and shunt losses are so related as to make the circuit 
quasi-distortionless. We therefore take, operationally, 


zi=pL+R=L(p+d) (230) 
1/2.=pC+G=C(p+nh) 


where \=R/L=G/C. 
We then have 
mee=L/C, 
2/22=LC(p+A)?, (231) 


L 
p= Coty)? 


1 See Bibliography. 
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Now by reference to formula (229) we see that A, is a function of 
(p+); thus . 


ek et See 5 es 
An= 7 pn ais ceawacrsy 


Now write 
Ik 
Av= . 
" Znl(p) 
It follows at once from reference to theorem V that 
t 
A= (1 +n au) et (232) 
0 


so that the problem is reduced to the solution of the operational 
equation for A2. Writing w.=2/+/LC, we have 


C il ee —2n 
HOS hE ek ae VIF (@p/acl+p/e| (233) 
L V1+(p/w-)? 
_ [Ce [ee tlie 
ND Vp? +o We 
Now refer to formula (mn) of the table of integrals; writing 
./L/C =k, we see by Theorem II that 


il wt 
Aga) |, Tas (234) 
0 


where Jon(r) is the Bessel function of order 2n and argument 7. 
We note also that this is the indicial admittance of the non- 
dissipative low pass wave-filter; that is, the current in the n‘” 
section in response to a unit e.m.f. applied to the initial section 
(n=0). From (232) and (234) it follows at once that 


wt 
Ay=et | J on(7) dr 
‘Jo 


k 
NN t 67 
+R dre | J on(71)dr1. 
0 


Integrating the second member by parts and noting that 
A? (0) =0, this reduces to 


il wet ese 
An= i ace hated (235) 


which is the indicial admittance formula for the quasi-distortion- 
less low pass filter, or artificial line. 
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Before discussing these formulas, it is of interest to derive the 
formula for A? by power series expansion. Formula (233) can be 
written 


ate weit 1 y 
“ Kp} V1+(o./p)?1+-VI1+ (@./p)?)” 
This can be expanded in a series in inverse powers of p; thus 
ee 1 (eye (yr 
i cs Daite\en 


(2n+3) (2n+4) /wo.\ 25 
x ane ) sed. 


Replacing 1/p” by t”/n! in accordance with the Heaviside Rule we 


get 
Ag=7i el Web ow Qn+2 /at\2r+3 
" ki (2n+1)'\ 2 1!(2n+3)!\ 2 


Qn+-3) (Qn-+4) /eet\245 
+ See 1a = ee) 


This can be recognized as the power series expansion of (234). 

The artificial cable is also of interest and practical importance. 
In this structure the series impedance is a resistance FR and the 
shunt impedance is a capacity C, so that 


ai=h, 1/ze=oC, 
212,.=B/pC, 21/22=pRC, (236) 
p=pRC/4. 


Now let us return to formula (229), and expand in inverse powers 
of p: we get 


it 1 2n-+-2 1 
A ee ale 
PA Zep 
(2n+3)(2n+4) 1 
Cnt ayGet) tan | G37) 


: R 
Now since ~/ p21 22 = 4’ We have 


poe BACON i Ee 
*2R\\RCp} 2.1! \RCp 


(Qn-E3)(nta) (2.2 
ery (a) 
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Replacing 1/p” by t*/n! we get finally 


Yen eet Ce Cierfet 
yD Fa Wie) 2.1!(n+1)!\RC 


wee 2t i eee |. (238) 


2.2\(n+2)! \RC 


For large values of n and ¢ this series is difficult to compute or 
interpret. It can, however, be recognized as the series expansion 


of the function 
2t 
eee #l.( ec) (239) 


where [,(2¢/RC) is the Bessel function J of order n and argument 
(2t/RC). This solution, it may be remarked, can be derived 
directly by a modification of the integral formula (7). 

It is beyond the scope of this paper to consider other types of 
artificial lines and wave-filters;for a fairly extensive discussion the 
reader is referred to “Transient Osillations in Electric Wave- 
Filters,” B.S. T. J., July, 1923. The low pass wave-filter, how- 
ever, both in its own right and on account of its close relation to 
the periodically loaded line, deserves further discussion. 

For the non-dissipative low pass wave-filter, we have 


Acai {| te caae (234) 
k 0 
while for the quasi-distortionless low pass wave-filter 
A,= : |, al Ton (r)dr (235) 
0 


where 4 =A/w, = R/Lw, = R/2QvL. 

Computation and analysis of these formulas involve an ele- 
mentary knowledge of Bessel functions. The properties neces- 
sary for our purposes are briefly discussed in an appendix to this 
chapter. 

The indicial admittances for the non-dissipative low pass filter, 
that is, the current in response to a steady unit e.m.f. applied at 
time t=0, are shown in the curves of Figs. 20, 21 and 22, for the 
initial or zeroth, the 3rd and the 5th sections, respectively. These 
curves together with the exact and approximate formulas given 
above are sufficient to give a reasonably comprehensive idea of 
the general character of these oscillations and their dependence on 
the number of sections and the constants of the filter. 
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, It will be observed that the current is small until a time approx- 
imately equal to 2n/w,=n+/L,C, has elapsed after the voltage is 
applied. Consequently the low pass filter behaves as though cur- 
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20 


Fie. 20.—Low pass wave filter. Indicial admittance of initial section (n =0). 
Multiply ordinates by V/C/L. 
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Fig. 21.—Low pass wave filter. Indicial admittance of third section (n =2). 
Multiply ordinates by /C/L. 
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14 16 18 20 22 2 26 28 30 32 34 36 38 40 
Values of wb 
Fig. 22.—Low pass wave filter. Indicial admittance of fifth section (n =4). 
Multiply ordinates by /C/L. 


t5 


10 


rents were transmitted with a finite velocity of propagation 
w,/2=1/+/LC sections per second. This velocity is, however, 
only apparent or virtual since in every section the currents are 
actually finite for all values of time >0, 
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After time t=n+/LC has elapsed the current oscillates about 
the value 1/k with increasing frequency and diminishing ampli- 
tude. The amplitude of these oscillations is approximately 


1/k y 5) 
AL (2n/w.t)? TwWet 
and their instantaneous frequency (measured by intervals between 
zeros) 


av = (2n/wt) 2 


The oscillations are therefore ultimately of cut-off or critical 
frequency w,/27 in all sections, but this frequency is approached 
more and more slowly as the number of filter sections is increased. 


10 200~—«220 240.2~«C«60~SCSC*« 300 +320 +340 360 380 
Fia. 23.—Low pass wave filter. Indicial admittance of 100th section (mn =99). 
Multiply ordinates by V/C/L. 


Figs. 23, 24, 25, give the indicial admittance in the 100th, 500th 
and 1000th section of the low-pass wave-filter. The filter itself 
seldom embodies more than 5 sections. The case of a large 
number of sections is of interest, however, because it represents 
a first approximation to the periodically loaded line. While the 
non-dissipative line is ideal and unrealizable, its study is of practi- 
cal importance because in this type of line the effect of the dis- 
continuous character of the loading of the periodically loaded 
line is isolated and exhibited in the clearest possible manner. 

The dotted curves represent the current in the corresponding 
smooth line. For the smooth line, the current, as we have seen, is 
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discontinuous, being identically zero for a time vt=n and having 
an instantaneous jump to its final value »/C/L at vt=n. The 
current in the artificial or periodically loaded line differs from that 


FoHNAREEEE 
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Values ofw,t 


0 me 
960 980 {000 {020 {040 {060 {080 1100 1120 1140 1160 
Fig. 24.—Low pass wave filter. Indicial admittance of 500th section (n =499) 


Multiply ordinates by V/C/L. 


0 rate 
{960 1980 2000 2020 2040 2060 2080 2100 2120 2140 2160 
Fig. 25.—Low pass wave filter. Indicial admittance of 1000th section (n =999). 


Multiply ordinates by /C/L. 


in the corresponding smooth line in three important respects: (1) 
the absence of the abrupt discontinuous wave front, (2) the 
presence of superposed oscillations, and (3) the absence of a true 


128 ELECTRIC CIRCUIT THEORY 


finite veloctiy of propagation. It will be observed, however, 
that the current in any section is negligibly small or even sensibly 
zero until vt=n, so that the current is propagated with a virtual 
velocity 1/+/LC per section. The presence of a well marked wave 
front is also evident although this is not abrupt, as in the smooth 
line. The effective slope of the wave front becomes smaller 
as the current wave travels out on the line, decreasing noticeably 
as the number of sections is increased. When the number of 
sections becomes large, however, the decrease in the slope is not 
rapid, being in the 500th section about 60 per cent. of that in the 
100th section. 

The superposed oscillations are of interest. These are initially 
of a frequency depending upon and decreasing with the number of 
sections, n, but in all sections ultimately attaining the frequency 


1 y 


rJVLC 7 


which is the critical or cut-off frequency of the line, above which 
steady state currents are attenuated during transmission and 
below which they are unattenuated. When vt is large compared 
with n the amplitude of these oscillations becomes +/1/zvt so 
that they ultimately die away and the current approaches the 
value 1/C/L for all sections. The current in the loaded line is 
thus asymptotic to the current in the corresponding smooth line 
and oscillates about it with diminishing amplitude and increasing 
frequency. 

Since the abscissas of these curves represent values of 2vt= 
2t/+/LC, and the ordinates are to be multiplied by »/C/L to trans- 
late into actual values, the curves are of universal application for 
all values of the constants L and C. 

The investigation of the building up of alternating currents in 
wave filters and loaded lines is very important. It depends for 
the non-dissipative case on the properties of the definite integrals 


wt 
i sin wr Ja(r)dr, 
0 
rat 
| cos wr J_(r)dr. 
0 


where w=w/w, and w=2mr times the applied frequency. The 
mathematical discussion is, however, quite complicated and will 
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not be entered into here. The reader, who wishes to follow this 
further, is referred to Transient Oscillations, Trans. A. I. E. E., 
1919 and Transient Oscillations in Electric Wave Filters, 
Brow eee July. 1923, 


APPENDIX TO CHAPTER VIII 


Note on Bessel Functions.—The Bessel Functions of the first 
kind, Jn(x) and I,(x), are defined, when n is zero or a positive 
integer, by the ge convergent series 


Jn(x) = os 7a | engin Tent NOAH 
~2.4.6(2n+2) CATT “ | j 

L,(2)=5— Beant 2 Tee 
1246. ase) : ; 


In the following discussion of the properties of these functions it 

will be assumed that the argument z is a pure real quantity. 
For large values of the argument (x large compared with n), the 

behavior of the functions is shown by the asymptotic expansions: 


An?—1 , (4n?—1)(4n?—9) 


Be 
I= | 1 Tees + 2(8a)2 
_ (4n?—1) (4n?—9) (4n?—25) , | 
3! (8x)8 P 
J,(«) = Zi? cos (p70) Q,, sin (2— er r 
where 
LA mee = ane 0) 
ine 2! (8x)? 
(4n?—1)(4n?—9)(4n?—25)(4n?—49) 
a 4! (8x)4 
_4n®—1_ (4n?—1)(4n?—9)(4n®—25) 
aaa oe 31(8x)° 


We thus see that J, increases indefinitely and behaves ulti- 


mately as 
e* 


V/ 20m 


130 ELECTRIC CIRCUIT THEORY 


The function J,(x), however, is oscillatory and ultimately 
behaves as 


For all orders of n 


[sated a1) 
0 


The properties of Jn(x) may be described qualitatively as 
follows :— 

When the argument is less than the order (0S 2<n) the func- 
tion is very small and positive, and is initially zero (except when 
n=0). In the neighborhood of «=n, the function begins to 
build up and reaches a maximum a little beyond the point «=n. 
Thereafter the function oscillates with increasing frequency and 
diminishing amplitude, and ultimately behaves as 


When n=0, the initial value is unity, but the subsequent behavior 
of the function is as described above. 

A more precise description of the function is gotten from the 
following approximate formulas. ! 

Jal@) = Bale) cod On(e). «for 21 

where . 
if 1 
=e wi Sane 1/4’ 
TX (1 m 43 m 1 ) 


~ g2 12 a (1 —m?/x?)? 


_ mm. m m? 1 2n+1 
Q,(x) = = =i (le = eres eee 
(2) “ \! oe a cae (2) 4x4 (1 —m?/x?) | aoe 


0,,(2) = £4(2), 


B, (x) = 


m 3m 1 
= V 12 8 of mia? 


1 These formulas are due to Mr. Zobel (See ‘‘Transient Oscillations in 
Electric Wave-filters”” Bell System Technical Journal, July, 1923.) Some- 
what similar formulas had previously been developed by Graf and Gubler 
(Kinleitung in die Theorie der Besselschen Funktionen), and by Nicholson 
(Phil. Mag., 1910, p. 249). 
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and 
m=n?—1/4. 
This approximate formula is valid only where «>n, its accuracy 
increasing with x and with n. For all orders of n it is quite 
accurate beyond the first zero of the function. 
The ‘instantaneous frequency’”’ of oscillation is approximately 


Lie, 1 ie eee 1 

Iq Pa(2) a ae oe? > ax (1 —m?/x2)? 
By this it is meant that at any point 2(#> 7) the interval between 
successive zeros is approximately 7/Q’(x). Otherwise stated, in 
the neighborhood of any point 2, the function behaves like a 
sinusoid of amplitude B,(x) and frequency w/27 where w= ’,(). 

The following approximate formulas, while not sufficiently 

precise for the purpose of accurate computation except for quite 
large values of 2, clearly exhibit the character of the functions 
for values of the argument x>~n, and of the order n> 2. 


J,(2) SAE cos, 2—0,,), 
ee Q 
Jn(@) =—Qn aA sin (¢n2—9n), 


; ne PR 
f. J,(a)dx=1 oye SIN (dnX— On), 


where 


= oO n sin—! (n/2). 


CHAPTER IX 
THE FINITE LINE WITH TERMINAL IMPEDANCES 


So far in our discussions of wave propagation in lines and wave- 
filters, we have confined attention to the case where the impressed 
voltage is applied directly to the infinitely long line. We have 
found that, by virtue of this restriction, the indicial admittance 
functions of the important types of transmission systems are 
rather easily derived and expressible in terms of well known 
functions, and the essential phenomena cf wave propagation 
clearly exhibited. In practice, however, we are concerned with 
lines of finite length with the voltage impressed on the line 
through a terminal impedance Z, and the distant end closed 
by a second terminal impedance Z,. Wenowtakeup the problem 
' presented by such a system. 

Let K=K(p) denote the characteristic operational impedance 
of the line, and y=y(p) the operational propagation constant of 
the line. We have then 


V =Ae-*-+ Be, 


1 il 
I =—Ae-”* ——Be"* 

re Kio? 
where A and B are so far arbitrary constants. To determine 
these constants we assume an e.m.f. # impressed on the line at 
x=0 through a terminal impedance Z, and the line closed at x=s 
by a second terminal impedance Z. At «=s we have therefore 


(240) 


Zol=V 
whence from (240) 
Fe er A — 2 ev B = Ae-v*+ Bev 
and 
= — LS p2 o— 273 ¢ 
[ees risen: A, (241) 

where p2=Z2/K. 

At x=0 we have 

V=E-Zil 


132 


ry 
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whence 
ae ee 
A+B=E Kit RK? 
(rep A Usp b =z, (242) 
where pi=Z,/K. 
From (241) and (242) we get 


wt 1+p2 E 
(1-1) (192) (191) pa e 2” 
oe (So p2)e~*7 _E 
~ (I+en( + p2) — (1—p1)(1—p2)e 2" 


and finally 
e7ve —v(2s—=x) 
eo Hi Thee dp 
K+Z,,_1- pil P2 tne 
l+p1 1+ 2 


If we replace # by a unit e.m.f. we get the operational formula 
for the indicial admittance A,; thus 


aN e-¥*-+ yee y(2s—=z) all 
K 1-mme?* — Z,(p) 


i= (243) 


A,= (244) 
where 
A=K/(K+2Z,), 
es eS Sh 
Mm“ Tte K+Z, 
wl pe ae, 
a 1+): eo, 
K,y,Z1,Z 2,41 and p are, of course, functions of the operator p. 
The integral equation corresponding to the operational formula 


(244) is 
SAGE ii e-?t A, (t)dt. (245) 
« 0 


Now by (244) we can expand 1/Z,(p); it is 


1 e7ve ey (2s—=) 
Zion tne eK 
e~y(2s+2) Pert, 
+Xp1 Me a oe? Te - ov: 
e7v(4s+2) ( ) 


ST Hs erat 
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Now we observe that e~’*/K is simply the operational formula 
for the indicial admittance at point «x of an infinitely long line 
with unit e.m.f. impressed directly on the line at c=0. This will 
be denoted by az(t). Similarly e~’°*-* /K is the operational for- 
mula for the indicial admittance at point (2s—2) with unit e.m-f. 
impressed directly on the line at e=0. This will be denoted by 
Qos—z(t), ete. 

Recognition of this fact allows us to derive a formal solution in 
terms of a series of reflected waves. For let a set of functions 


V0,V1,V2,03, . . . Satisfy and be defined by the operational 
equations 
Vo—= Ap) =r 
V1=Ape 
247 
V2 = Ai be ( ) 


2 
U3 =i, etc. 


It then follows from the preceding and theorem II that 


Vo(t—7)az(7) +01(t—7)a2s_2(7) } 


(248) 
+v2(t—7)dos42(7) +: - 


t 
AS) =5 dr | 
0 


If, therefore, we know the indicial admittance of the infinitely 
long line with unit e.m.f. directly applied and if we can solve the 
operational equations (247), then A,(¢) is given by (248) by inte- 
gration. This solution may well present formidable difficulty in 
the way of computation. It is, however, formally straightforward 
and the numerical computation is entirely possible, the only ques- 
tion being as to whether the importance of the problem justifies 
the necessary expenditure of time and effort. Without any 
computations, however, the solution (248) admits of considerable 
instructive interpretation by inspection. The first term repre- 
sents the current at point x of an infinitely long line in response 
to a unit e.m.f. impressed at a=0 through en impedance Z;; 
Vo =vo(t) is the corresponding voltage across the line terminals 
proper. The second term is a reflected wave from the other 
terminal due to the terminal irregularity which exists there. The 
third term is a reflected wave from the sending end terminal, ete. 
The solution is therefore a wave solution and is expanded in a 
form which corresponds exactly with the sequence of the 
phenomena, which it represents. 
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The solution takes a particularly instructive form when 
Zy=khK and Z,=k,K where k; and ky are numerics. Then 


(249) 


and 
1—k 

a2(t) ote ik, gs—z2(t) 

A ,(t) =—— . (250) 

so 2) (aes eds ag ad 
Te ei eae 

If k,=0, ky =1 we have the case of the e.m.f. impressed directly on 
the sending end of the line and the distant end closed through its 


characteristic impedance; the solution reduces to 
A,(t) =az(t) 


as, of course, it should by definition. 
If k;=0 and k2= ™, we have the case of the line open-circuited 
at the distant end, and the solution reduces to 


A,(é) = {az(t) —2.—2(t) — Aos42(t) + A4e—2(E) ain eyes } ‘ (251) 
Finally, if both k; and kz are zero, the line is shorted and 
A(t) = { de(t)+ass—2(t) +264 a(t) +d4e-2(t) +--+}. (252) 


The operational equations (247) admit of further interesting 
and instructive physical interpretation without computation. 
Consider a circuit consisting of an impedance Z; in series with an 
impedance K. Let a unit e.m.f. be applied to this circuit and let 
vo be the resultant voltage across the impedance K. Then, 
operationally, 

VE 


mz 


V0 
so that vy, thus defined in physical terms, is the vp of equations 
(247). 

Now let this voltage be impressed on a circuit consisting of 
an impedance 2Z, in series with an impedance K —Z;, so that the 
total impedance is K+Z.. Let the resultant voltage drop across 
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the impedance element K —Z, be denoted by 11; then operation- 


ally . 
igs ae 
OS Kk eae 


which agrees with v1 as given by equation (247). 

Similarly if voltage v; is applied to a circuit consisting of an 
impedance 2Z, in series with an impedance K — Z; and if v2 denote 
the voltage drop across impedance K —Z;, then 


Vo=Apipbe 


We can thus see physically what the voltages v9,v1,v2, . 
mean in terms of simple circuits consisting of K and Z,; in 
series and K and Z, in series respectively. 

I shall now work out a specific problem exemplifying the 
preceding theory. The example is made as simple as possible for 
two reasons. First because its simplicity makes it more instruc- 
tive than when the phenomena depicted and the essentials of 
the mathematical methods are obscured by complicated formulas 
and extensive computations. Secondly while the general method 
of solution illustrated is thoroughly practical we cannot hope to 
arrive at the numerical solutions of the complicated problems 
without a large amount of laborious computations. Problems 
involving transmission lines with complicated terminal impe- 
dances are among the most difficult, as regards actual numerical 
solution, of any which present themselves in mathematical 
physics. On the other hand, the formal solution (248) gives at a 
glance the essential character of the phenomena involved. 

The specific problem we shall deal with may be stated as follows: 
A unit e.m.f. is directly impressed on the terminals of a transmis- 
sion line of length s, the distant end of which is closed by a con- 
denser Cy. The line is supposed to be non-dissipative, its 
constants being inductance L and capacity C per unit length. 
Required the current at any point a(7<s) of the line. 

Wewrite WL/C =k, 1/+/LC =v: then by virtue of the preceding 
analysis of transmission line propagation the indicial admittance 
a, of the infinitely long line is given by 


a,=0, for t<a/v, 
L 
k 


POL bw vi 
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The operational characteristic impedance is, of course, k=+/L/C, 
and the terminal impedances Z; and Z, are given by 


Zi = 0, 
Zo = 1/pCo. 
Referring now to equation (244) we have: 
ve il. Ly iF 


_k-1/pCo_kCop-1 
ie EV eh eRe 


Consequently, referring to equations (247), we have, operationally, 


Vo=1 


peo eetine | 


2 kCop+1 
ee, =a aa)) 
: : kCop+1 
ee ea 
a Ca 


In order to determine these functions we have therefore to solve 
the general operational equation 


ay CT, 
tae ey) 


where JV,, denotes either von—1 OF Ven. 
In order to eliminate the coefficient kCy we make use of theorem 


VI, and write 


In accordance with that theorem 
Valt) = bn(t/kCo). 
We therefore start with the operational equation 
Woes 1 n 
(at) 
Now the solution of this operational equation is very easy and can 
be expressed in a number of ways. We require it expressed in the 
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form most easily computed. The following appears best adapted 
for our purposes. Consider the auxiliary operational equation: 


n n)(n—1 
== (o—25 prt-g2 2) “ c 1) pias 


1 
Some — | nn 23 
qpoeea(=1) Me 


The explicit solution is gotten by placing 1/p” by t”/n! and p” by 
d" /dt”, whence 


d” d” 
on(t) = & Pare ee Thee 12") 


tite ee ea 


(Die Ne 
ag ( p ) ~ A(p) 


iia. pi ‘ a 

So \ net) eae 1) 
Dal a 

pS pa Lippe 1) 

me 1 Died want 

= (te ee aay 


Referring now to theorem V, we see that 


t 
dn(t) =(1+ [ at)ont - eT 


Since we have already solved for on(t), this determines ¢,(t) 


But writing 


it follows that 


and hence V,(t). The functions v,v1,2 . . . are therefore 
determined. 

Now refer back to equation (248) giving the required current in 
terms of vo,2 . . . and the admittances az(t),dos—2(tn), . 


It follows at once by substitution of the preceding that 


A.(t)=% v(t 8) ¢0i( 2) (1242) 4. o. | 


the functions 9,v:,v2 being zero for negative values of the argu- 
ment. This result may possibly require a little explanation. 
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Consider the expression 


7 ices SORE 


where 1(¢) denotes a function which is zero for ¢<¢) and unity for 
t>t. It is evidently identical with the admittance a,(t) pro- 
vided the proper value is assigned to fp. 

Now since 1(t) =0 for t<f and unity for t> ty the preceding may 
be written as zero for t<¢) and 


call f(t—7r)dr for t >to 


which is equal to f(t—to). 
If we set x=0, we get the current entering the line; thus 


srt 2) 
mete 
= 142vi(t— 7) 427; (-?) 
+20.) +- is ; 


This has been computed for the case where 1/sLC,=10 and is 
shown in Fig. 26. Referring to this figure we see that the current 


10 


08 


06 


sl = Total Inductance of Line 
sC= » Capacity 

C= TerminalCapacity 
1VsLCo = 0.4 
Multiply Ordinabes by VC/L | 


a4 


02 


Values of t/nVSLC, 
S 0 Of 02 03 04 05 06 07 08 09 


10 
Fig. 26.—Current entering non-dissipative line terminated by capacity Co unit 
e.m.f. applied to line. 


jumps at ¢=0 to the value »/C/L=1/k, and keeps this constant 
value for a time interval 2s/v. At this instant the first reflected 
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wave arrives and the current takes another jump of 2/k. There- 
after it begins to decrease very slowly until time t=4s/v at which 
time it takes another jump of 2/k. Thereafter we have a series 
of jumps of 2/k at time intervals 2s/v, the current decreasing 
between successive jumps. The smooth curve is the indicial 
admittance of an oscillation circuit consisting of an inductance sL 
in series with a capacity Co. We see, therefore, that the current 
in the line oscillates with discontinuous jumps about the current 
in the corresponding oscillation circuit. Since the whole circuit 
contains no resistance, the oscillations never die away, but con- 
tinue to oscillate, as shown, about the curve 


Nia (Jez) 


which is the indicial admittance of the corresponding oscillation 
circuit. 

I shall now discuss a method of solving circuit theory problems, 
quite generally applicable to complicated networks, and particu- 


Network 


Fiq. 27. 


larly useful in dealing with transmission lines terminated in 
impedances. I have found it particularly useful in arriving at 
numerical solutions where other methods prove far more laborious. 
It is also of mathematical interest, as it applies another type of 
integral equation to the problems of electric circuit theory. 

Suppose that we have a network with two sets of terminals as 
shown in Fig. 27.1. Now suppose that terminals 22 are short 
circuited and a unit e.m.f. inserted between terminals 11. Let 
the resultant current flowing between terminals 11 be denoted by 
Su(t)=Si. and that between terminals 22 by So:(t) =Sa. 
Siu is the driving point indicial admittance with respect to 
terminals 11 and S.; the transfer indicial admittance of terminals 
22 with respect to 11 under short circuit conditions. 

Similarly if terminals 11 are shortcircuited and a unit e.m-f. 
inserted between terminals 22 the current flowing between termi- 

‘Regarding conventions as to signs, see the Appendix to this chapter. 
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nals 22 is denoted by S22(t)=S22 and that flowing between ter- 
minals 11 by Si2(¢)=Si2. If the network is passive, i.e., contains 
no internal source of energy, it follows from the reciprocal theorem 
that So.=Si. As far as the two sets of terminals are concerned, 
the network is completely specified by the indicial admittances 
S11,S22,Se1 = Sie. 

Now let a voltage V(t) = V1 be inserted between terminals 11, 
and a voltage V.2(t) =Ve2 between terminals 22. The current 
flowing between terminals 11, denoted by J, is 


at a 
I=" i} Vilr)Sult—r)de +5 |, ACR Cie On 
0 0 
while the corresponding current between terminals 22 is 
Orie a [* 
I) =4, {| Va(r)Sox(t—1)dr+ 4, il Vi(r)Sa(t—7)dr. (254) 
0 0 


Now consider two networks of indicial admittances S11,S22,Si2 
= Sp», and 111,722,712 = T21 respectively and let them be connected 


Fie. 28. 


in tandem as shown in Fig. 28 to form a compound network. 
We require the indicial admittances of the compound network in 
terms of the indicial admittances of the component networks. 

Short circuit terminals 22 of the compound network and insert 
a unit e.m.f. between terminals 11. Let V3i(¢) denote the result- 
ant voltage between terminals 33 measured in the direction of the 
arrow, and J; the current flowing between the networks. We 
have then the two following expressions for the current J;. 


t 
L=Sn—%y | V31(7)Se0(t —7)dr (255) 
0 
and 
d t 
heal V31(7) Ti1(t—7)dr. (256) 
0 
Equating, we get 


iu, Vaxlr)[Saa¢—7) + Tuut—7)|dr=Sar(t). (257) 
‘J0 
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By precisely similar reasoning, if terminals 11 are short circuited 
and a unit e.m.f. inserted between terminals 22, and the corre- 
sponding voltage across terminals 33 denoted by V32, we have! 


dt _ 
£ il Va Sn@ pe tera (258) 


Equations (257) and (258) are integral equations of the Poisson 
type which completely determine V3: and V3.2 in terms of the 
indicial admittances S and 7. We shall discuss the solution of 
these equations presently. 

If U11,U22,U21=U12 denote the indicial admittances of the 
compound network we have at once 


it : 
Uu= Su | V31(7r)Si2(¢—7) dr (259) 
‘ 0 
1 t 
Up=Tat) a Veo(r) Taslt—1) dr (260) 
0 


and 


di fe 
Uae = U2 eal V 31(7) T2i1(t—r)dr 
0 (261) 
= + | Vse(r)Si2(¢—7)dr. 


If, therefore, equations (257) and (258) are solved for V3; and 
Vso, the required indicial admittances of the compound network 
are given by (259), (260) and (261) in terms of the indicial admit- 
tances of the component networks. 

A simple example will now be worked out illustrating the 
method of solution just discussed. Suppose that a unit e.mf. is 
impressed on a transmission line (infinitely long) of distributed 
constants #,L,C, through a terminal resistance Ro. Required 
the terminal line voltage V. 

The operational equation of this problem is gotten in the usual 
manner. The current entering the line is 


Cp 


YNIp +R 


' Vso being opposite to Vs; in direction. 
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It is also obviously equal to p > (1- V): equating the two expres- 


sions, and rearranging we ae 


Writing R/2L =p and setting R)=+/L/C, this becomes 


V/1+2p/p 
14 1+29/p 


This operational equation can, of course, be solved in a number 
of ways, though, as a matter of fact, its numerical solution is 
quite troublesome. This point will be returned to later; we shall 
first formulate the problem in accordance with the method just 
discussed. 

The indicial admittance of the line is known; it is 


af et Talot =Ai(é). 


(262) 


Consequently the current entering the line is explicitly 


d 
£ alk Va gs 


But the current is also equal to p * (= V(t)); equating, we get 
V()=1-Ros mm ave one 
t Jo 
Performing the indicated differentiations 
t 
V(t) =1-—R, A(0)V@) -Rs{ V(r) A’(t—r)dr. 
0 
Now A(0) = Re and 
, 7 C 
A’ (t) = pe [Ii (pt) Sie 


and Ro=+~/L/C; therefore the equation becomes 
t 
viy=5+8 i} V(t—1)o(or) —i(or)]e~*r dr. 
0 
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As a matter of convenience we change the time scale to pt, and get 
t 

V@)=5+5 {| V(t=1)llols) Liye“ ar (263) 
0 


ke _—- 
sara drV (t— 1) £07" Lo(r) 
where it is understood that ¢ is actually pt. This is the integral 
equation of the problem and is in the canonical form of Poisson’s 
integral equation. 

Before solving this equation numerically I shall show how a 
simple approximate solution is obtainable immediately; an advan- 
tage often attaching to this type of integral equation. 

The function a e~'I,(t) is equal to —1 for t=0 and converges 
rapidly to zero. V(t) has, as we know from the operational equa- 
tion, the initial value 1/2 and the final value 1. Neither function 
changes sign. It follows from the mean value theorem that the 
equation can be written as 


a a) 
iG =>—— i EU 
V(t) , a! o{ die I, (t)dt 
where aSl. Integrating 


V(@)=5-5V Wlem*'Zo(ot) 1) 


and 
iL 


mrs ae 


(264) 
The correct initial and final values of V(¢) result for all final 
values of a1; so that approximately 


ee ene 
VOTH TO 
This equation, while not exact, except for t=0 and ¢ very large, 


shows faithfully the general character of V(¢) and the way it 
approaches its final value unity. For large values of t 


e-*Io(t) =1/+/2nt 


whence 
] 


clara yavse tS aie 
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Approximations of the foregoing type are not always possible 
and may not be of sufficient accuracy. I shall therefore give 
next a method of numerical solution which is generally applicable 
to integral equations of this type and works quite well in practice. 
We shall write the integral equation in the more general form 


AG Sion { TCO (265) 


where f(x) and k(y) are known and u(x) unknown. The method 
depends on the numerical integration of the definite integral. 
Let us divide the z scale into small intervals d and for convenience 
write 


u(nd) = Un 
S(nd) =fn 
k(nd) =ka. 


Now from the integral equation we have at once 
u(0) = Uo =fo 
d 
u(d) =u =it| u(d—y)k(y)dy. 
0 


Now if d is taken sufficiently small 


d 
[e@—wpkaday =$ tas herb u0ke 
0 


whence 
m=fits [wi kot+uoksl 
and 
Vet ae [fituokid/2) 
which determines uw; since wu is known. Similarly 
U2=fotd[d uo ko+turkit3 uz ko] 

which determines wz. Proceeding in the same manner 

Us =f3tdl¥ uo ks tus ket U2 ki+3us ky], ete. 


In this way we determine the value of u(x), point by point from 
the recurrence formula 
u _fntadlzuokn Ht hint Ue Kin—9F > FUn—1 ka), (266) 
“ 1—$kod 
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The result of the application of numerical integration, in 
accordance with formula (266), to the integral equation (263) is 
shown in Fig. (29). The dotted curve is a plot of the approxi- 


Values of pt 
0 {0 20 30 40 5.0 6.0 70 80 90 100 


Fie. 29.—Line terminal voltage unit e.m.f. impressed on line through resistance. 


Ro = V/L/G. 


mate solution as given by equation (264), for a=1. We see 
that the voltage starts with the value 1/2 and slowly reaches its 
ultimate value, unity, its approach to unity, for large values of f, 
being in accordance with the formula 


ai 
my 1+1/./2nt 

The application of the foregoing method to the transmission 
line problem proceeds as follows. Let Si:(¢), Soo(é) and Sy.() 
be the short indicial admittances of the line. S1:(¢) is the current 
entering the line (at «=0) with unit e.m.f. directly impressed 
and the distant end short circuited. Sj2(t) is the current at 
x=s under the same circumstances. Consequently from (252) 


S(t) = ao(t) + 2ao5(t) +2a4,(t) +: - - 
Sie(t) = 2{ a,(t) + as.(é) +as5s(t) -+- Snel } : 
(267) 
Soo 1s clearly equal to Si; by symmetry. 

Now suppose that an e.m.f. #=E(t) is impressed on the line 
at x=0, t=0, through a terminal impedance Z;, and the distant 
end (x=s) closed through an impedance Zp». We suppose these 
terminal impedances and the actual impressed e.m.f. replaced 
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by the actual line voltages V; and V2, impressed directly on the 
line at ~=0 and at x=s; the terminal currents are 


Lol) fg | Sit—2)V lo) 
0 


d t 
a i RG Vade 


(268) 
Tyee S i SCENE 
0 
ae (269) 
“ai fe Eels 


But the current at x=s is also equal to the current in the terminal 
impedance Z, in response to the terminal voltage V2: denoting 
by a(t) the indicial admittance of Z, it is 


L(t)= S ii Vad (270) 
0 


Similarly the current entering the line at =O is the current 
flowing in the terminal impedance Z, in response to the e.m-f. 
E—V'. Denoting by a:(t) the indicial admittance of Z:, it is 


In) =5, | i eViGae = ameeTD 
0 


Equating equations (268) and (271) and (269) and (270) we 
eliminate J,(t) and J,(¢) and get 


[tsud—a) exten IViear— fSilt—2) Valea 
0 0 
= | ault—2) Beer, (272) 
0 


= [ Sult—nViG)drt | [Sat—2) —aalt—))Valdr=0. (273) 
0 0 


These two equations are simultaneous integral equations of the 
Poisson type in Vi and V2, which completely determine these voltages 
provided the admittances and the impressed voltages are known. 
They therefore represent the application of a new type of integral 
equation to the problem of electric circuit theory. 

The numerical solution of the general case, either by (248) or 
(272-273) is necessarily laborious when the terminal impedances 
are complicated and is only justified when the technical impor- 
tance of the problem is considerable. I wish, however, to empha- 
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size two points in this connection: the numerical solution is 
always entirely possible and, compared with other and older 
forms of solution, enormously simpler. One has only to inspect 
the classical forms of solution of problems of this type to realize 
the truth of this last statement. 

I shall now give two applications of equations (272-273) to 
specific problems, in one of which an approximate solution of the 
integral equation can be gotten, and in the other of which numer- 
ical integration is applied. 

Problem I. Given a non-inductive cable of distributed con- 
stants C and FR and length s, with unit e.m.f. applied at x=0, 
while at «=s the cable is closed by a condenser Co. Required 
the terminal voltage V(t) across the condenser Co. 

We first write down the short-circuit indicial admittances of 
the cable; from equation (168) of a preceding chapter and equation 
(267) they are:— 


Sit) = Soo(t) 


C 246) sues 868 
malay Tt2e7 +27 $20 + |, TA) 
Sro(t) = Sei(t) 
en Vos) a ey eee 212) ) 
m2) 5 [evete tte et} (275) 


where B=s?RC/4. 
Now the current at «=s is equal to 


t 
Sie(t) Sema, cdl V(r) Soa(t se: t)dr. 


It is also the condenser current due to the voltage V(é); that is 
d 
Co qv 
Equating the two expressions and integrating we get 
t t 
CoV (t) =| Su(r)ar— f V (7) Soo(t—1) dr (276) 
0 0 


which is the integral equation of the problem. In order to get 
an approximate solution without detailed computation we assume 
that the cable is long. In this case the leading terms of (274) 
and (275) are large compared with the terms following: Further- 
more Si2(t) builds up very slowly while S20(¢) is a rapidly varying 
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function. A good approximation therefore results if we take 
V(r) outside the integral sign in (276) and write 


CoV (t) = [ sul dr — V(b J Sule)ar 


| Sloat 
Vee (277) 
Mi Toe cae Sea(t)at 


This approximation is quite good for long cables and shows the 
way V(t) builds up quite truthfully. We see that V is initially 
zero, and builds up ultimately to unity. For large values of ¢, it 
becomes 


whence 


Voj=—_. (278) 


This is the approximate formula also for the open circuit voltage, 
as may be seen by setting C)=0 in (277). 

In electric circuit problems, it is often sufficient, as implied 
above, to know qualitatively the behavior of an electric system 
without going through the labor of detailed computation. For 
this purpose the formulation of the problem as a Poisson Integral 
Equation is particularly well adapted. A simple example 
will be given, which can be checked from the known solution. 
Suppose that we required the voltage V at point x of an infinitely 
long transmission line (L,R,C) in response to a unit e.mf. 
impressed at x=0. This is, of course, known from formula 
(211-a): we shall here be concerned, however, with the approxi- 
mate solution from the Poisson integral equation of the problem. 

If a,(t) denote the indicial admittance of the line at point z, 
then the current at point x is simply a(t), which is given by form- 
ula (210-a). But if V(t) is the voltage at point z, the current is 
also given by 


d T)AQ\’U—T AT. 
| vera )d 


Equating these two expansions, we get the integral equation of 
the problem 


‘ f Vea 7dr =a), 
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RI¢ 
Nowif we write 7’ = pt—A where p=R/2LandA=",/°, then 


a.=4/Se-TH® Ion TOF IA), T>0, 


and in terms of the relative time 7’, the integral equation is 
reducible to 
a 
dT 
while the exact formula for V is by (211-a) 


[vane tiarm ent TV TO FRA) 
0 


Pies oy ae 
Ue ee sa | CECE, 
0 WV 7(7 +24) 
From the integral equation it is easy to establish superior and 
inferior limits for V(7); it is 


elicy? Te ey 
Vise. * TT) 


[le rncmv nar 


= rT 
ib e- Ti yar 
0 


Both formulas give the correct initial and final values of V; 
namely e~“ and unity. Since V lies between V. and Vs, the 
mean value (Va+V,)/2 also has correct initial and final values 
and should be a better approximation than either. The table 
given below shows the orders of approximation obtainable from 
the case where A=3. It is evident from this table that the 
foregoing approximate formulas exhibit the form of V(7’) qualita- 
tively in a quite satisfactory manner. 


ez Lae 


=V,(T). 


7 Va Vo 2(VatVo) va 

0 0.05 0.05 0.05 0.05 
2 0.25 0.12 0.18 0.17 
4 0.39 O18 0.29 0.26 
6 0.50 0.23 0.36 0.32 
8 0.57 0.27 0.42 0.37 
10 0.64 0.31 0.47 0.41 
12 0.69 0.34 0.51 0.44 
15 0.74 0.38 0.56 0.48 
18 0.78 0.41 0.60 0.52 
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Problem II. Our second illustrative problem may be stated 
as follows:—A unit e.m.f. is impressed on a transmission line of 
length s and distributed constants L,R,C. At x=s the line is 
closed by a resistance Fy in parallel with an inductance Lp. 
Required the current in the terminal resistance. If V(t) denotes 
the terminal voltage, the current at «=s is given by 


Sio(t) -4 V (r)Se2(t—7)dr. 


It is also equal to the current flowing into the terminal impe- 
dance; that is 


VO+p | vend 
Ry ise F Te Ane 


Equating and rearranging 


| A +50) |PO=Su— | VC] P4807) far (279) 


Now the short circuit admittances Sez and Sizare given by formula 
(210-a) of a preceding chapter, and S.(0)=+/C/L. In order 
to apply numerical integration to (279), numerical values must 
be assigned to the constants. We take 


Ro=~V/L/C =1935 ohms, 
I, =0.4 henry, 


R 

oy, P= 292, 
Het A LO = 1105-105, 
s== 0). 


The results of the numerical evaluation of equation (279), with 
these values inserted, is shown in Fig. 30. The voltage is identi- 
cally zero until v6 =100; t=100/v is the time of propagation of 
the line. At that instant it jumps to the value e-"=e71%/» 
and then begins to die away rapidly due to the draining action 
of the inductance. The effect of secondary reflection is insignifi- 
cant and therefore not shown. The current in the terminal 
resistance is V/Ry so that it is given by the same curve. 

I have reserved until the last the exposition of the expansion 
theorem solution as applied to transmission lines with terminal 
impedances, for the reason that it is the least powerful and the 
most restricted, although most closely resembling the classical 
form of solution. Furthermore, it does not represent the 
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sequence of physical phenomena, in fact it is not a wave solution, 
but a solution in terms of normal or characteristic vibrations. In 
practical application its usefulness is restricted to the non- 
inductive cable. 
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Fig. 30.—Voltage across terminal impedance on smooth line. 


It will be recalled that the expansion theorem solution is 
formulated as follows :— 
If 
A=1/Z(p) 


is the operational equation of the problem, then the explicit 
solution is 


n 


ere 6 
ZO) de Z' (pe) 
a 
where pi,p2 . . . are the roots of the equation Z(p) =0. 


Let us apply this formula to the case of a line of length s, with 
unit e.m.f. directly applied at «=0, and line short circuited at 
x=s. Referring to equation (244) and putting \=4,;=p2.=1 
we get 


ee 1 cosh 7 (Sau) eal 


K sinhys = Z,,(p) 
as the operational formula of the problem. This can be written 
as 


(280) 


cosh y(s—«)_ 1 
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where in the general case, 
y= (Lp+R)(Cp+6@). (282) 


The values of y for which Z,(p) vanishes are the roots of the 
transcendental equation 


sinh ys=0 


excluding zero. These roots are infinite in number: Let ym be 
the m™ root; then 


Ym=t) m=1,2,..., 0, (283) 
The corresponding values of p» are then gotten by solving (282) 


for p and writing y = ym. 
The explicit solution of the operational equation (281) is then 


cosh Ym(S—%) 
Aa ray tote.) cay 


SYm ” ID. cosh Ym 


(284) 


h Ym 
= 70) + DCH E/ Dm) Fe oom 


™ 


colt dp 


Let us apply this to the non-inductive, non-leaky cable in which 
L=G=0 and y =~/RCp, so that 


Op ep ee 
Baa nl BOS RG 
and 
im ee 
UCT ee 
Also Z,(0)=sR. We thus get 
A j=. S ccs oo #ne' (285) 
Te aR iee aie Sta ; 
m=1 


This is a thoroughly practical formula for computation, owing 
to the rapid convergence of the series. In fact, for this particular 
line termination chosen, it is probably the simplest and most 
easily computed form of solution. These advantages depend, 
however, strictly on two facts. First, the fact that the line is 
taken as non-inductive and secondly, that the terminations 
chosen are those of a short circuit. In fact, as we shall see, it is 
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only in the case of the non-inductive cable that this type of 
solution is of any practical value. 

There is one other point which should be carefully observed 
in connection with this solution (285). This is that it is not 
expressed in terms of a series of direct and reflected waves, 
corresponding to the sequence of physical phenomena, but in 
terms of normal or characteristic vibrations. This point will be 
returned to later. 

Let us now attempt to apply this type of solution to the 
transmission line, L,R,C,G. Writing 


Ip 1G 
Shoe 
Lia! 
oS 2020 
v=1/VLC 
we have 
il 
Be Page) ae! 
whence 


BE 
Pin =P. tole 


_ |/mm? o? 
=—ptw ===) m=12.... 
8 vy? x 
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ns 


Dk cies pene dim 


Setting G=0 for simplicity and substituting in (284) we get, 
after easy simplifications, 


nO cos (& r) aN TS 
A,(t)= sR Rt aS a ae a (4/22) eee coe 


If we write 
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(286) can be written as 


mr s§ 


2th ee et ff | mmr es mr 
OSS ae te (Um vt ) + Sin~—(Hm vt +x) : 
S 


™ 


(287) 


This type of solution is often referred to as a wave solution and the 
component terms of the series regarded as travelling waves. As 
a matter of fact it is a solution in terms of normal or characteristic 
vibrations, each of which is to be regarded as instantaneously 
produced at time t=0. The solution in terms of true waves has 
been fully discussed in the preceding. 

Formula (287) is practically useless for computation on account 
of the slow convergence of the series (the series are only condi- 
tionally convergent), and cannot be interpreted to bring out the 
existence of the actual direct and reflected waves and the physical 
character of the phenomena it formulates. In fact, as stated 
above, this form of solution is useful only in connection with the 
non-inductive cable. 

In the cases considered above we have taken the simplest 
possible terminations—these of short circuits in which case the 
roots of Z(p) are easily evaluated. If, however, the line is 
closed by arbitrary impedances, the case is quite different, and 
the location of the roots becomes, except for simple impedances, 
and then only in the case of non-inductive cable, practically 
impossible. While, therefore, the expansion theorem solution 
can be formally written down, its actual numerical evaluation is 
a practical impossibility, except in a few cases. For this reason 
it will not be considered further here. 

The physically artificial character of the expansion solution, 
as applied to transmission lines, may be seen from the following 
considerations. When a wave is sent into the line, for a finite 
time equal to the time of the propagation of the line, it is inde- 
pendent of the character of the distant termination. Yet in the 
expansion solution every term involves and is dependent upon 
the impedance constituting the distant termination. Evidently, 
from physical considerations, the series of component vibrations 
making up the complete solution must therefore so combine as 
to annihilate the effect of the distant termination for a finite 
time. The solution is, therefore, mathematically correct but 
physically artificial. 
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Note on Integral Equations. 


An integral equation is defined as an equation in which the 
unknown function occurs under a*sign of integration; the process 
of determining the unknown function is called solving the 
equation. 

Integral equations are of great importance in mathematical 
physies and in recent years very considerable work has been done 
on them from the standpoint of pure analysis. 

The types of integral equations with which we are concerned 
in the present work are Laplace’s equation 


F(p) =|} eP f(t)dt 


and Poisson’s equation 


ROLOL i Ok G nar 


But little work has been done on Laplace’s equation from the 
standpoint of pure analysis; its most extensive and useful appli- 
cations appear to be in connection with the operational calculus. 
Practical methods of solution are extensively discussed in the 
text. 

We shall now briefly discuss the solution of Poisson’s equation. 

The formal series solution, which is absolutely convergent, 
is obtained by successive substitution. Thus suppose we write 


(0) = Gol) + 41(2) +onl@)+- 

and define the terms of the series in accordance with the scheme 
bo(x) =f(x), 
62) = [dK @—way, 


al I Moiaeanen aie 


the resulting series satisfies the integral equation and is absolutely 
convergent. It is, however, practically useless for computation 
or interpretation. 
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A power series solution, when it exists, can be gotten by 
repeated differentiation; thus 


$(0) =f(0), 
$' (x) =f" (x) + (0) K(x) + {| ¢'(x—y) K(y)dy, 
$'(0) =f’(0) ++ 6(0) K(0) 


In this way all the derivatives at r=0 are calculable; let them 
be denoted by ¢o, $1, ¢2 . . . Then 


2 
d= dotdiz torts: 


This form of solution, also, is of limited practical usefulness, 
except for small values of z. 

A number of mathematicians, including Whittaker and 
Bateman, have studied the question of numerical solution and 
suggested other processes. After quite extensive study of the 
question, however, the writer is of the opinion that point-by- 
point numerical integration like that discussed in the text is, in 
general, the most practical, rapid and accurate method of numer- 
ical solution. This judgment is confirmed by G. Prasad who, in 
a paper on the Numerical Solution of Integral Equations 
before the International Mathematical Congress (Toronto, 
1924), discusses the whole question and arrives at the same 
conclusion. 

In the text, numerical integration is carried out in accordance 
with Simpson’s Rule. It is possible, of course, to employ more 
complicated and refined formulas for approximate quadrature. 
It is the writer’s opinion that this is hardly justified in practical 
problems and that the required accuracy is more simply obtained 
by employing smaller intervals. 


APPENDIX TO CHAPTER IX 
Note on Conventions as to Signs in Networks. 


In the network shown on page 140 the arrows indicate the 
directions chosen as positive in the network itself, quite regardless 
of the presence of any e.m.fs. and currents. 

The sign attributed to a current, an e.mf., or a voltage is 
positive if the current, e.m.f., or voltage is in the positive direc- 
tion; otherwise the sign is negative. 
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Stated more fully: 

A current at a specific point (at a specific instant of time) is 
positive if it is flowing in the positive direction; negative if 
flowing in the negative direction. 

An e.m.f. or a voltage between two points is positive if the 
potential increases in the positive direction between the two 
points; negative if the potential increases in the negative direc- 
tion. (It may be noted that this convention makes the sign of a 
voltage the same as the sign of that e.m.f. which could be inserted 
between the two points without producing any effects in the 
network.) 


CHAPTER X 


INTRODUCTION TO THE THEORY OF VARIABLE 
CIRCUITS! 


In the preceding chapters it has everywhere been assumed 
that the networks are invariable: that is to say, that the constants 
and connections of the network do not vary or change with 
time. In many important technical problems, however, we 
wish to know, not merely what happens when an electromotive 
force is applied to an invariable network, but the effect of sud- 
denly changing a circuit constant or of introducing a variable 
circuit element. In the present chapter we shall show that this 
type of problem can be dealt with by a simple extension of the 
methods discussed in the preceding chapters. 

The simplest and at the same time one of the most technically 
important problems of this type is the effect of sudden short 
circuits and sudden open circuits on an energized network or 
system. This type of problem will serve as an introduction to 
the more general theory. 

The Sudden Short Circuit.—Consider the network shown in 
Fig.31. This network, which for generality is supposed to consist 


Fie. 31. 


of two parts A and B, indicated schematically, is energized by an 
electromotive force E(t) which produces a voltage V(t) between 
the points ab. The voltage V(é) is calculable by usual methods 
from H(t) and the constants and connections of the network, 
supposed to be specified. 

1 The material in this chapter is largely taken from a paper by the writer on 
“Theory and Calculation of Variable Electrical Systems,” Phys. Rev. 


Feb. 1921. 
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We now suppose that, at reference time t=0, a short circuit is 
suddenly placed across ab; and require the effect of this short 
circuit on the distribution of currents in the network. The 
solution of this problem is based on the following proposition: 

The effect of the short circuit is precisely the same as the insertion 
at time t=0 of a voltage— V(t), equal and opposite to V(t), between 
points a and b. 

The resultant currents in the system for t2 0 are then composed 
of two components :— 

(1) The currents which would exist in the invariable network, 
in the absence of the short circuit, due to the impressed source 
E(t). These are calculable by usual methods. 

(2) The currents due to the electromotive force V(#) inserted 
at time t=0, between the points a and b. These are also calcu- 
lable by usual methods, since V(¢) is itself known from the 
primary distribution of currents and charges. 

By the preceding analysis we have succeeded, therefore, in 
reducing the problem of a sudden short circuit, to the determina- 
tion of the currents in an invariable network in response to a 
suddenly impressed electromotive force: that is, the problem 
to which the preceding chapters have been devoted. 

The Sudden Open Circuit.—The problem of a sudden open 
circuit in any part of a network can be dealt with ina precisely 


analogous manner, although the actual calculation of the result- 
ant current and voltage distribution is mathematically more 
complicated. Consider the network shown in Fig. 32. Here 
the network is supposed to be energized by an electromotive 
force E(t) which produces a current I(t) in the invariable network 
in branch ab. We require the effect of suddenly opening this 
branch. The solution of this problem depends on the following 
proposition. 

The effect of opening branch ab at reference time t=0 is the same as 
suddenly inserting at time t=0, a voltage V(t) which produces in 
branch ab a current—I(t) equal and opposite to the current which 
would exist in the branch in the absence of the open circuit, 
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While this proposition is precisely analogous to the correspond- 
ing proposition in the case of a sudden short circuit, it does not 
explicitly determine the voltage V(t), which must be calculated 
as follows: 

Let the driving point indicial admittance of the network, as 
seen from branch ab be denoted by Aa(t). Then, from the 
preceding proposition, it follows at once that V(t) is given by 


d t 
oF {| V7)A a(t —7)dr= —I(t),t >0. 


This is a Poisson integral equation in V(t), from which V(é) is 
calculable. With V(t) determined, the currents in any part of 
the network are calculable by usual methods, and consist of two 
components :— 

(1) The current distribution in the network due to the 
impressed source H(t) in the absence of the open circuit. 

(2) The current distribution due to the electromotive force 
V(t) inserted in branch ab at time t=0. 

As in the case of the sudden short circuit, we have thus reduced 
the problem of a sudden open circuit to the determination of the 
current distribution in an invariable network in response to a 
suddenly impressed electromotive force. 

Variable Circuit Elements.—-In the preceding cases of sudden 
open and short circuits it will be observed that the network 
changes discontinuously from one invariable state to another. 
A more general case, and one which includes the preceding as 
limiting cases, is presented by a network which includes a 
variable circuit element: that is, a circuit element which varies, 
continuously or discontinuously, with time. A network which 
includes such a variable circuit element will be called a variable 
network. Variable circuit elements of practical importance are 
the microphone transmitter, which consists of a variable resist- 
ance, varied by some source of energy outside the system; the 
condenser transmitter, which consists of a condenser of variable 
capacity; and the induction generator, in which the mutual 
inductance between primary and secondary, or stator and rotor, 
is varied by the motion of the latter. The case of a variable 
resistance will serve as an introduction to the general theory 
of such variable networks. 

Consider a network, energized by a source H(t) in branch 1, 
and containing a variable resistance element r(¢) in branch n. 
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The functional notation r(t) indicates that the resistance 7 varies 
with time. Let J, (t) denote the current in branch n, and assume 
that the network is in equilibrium prior to the reference time 
t=0. The mathematical theory of this network depends on the 
following proposition :— 

The network described above can be treated as an invariable network 
by eliminating the variable resistance element r(t) and inserting an 
electromotive force —r(t)I,(t): that zs, an electromotive force equal 
and opposite to the potential drop across the variable resistance 
element. Consequently the current in the variable resistance 
branch is determined analytically by the integral equation 


meals Bry eel ee de 
L@=5 i} E(7)Aw(t—1)dr—— i r(r)In(t)Ann(t—7)dr. (288) 


The first component is simply the current [,(t) which would exist 
in the variable branch if the variable element were absent; 
hence, dropping the subscript n for convenience, the current in 
the variable branch is given by the integral equation 


I) =Io(t) - 4 i} “FOIA eae (289) 


and the voltage across the variable element by 
v(t) =r) )- (290) 
Having determined /(t) and v(t) from this integral equation, the 
distribution of currents in the network is calculable as that due to 
a source H(t) in branch 1 and a source v(t) in branch n of the 
invariable network: that is, the network with the variable resist- 
ance element eliminated. 
A very simple example will serve to illustrate the foregoing :— 
Into a circuit of unit resistance, and inductance L=1/a, in 
which a steady current J, is flowing, a resistance r is suddenly 
inserted at time ¢=0: required the resultant current J(t). In 
this case we have: 
A(t) =indicial admittance of unvaried circuit 
= 1-e-*!, 
ri) =r, 


and the integral equation of the problem is: 


*t 
1 = tr | -eyit—w)ay 
t 0 


t 
=Iy—ra | I(¢—y)e—* dy. 
Jo 
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If the solution is carried out as indicated by (291) below, and if 
the notation at=x is introduced, we get without difficulty 
I(t) =Iy{ 1—r(1—e,(x)e-*) +1r?(1 —e9(x)e-*)—r3(1 —e3(x)e-*) + 


where the function e,(x) is defined as: 


€n(v) =1+a/1!+4?/2!4+-43/3!4+---+a7-1/(n—1)! 


=first n terms of the exponential series. 
For all finite values of the resistance increment r the series can be 
summed by aid of the identity 
l—¢,(z)e*= {| dxe—*x"—-1/(n—1)! 
0 


Substitution of this identity gives 


I(t) =1(1 -r{ ctor) 
0 


I 1+re-“Gtn) 
= 1 g— ie —- 


Equation (289) is an integral equation of the Volterra type, 
which includes the Poisson integral equation as a special case. 
Its formal series solution is obtained as follows:—Assume a series 
solution of the form 


I) =1@)-hO+h®)-hO®)+-:-- (291) 
and define the terms of the series by the scheme 
Lo=4, i ads 
t Jo 
Be ee es ie en haa (292) 
d t 
Te O= iy | TOR(YA (oar 
0 


Direct substitution shows that this series satisfies the integral 
equation. Furthermore, it is easily shown that it is absolutely 
convergent. 

While this series solution is not, in general, well adapted for 
numerical calculations, it throws a good deal of valuable light on 
the ultimate character of the oscillations in the important case 
where H(t) and r(t) both vary sinusoidally with time. In this 
case, if the frequency of the applied e.m.f. be denoted by F and 
that of the resistance variation by f, it is easy to show that the 
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current I,(t) in the unvaried circuit is ultimately! a steady state 
current of frequency F. This follows from the fact that the 
definite integral which defines the current Jo(é) is resolvable 
into the ultimate steady state current corresponding to an 
applied force of frequency F, and the accompanying transient 
oscillations which ultimately die away. ‘The fictitious e.m.f. 
which may be regarded as producing the component current 
I,(t) is r(t)Io(t); this is ultimately the product of the two fre- 
quencies F and f, and therefore resolvable into two terms of 
frequency F+f and F—f respectively. Carrying through this 
analysis, it is easy to show that each component current is 
ultimately a steady state but poly-periodie oscillation, as indi- 
cated in the following table: 


Component Current Frequency 

Di ee eh ee F, 

y Ea meta rae eer A F+YS,F—-f, 

Site RAC nee Aa oN F+2f,F,F —2f, 

Lge yen eset re eee: F+3f,F+/,F—f,F—-3f, 

Ye res ae aera een ate eee F+4f,F+2f,F,F—2f,F—4f. 


It is of importance to observe that the component currents 
involve, from a mathematical standpoint, multiple integrals of 
successively higher orders, the nth component /,(¢) involving a 
multiple integral of the nth order with respect to J)(t). Conse- 
quently the successive currents require longer and longer inter- 
vals of time to build up to their proximate steady state values, 
so that the time required for the resultant steady state to be 
reached cannot be inferred from the time constant of the unvaried 
circuit. 

From the preceding table it will be seen that the ultimate 
steady state current is obtained by rearranging the series J>-+J,-+ 
I, and is of the form 

+0 
= An cos (Q+nw)t +B, sin (Q+nw)t 
n=— oO 
where Q=27F and w= 2rf. 


1Tt hardly seems necessary to remark that the reference time {=O is 
purely arbitrary and that the resistance variation may start at such a time 
thereafter that Jo(t) may be regarded as steady state during the entire 
time interval in which we are interested. Going farther, if we confine our 
attention to sufficiently large values of ¢, the whole process may be treated 
as steady state. 
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It is interesting to note that this series comes within the defini- 
tion of a Fourier series only when F =0 or an exact multiple of f. 
The steady state solution is of very considerable importance and 
is considered in more detail in a succeeding section. 

From the foregoing we deduce an outstanding distinction 
between the variable and invariable networks. In the latter the 
currents are ultimately of the same frequency as the impressed 
e.m.f., whereas in the former they are ultimately of an infinite 
series of frequencies. 

In the preceding example, the variable impedance element is a 
resistance r(t). If the variable element is taken as an inductance 
\(t) the voltage, corresponding to equation (290) is 


d 
ql. 


The case of a variable capacity element is handled as follows: 
Let 1/C=S and assume that S is variable; thus, S=So+o(¢). 
The drop across the variable condenser element is then 


v(t) =o(t) | Todt 


Similarly a variable mutual inductance p(t) between branches m 
and n produces the voltages 


d 
FAOLEO 
in branch m, and 


SHOIn(l 


in branch n. This case may be illustrated by: 


The Induction Generator Problem 

In a sufficiently general form, this problem, which includes the 
fundamental theory of the dynamo, may be stated as follows: 

Given an invariable primary and secondary circuit with a 
variable mutual inductance Mf(t) which is an arbitrary but 
specified time function, and let the primary be energized by an 
e.m.f. H(t) impressed on the circuit at the reference time ¢=0: 
required the primary and secondary currents. 

In operational notation the problem may be formulated by the 
equations: 

Zili—pMfOl2= Lb), 
—pMf(t)I1+Z 22 I, = 0, 
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in which Z,; and Zo. are the self impedances of the primary and 
secondary respectively; Mf(t) is the variable mutual inductance; 
E(t) is the applied e.m.f. in the primary; and p denotes the 
differential operator d/dt. By aid of the fundamental formula 
these equations may be written down as the following simultane- 
ous integral equations: 


hO= Hl WAnt—v (BO) +My, IOeW)) 


as (ee d 
I,(t) = Maj ‘ dyA .2(t—y) ay [fwliy)l. 


In these equations, Ai:(¢) and A»2(t) denote the indicial admit- 
tances of the primary and secondary circuits respectively (when 
M=0): that is, the currents in these circuits in response to a 
unit e.m.f. (zero before, unity after time f=0). We assume, of 
course, that they are known or can be determined by usual 
methods. 

It follows at once that the formal solution of these equations is 
the infinite series: 


hd) =Xot)+X2@)+Xst)+---+Xoant)t+-:-:: 
I.) =Yi+Y3;0+Y,0)+--- 


in which the successive terms of the series are defined as follows: 


5 Oates 
Xo =F [dvd wll NBW) =H, 
0 
Vite Me apdnG i eee 
1(t) = tgp ee YW) gy Fu Xow)l, 
iy || arlene ee ee 
A? = dt " YA y dy IY 1 y)I; 


Gale d a 
YQ= MG | dyAnlt—y) 7, UX), ete 


In the light of formula 
d t 
10=7 | SWAC—way 
‘Jo 


the physical interpretation of the series solutions follows at once: 
thus, NXo(t) is equal to the current /(t) flowing in the zsolated 
primary in response to the applied e.m.f. H(t); the first component 
current Y,(¢) in the secondary is equal to the current which 
would flow in the isolated secondary in response to the applied 
e.m.f. M(d/dt)f()Xo(t); Xe(t), the second component current 
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in the primary, is equal to the current in the isolated primary in 
response to the applied emf. M(d/dt)f(t)Yi(t); ete. The 
resultant currents are thus represented as built up by a to-and- 
fro interchange of energy between primary and secondary, or 
by a series of successive reactions. In the important case where 
the applied e.m.f. and the variation of mutual inductance are both 
sinusoidal time functions, of frequency F and f respectively, it is 
easy to show that each component current becomes ultimately 
equal to a set of periodic steady state currents. Thus the com- 
ponent Xp, is ultimately singly periodic, of frequency F; Y, is 
ultimately doubly periodic, of frequencies F+f and F—f; Xe 
triply periodic, of frequencies F+2f, F and F—2f; Y3 quadruply 
periodic, of frequencies F+3f, F+f, F—f, F—3f; ete. 

The Solution for the Steady State Oscillations——For the very 
important case of periodic applied forces and periodic variations 
of circuit elements we are often concerned exclusively with the 
ultimate steady state of the system, and not at all with the 
mode in which the steady state is approached: that is, attention 
is restricted to the periodic oscillations which the system executes 
after transient disturbances have died away. In this case, if 
the periodic variations of circuit elements are sufficiently small, 
the required steady state is obtained in the form of a series by 
replacing each term of the complete series solution by its ultimate 
steady state value; a process which is very simple in view of the 
physical significance of each term of the latter series. The 
appropriate procedure will be briefly illustrated in connection 
with the variable resistance element. In view of the fact that 
we are concerned only with the ultimate steady state oscillations, 
we can base the solutions on the symbolic equation 


fener (293) 
Z 
Here r(t) is the variable resistance element; J is the current which 
would flow in the absence of the resistance variation; and Z is a 
generalized impedance of the network, as seen from the variable 
branch. Its precise significance and functional form is given 
below. 
We now suppose that /, is given by 


ae (real part) (294) 
=1(Jy Cn een) (295) 
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where the bar indicates the conjugate imaginary of the unbarred 
symbol, so that (295) is entirely real. Correspondingly the 
variable resistance will be taken as 


\ 


r(t) =) (gist + @—tet) 


= re (real part) (296) 
=r COS wt. 


Here r is taken as a pure real quantity, which fixes the size of the 
resistance variation. No loss of generality is involved in this, 
since it merely involves referring the time scale to the zero of the 
resistance variation. 

The symbolic impedance Z, as employed in the theory of alter- 
nating currents, will depend on the frequency and is, in general, a 
complex quantity. Its value at frequency 2/27 will be denoted 
by 
while its value at frequency (Q-++nw)/27 will be written as 

Z(i(Q+nw)) =Zn. 
We now assume a series solution of (293) of the form 


ee tee 


where the terms of the series are defined by the symbolic 
equations 


r(t) 
I,= y vr 
pene (297) 
_ _ rt) 
Ini = F LF 


Substitution shows that this series formally satisfies the equation. 
Starting with the first of (287) and substituting (295) and (296) 
we get 


Ne ie (etotte—tat) (Fg et@tt Jy e— tat) 


= 47 FogiOtott Fye Karat Fy glll=alt-+ Jyg-Ha—a)} » (298) 


or 
et (2+w)t ei(Q2Q—w)t } 


oe = — ee 
. Z1 + pipe 


5 Jo (299) 
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In (299) it is to be understood that the real part is alone to be 
retained. 
Proceeding in a similar way with the equation 


we get me ' 
r\?2 ei(Qt2w)t pi(N—Bw)t giMt / | fi 
h=(5) Ze | Lai) BSAA Gtz) 

In this way the steady state series solution is built up term by 
term, the component currents being poly-periodic as indicated 
in a previous table on p. 164. 

For sufficiently small impedance variations this method of solu- 
tion works very well, and leads to a rapidly convergent solution. 
In other cases, however, the solution so obtained may be diverg- 
ent, even when the complete series solution from which it is 
derived is absolutely convergent. The explanation of this lies 
in the fact that the steady state series so obtained is the swm of 
the limits (as t approaches infinity) of the terms of the complete 
series solution, whereas the actual steady state is the limit of the 
sum. These are not in general equal; in particular the former 
may be and often is divergent when the latter is convergent. 

In view of the foregoing considerations it is of importance to 
develop another method of investigating the steady state oscilla- 
tions which avoids the difficulties in the formal series solution. 
The following method has suggested itself to the writer and works 
very well in cases where the previous form of solution fails. It 
should be stated at the outset, however, that the absolute con- 
vergence of the solution to be discussed, while reasonably 
certain in all physically possible systems, has not been established 
by a rigorous mathematical investigation, which appears to 
present very considerable difficulties. 

We start with the problem just discussed and, in view of the 
results of the formal series solution there obtained, assume a 
solution of the form: 


(300) 


N 
Hl z = : 
te 4(2-+-mw)t —i(2+mo)t 
I 9 S Ane +Ame (301) 
—N 


N 
i Dee hee (real part). (302) 
-N 
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Here the series is supposed to extend from m=+WN tom=—N. 
Ultimately, however, N will be put equal to infinity. As before, 
the bar indicates the conjugate imaginary of the unbarred symbol 
and (301) is therefore entirely real. 

If we now substitute (301) in the symbolic equation (293) we 
get, by (295) and (296), 


; = | Ametabme Aa orien | =1J ye M*1L J, e-imt 


= spe he) = | Apne te Ae ae aUe) 


Simplifying this equation and dropping the conjugate imagi- 
naries gives :— 


; . if c 
> Anette Toe SA, eil@+(m+1)e)t (304) 


if p 
“F > A menat (m—1)e)t, 


Z(4(Q+-mw)) = Zn 


Finally, if we write 


and 
} Lan = lin (305) 

and equate terms of the same frequency on the two sides of the 
equation, we get 

PAG =— los A Neat 

Aan —Rm( Ami A min) 0<\m <N (306) 

Agp=Jo—hp(A_,+A)1). 
It will be observed that, by (805), starting with Ay each coeffi- 


cient is determined in terms of the coefficient of the next lower 
index. Thus: 


A v= hy A gal 
Ay-1 = hy-1(Ay-2+Ay) 
— —hy-1Ay—2 
x 1— hy—thy 0G 
Similarly 
_ __hw-2An-s 
Se Be FEA Sc ne 


1—hy_ihy 
Continuing this process it is easy to show that, for positive indices 
(m positive), 

A n= h m Cie A m—1 (308) 
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where C,,, designates the continued fraction 
1 


1-— hee ye 1 


iL 


Piri lran2 


———— Ves 
1—hy—1hy 
The procedure for the coefficient A_,, is precisely similar. For 
convenience we write A_,=A'n, Z-m=Z m, and T/Z_m=h'm. In 
this notation we get by precisely similar procedure 


Al Cin ind (309) 
where C’,, designates the continued fraction 


1 


a 1 
il i ton tre aa fy a fy ce 
e e ails Lh! mph’ myo. 


ese 
1—h'y—th'y 


We now put the index N equal to infinity and the continued 
fractions C,, and C’,, become infinite instead of terminating 
fractions. 

Collecting formulas we now have 


Ap = Nig Cm Amat (310) 
AG = San aA pal 
and 
Ag=Jo— (ho Art+ho'A's) 
whence 
— Jo . 
~ L=hohiCi—ho'hi’ Cr’ 
The coefficients are thus all determined in terms of Jo. 

The practical value of this method of solution will depend, of 
course, on the rate of convergence of the continued fractions. 
While no rigorous proof has been obtained, it is believed that they 
are absolutely convergent for all physically possible systems, but 
this question certainly requires fuller investigation. Neverthe- 
less any doubt regarding the convergence of the solution need 
not prevent the use of the method in a great many problems 
where physical considerations furnish a safe guide. or example 


Ao 
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this method of solution, when applied to the problem of the 
induction generator, discussed above, leads to the usual simplified 
engineering theory of the induction generator and motor, besides 
exhibiting effects which the usual treatment either ere or 
fails to recognize. 

Non-linear Circuits.—In the previous examples discussed, the 
variations of the variable circuit elements are assumed to be 
specified time functions, which is the same thing as postulating 
that these variations are controlled by ignored forces which do 
not explicitly appear in the statement and equations of the problem. 
We distinguish another type of variable circuit element, where 
the variation is not an explicit time function but rather a func- 
tion of the current (and its derivatives) which is flowing through 
the circuit. For example, the inductance of an iron-core coil 
varies with the current strength as a consequence of magnetic 
saturation. The equation of a circuit which contains such a 
variable element (provided it is a single valued function) may be 
written down in operational notation 


ZI+$(1) =E(), 
or 


ZI=E()—4[1 (0). (311) 


In this equation Z is, of course, to be taken as the impedance of 
the invariable part of the circuit, the indicial admittance of which 
is denoted by the usual symbol A (t). 

Equation (311) may be interpreted as the equation of the 
current /(¢) in a circuit of invariable impedance Z when subjected 
to an applied e.m.f. E(t)—¢[I(4)]; consequently, by aid of our 
fundamental formula, /(t) is given by 


aia Cals 
1) =4 | AG—wE@day—4 | 'AG—woll ley 
0 0 
The first integral is simply the current in the invariable circuit of 


impedance Z in response to the applied e.m.f. H(t); denoting this 
by I(t), we have 


1()= 1) —4 [Aol 


This is a functional integral equation, the solution of which is 
gotten by some process of successive approximations. For 
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example, provided the sequence converges, J(¢) is the limit as n 
approaches infinity of the sequence 
T(t), Ii), I,(t),. ery I,(t), 


where the successive terms of the sequence are defined by the 
relations: 


t 
Iy(t) = I(t) — All A(t—y)@[Io(y)|dy. 


Tni(t) = Io(t) -#/ A(t—y) o[In(y)]dy. 


We shall not pursue the discussion of non-linear circuits 
further, in view of their mathematical complexity and their 
relatively specialized technical interest. The reader who is 
interested may, however, consult the writer’s paper on Variable 
Electrical Systems,! for a fuller treatment of the subject. 


1 Phys. Rev. Feb., 1921. 
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CHAPTER XI 


THE APPLICATION OF THE FOURIER INTEGRAL TO 
ELECTRIC CIRCUIT THEORY 


The application of Fourier’s series in electrotechnics is a com- 
monplace; the use of the Fourier integral, however, has largely 
remained in the hands of professional mathematicians. An 
outstanding distinction between the series and the integral, 
from which the greater power of the latter may be inferred, is 
that the series represents only a periodic regularly recurrent 
function, whereas the integral is capable of representing a non- 
periodic function: in fact all types of functions, subject to certain 
mathematical restrictions which are usually satisfied in physical 
problems. 

Before taking up the application of the Fourier Integral to 
Electric Circuit Theory, we shall very briefly review the elemen- 
tary mathematics of the series and integral; for a fuller treatment 
the reader is referred to Byerly, Fourier’s Series and Spherical 
Harmonices.! 

Consider a function ¢(¢), which in the region 0<t< T is finite, 
single-valued and has only a finite number of discontinuities or 
of maxima or minima. In this region it can then be expressed 
as the Fourier series 


p(t) 4+ > A, COs ri)+B, sin( 27") (312) 
1 


ae Q27rn 
Asal o(t)- cos G at 


yd fe . (2Qrn 
B= All o(t)+ sin ( ia tat 


‘In this chapter the Fourier integral is approached from the view-point of 
its physical application and no completeness or rigour is claimed for the 
treatment. The mathematical theory of the Fourier integral is, of course, 
completely developed in treatises on the subject. The object of this chapter 
is merely to outline some of its applications. 

174 


where 


(313) 
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An equivalent series is 


, 2 
g(t) =4F + ye Coe (2r"1-0,) (314) 
1 
where 
F, = VA2+B?, (315) 


Oo; = tan (B,/ A’). 


This expansion is valid in the region 0<t<T7, irrespective of 
the form of the function elsewhere. Let us, however, assume 
that the function repeats itself in the period 7’: that is 


OU) y= ojos b= 1 2. da ais 
Then the expansion represents the function in the region 
—NT<t<(N+ 1)T. 


Finally if N is made infinite, the function is truly periodic and 
the Fourier series represents it for all positive and negative 
values of time. 

It follows from the foregoing that, if the Fourier series repre- 
sents the function for all positive and negative values of time, 
the function must be periodic for all positive and negative values 
of time; otherwise the expansion is valid only over a restricted 
range of time. ; 

Now let us suppose that ¢(¢) isnon-periodic. For convenience, 
in connection with subsequent applications we shall suppose that 
it is zero for all finite negative values of time, that it converges 
to zero as t—o, and that 


{ |o(t) |e 
JO 


exists. Such a function obviously cannot be represented by the 
usual Fourier series for all finite positive and negative values of 
time; it can be represented, however, by the limiting form 
assumed by the series as the fundamental period 7’ is made 
infinite. That is, we can assume that the function is periodic 
in an infinite fundamental period and this will not affect the 
expansion for finite positive and negative values of time. Pro- 
ceeding in this way and putting the fundamental period T 
equal to infinity in the limit, the Fourier series (314) becomes an 
infinite integral and we get 


$(t) =) [Fe) - COS (wt—O(w)) dw (316) 
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where 


eae [eo Cone uf+| [eo ret ae (317) 


tan 6(w) = |, $(t) sin wt us [ ¢(t) cos wt dt. (818) 
0 0 


This is the Fourier integral identity of the function ¢(¢) and is 
valid for all finite positive and negative values of time. 

In physical applications, particularly those to electric circuit 
theory, it is often convenient to employ exponential instead of 
trigonometric functions. The required transformation follows 
easily from the relation 


e’ =cos 0+7 sin 8, t=+~/ —1. 
Thus if we write 27/7’ =w) the Fourier series (312) is easily 
reduced to the form 


and 


+o 
B= SF Ginremr! (319) 
where 
he : 
F(inwo) = All (rem neu dr. (320) 


In precisely similar manner the Fourier integral (316) can be 
written as 


o(t) = i F (iw)-e dw (321) 
25 sorydr f ase da. (322) 


Applications to Electric Circuit Theory.—Let us assume that 
at time t= —NT, an electromotive force E(t), periodic in funda- 
mental period 7’, is impressed on a circuit of complex impedance 
Z(iw), where w denotes 27 times the frequency. Required the 
resultant current J. 

For values of t>—NT the electromotive force (see formula 
(319)) can be expressed as the Fourier series 


E(t) = >) Finare™* 
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where 
ii T 
F (ina ) => all E(r)e~eordr. 
0 


The resultant current for t> —NT is therefore 


transient oscillations 
{= Se inet | initiated at time | 
Z(inwo)- i= —NT. | 


If we are concerned with the current for values of t>0, and if 
NT is made sufficiently large, the initial transients will have 
died away and the complete current for t>0, will be given by 


ebony 
Z (ine) 


—o 


einwot, (323) 


This formula implies the periodic character of E(t) for suffi- 
ciently large negative values of time. If, however, H(t) is zero 
for negative values of time, we can employ the Fourier integrals 
(821) and (322) in precisely the same way and get, as the com- 
plete expression or the current for positive or negative values 
of time :— 


Sie . 
[= J Fite (324) 


|" E(r)dr IL ies (325) 


The infinite integrals (3824) and (325) formulate the current in the 
network, specified by the impedance function Z(iw), in response to 
an electromotive force E(t) impressed at time t=0; they therefore 
mathematically formulate, by aid of the Fourier integral identity, 
the fundamental problem dealt with in the preceding chapters and 
solved by aid of the operational calculus. 

No attempt will be made here to discuss the solution of the 
infinite integral (325), which is usually a problem presenting formi- 
dable difficulties, even to the professional mathematician. The 
general method of solution is by contour integration in the com- 
plex plane and the calculus of residues. By this process it has 
been successfully applied to the solution of special problems, and 
also to deriving some general forms of solution such as the expan- 
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sion theorem solution.!_ Compared, however, with the opera- 
tional calculus, it has no advantages from the standpoint of rig- 
our, and lacks entirely the remarkable simplicity and directness 
of the Heaviside method. 

In the direct solution of circuit problems, therefore, it is 
believed that the application of the Fourier integral is attended 
by few if any advantages, and presents formidable mathematical 
difficulties. On the other hand, there are certain types of prob- 
lems encountered in circuit theory, where the Fourier integral 
is a powerful tool. These will be briefly discussed. 

The Indicial Admittance Expressed as a Fourier Integral.— 
Without recourse to the general formulas (324) and (825), it is 
a very simple matter to formulate the indicial admittance A (¢) 
(and the Heaviside Function h(t)) as a Fourier Integral. To 
this we now proceed, starting with the known identity 


2 (°% sin wt 
(44 |, = da) =1 i>0 (326) 
0 


T 
=0 t>0 


Evidently this expression represents the unit function, zero 
before, unity after, time t=0. 

Now let us suppose that the unit e.m.f., as formulated by 
(326), is impressed on a network of impedance Z(iw), and let us 
write 


1 E29 (w) , att 
yee, ale (327) 
so that a(w) and 8(w) are the two components of the steady state 
admittance at frequency w/2rz. 
il ; 
Now the first term, » of (326) is of zero frequency and pro- 


1 ; : 
duces a current 9 a(0). The integral term contains components 


of all frequencies, that lying between w and w+dw, being 
ore sin wt.dw 
TW 


and the corresponding current is 


2 1 
7 w|Z(iw)| 
1 Bush, “Summary of Wagner’s Proof of Heaviside’s Formula.’ Proc. 


Inst. of Radio Engineers. Oct., 1917. Fry. ‘The Solution of Circuit 
Problems.” (Phys. Rev. Aug., 1919.) 


sin (wt+0(w)dw. 
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Consequently the total current, which is the indicial admittance, 
is given by 


met 1 (“sin (wat +6(w)) 
A()=5a(0)+5- {| Gt (328) 
=5e eel 2 sin wt dw (329) 


II ke) cos wt dw. 
T 0 w 


The same process applies, obviously, to the Heaviside function 
h(t), representing current, charge, potential drop, or any of the 
circuit variables with which we may be concerned. 

Thus if the generalized steady state impedance is represented 
by H(tw), and if we write 


ra a =a(w)+1b(w), 


we have 


h(t) =30(0)+" | ate) sin wt dw (330) 
T JO 
+-[ oo) cos wt dw. 
Tv 0 Ww 


Formula (329) for the indicial admittance is then simply a 
special case of formula (330) for the Heaviside function h(t); 
that is, the response of the network to a unit e.m.f. impressed 
at time t=0. 

By (330) we have, changing the sign of ¢ 


h(—t)= 5000) —_ | sin wt dw (331) 


+ Bie) cos wt dw 
T JO w 


By subtraction and addition of (330) and (331) we get at once 
h(t) —h(—t) -"| a sin wt (332) 
T JO [63) 


nin +h(—=a(0)+- | Bie) cos wt, (333) 
rT Jo > @ 
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Now let us suppose that ¢>0; then since h(t) is identically 
zero for negative values of t, it follows that h(—1) is zero in equa- 
tions (332) and (333), and we have 


h(t) = al ey sin wt (334) 
0 
=a(0)+ =| BAC ke (335) 
CUS A) (6) 


provided t>0. It is therefore only necessary to evaluate one of 
the infinite integrals of equation (330). As an interesting corol- 
lary, the behavior of the network under all circumstances is completely 
determined if either the real or imaginary component of the complex 
steady state admittance is specified over the entire frequeny range. 

The analytical evaluation of h(t) from (334) or (835) is usually 
quite complicated, unless the infinite integral can be recognized, 
and involves advanced mathematical methods. Its numerical 
evaluation can sometimes be effected without prohibitive labor 
by the following process, applied to either (834) or (335). Select- 
ing (334), it can be written as 


a(t) = i “Je(*) eae (336) 


and when ¢ becomes very large, this approaches 


AOE i a Sn ae (337) 
T JO 


which is, of course, the ultimate steady value of h(t). Now 
write (336) in the form 


h(t) =" 2) va (*) sin \.dd ae 


12a 


where x may have any finite real value we please. Now having 
assigned some particular value to x, the first integral of (338) 
can be evaluated mechanically or numerically. The second 
integral can then be expanded asymptotically by repeated 
partial integrations as follows: 


ot ; [o/) 
JokeQ) sinxan=— ff sae ). d cos x 
== ha(4) 00 +f x dh a(*) an 
ro \ 7g) cose COs Ay 


THE APPLICATION OF THE FOURIER INTEGRAL 181 


Repeating this process we get for the second integral term of 
(338) the series 


2, ae. Oe i fay 
2 cosaj1—- ee... |da(*) (339) 


ee pe 
ax ax re t 


which may be used for computation for sufficiently large values 
of ¢ and x. 

This process for evaluating (338) is usually laborious but is 
often quite practicable. It will be observed that the computa- 
tation becomes easier for large values of t. On the other hand, 
the evaluation of A(t) by the operational calculus is usually 
simpler for small values of t. This suggests that in some prob- 
lems the numerical solution may often be advantageously 
effected by a combination of the operational calculus and the 
Fourier integral, using the former to calculate h(t) for small 
values of ¢, and formula (338) to calculate it for large values of 
t. In practical problems, therefore, the Fourier integral may 
be of great value in supplementing the operational calculus. 

The Building-up of Alternating Currents.—The preceding 
analysis can be generalized and extended to the case of sinusoidal 
impressed forces by employing the expressions 


Leilani 
cos of S42 oN (340) 


and 
sin wf 342 {| aoa (341) 
2 Tv 0 Xr 


Formula (340) represents a function which is zero for negative 
values of time ¢t, and equal to cos wt for positive values of time 
t, while (341) represents a function which is zero for ¢<0O and 
sin wt for t>0. They can therefore be applied to formulating and 
investigating the behavior of networks in response to suddenly 
applied sinusoidal electromotive forces. 

The application suggested by the preceding is of considerable 
practical importance because in long transmission systems the 
caleulaton of the building-up of alternating currents 1s extremely 
complicated and laborious. On the other hand an investigation 
of the mode of building-up is necessary in long communication 
systems since on the character and duration of the building-up 
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process, the excellence of a signal transmission system largely 


depends. 
To formulate the problem of the building-up of alternating 


currents, let us suppose that an e.m.f. # cos wt is suddenly applied 
at reference time ¢t = 0, to a net work or transmission system of 


complex steady state impedance 
Z (iw) =Z|\ (iw) |e? 
It is then always possible to write the resultant current in the 
form, 


ie : é 
NO AGE {1 +p) cos (wt— B)+o sin (wt — B)} (342) 
ey meen reer ne = e 
on ara Ziw) °° (wt — B(w) — 8) (343) 
where o=tan 


In these formulas p, o and @ are time functions which it is our 
problem to evaluate. Evidently p and o must be —1 and 0 
respectively for negative values of time t, and approach the limit 
1 and 0 respectively as [> ©. 

The function ay (1+ p)?+o: of equation (343) is proportional 
to the envelope of the oscillations. In an engineering study of 
the building-up process we are often not concerned with the 
instantaneous values of the current but merely with their 
envelope. In such cases an investigation of this function is 
sufficient. 

The mathematical formulation of p and o as Fourier integrals 
is based on formula (340). Starting with this formula as repre- 
senting a suddenly applied sinusoidal e.m.f., cos wt, of frequency 
w/2m, and proceeding by reasoning strictly analogous to that 
employed in the preceding section, it can be shown that p and 
o are given by the formulas 


»=2 {| Piper ne = (344) 
‘ 0) 
=| 12-0) -@.(-»} eosin 2 
cai (0.0)+0.(-») sina (345) 
0 


it? a 
+i [P.OQ)—P.(-0)} cos 0 O 
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where 
FLO nes a Cota (346) 
O.O0e: ih. boeeey tae (347) 


By virtue of the fact that p and o are identically —1 and 0 
respectively for negative values of time ¢, it follows at once (See 
derivation of (332) and (833) from (330), that we NEN also for 
positive values of time (t>0) 


I+p =3\| (Pa) +Pal—r)j sin tr = (348-a) 
‘JO 
I~p =| (Qa(A) —Qu(—A)} cos tr = (348-b) 
0 
and 
ee te red dn 
c= | (Q.0)+@.(—A)} sin’ \- — (349-a) 
0 
=f aN) ee Np COs: tN ss (349-b) 
La r 


For important types of transmission systems, particularly 
the periodically loaded line, these formulas have been success- 
fully dealt with and solutions of a satisfactory approximate 
character obtained. For further details the reader is referred to 
a paper on “The Building-up of Sinusoidal Currents in Long 
Periodically Loaded Lines”’ (Bell System Technical Journal, Oct., 
1924). 

The Distortion of Signals in Communication Circuits.—In 
electrical communication, the transmission system must be so 
designed that the received signal must be a more or less faithful 
copy of the transmitted signal; the departure from this require- 
ment is termed distortion. The Fourier integral provides a 
simple method of formulating the conditions for distortionless 
transmission, and the factors on which distortion depend. 

Let f(t) be the voltage impressed on the sending end of a trans- 
mission system of transfer impedance 


Z(iw) = Z\(iw)|e?, 


and let us suppose that f(¢) is expressed as the Fourier integral 


10 = {| Slee CaO (350) 
(0) 
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where w= 27 times the frequency. The received current is then 
given by 


dels (ae Ca) 
1o=5[ Zia ©98 (ot 8(s) — Bo) (351) 
Now f(t) of formula (350) contains components of all fre- 
quencies from zero to infinity; in other words F'(w) exists for all 
values of w. Practically, however, F(w) is negligibly small 
except in a range depending on the character of the signal. 
For example in slow speed telegraphy, the essential frequencies 
lie in a band between zero and 10 to 20 cycles per second. In 
telephony the essential frequencies occupy a band extending 
from a low frequency, say 50 ¢.p.s., to an upper limiting fre- 
quency of approximately 4000 e.p.s. 
Let us suppose, therefore, that the essential frequency range 
of f(t) lies between w;/2r and w./2x. Then equations (350) 
and (351) can be replaced by 


i == | “Roy don (EBON) (352) 
Le) =F ort (wt —6(w) — B(w)). (353) 


Now let us suppose that the transmission system is so designed 
that over the range w1<w<we 

|Z (iw)| = K 

B(w) =or. 


where K and 7 are constants. Formula (353) then becomes 
ioe it FG) cos Gor hay (354) 
~ 7h ie 


Comparison of (852) and (354) shows at once that 
1 
I(t) =ailt—7). (355) 


That is, the received current has precisely the same wave form 
as the impressed e.m.f.; it is simply changed in amplitude by 
the factor 1/K, and delayed in time by the “transmission time” r. 
The received signal is therefore a faithful copy of the transmitted 
signal and the transmission is distortionless. We thus arrive 
at the following important proposition. 
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The necessary and sufficient conditions for the distortionless 
transmission of signals in communication systems is that, over the 
essential range of frequencies contained in the transmitted signal, 
the transfer impedance of the transmission circuit be equalized both 
as regards amplitude and phase; that is, the amplitude must be 
constant and the phase angle linear in the frequency. 

The case where the amplitude alone is equalized over the 
frequency range corresponding to the spectrum of the signal is 
covered by the following proposition. 

The necessary and sufficient condition that the transmitted and 
recewed signals have the same relative energy content with respect 
to frequency is that the transmission system be equalized as regards 
amplitude over the range of frequencies contained in the signal. 

It is to be observed that this condition alone, however, does 
not insure preservation of wave form. In types of signalling 
systems, therefore, where preservation of wave form is essential, 
it is necessary to equalize both phase and amplitude. 

The Energy Absorbed from Transient Applied Forces.—In 
many technical problems, the complete solution for the instan- 
taneous current due to suddenly applied electromotive forces, 
although formally straight-forward, involves a_ prohibitive 
amount of labor. In yet others, the applied forces may be 
random and specified only by their mean square values. In 
such problems a great deal of useful information is furnished by 
the mean power and mean square current absorbed by the net- 
work, and to the calculation of these quantities, the Fourier 
integral is ideally adapted. Its application depends on the 
following proposition, due to Rayleigh (Phil. Mag., Vol. 27, 
1889, p. 466), and its corollary. 

Let a function ¢(t), supposed to exist only in the epoch 0<t< T, 
be formulated as the Fourier integral 


o=_ | ise) . cos [wt — 0(w)|dw 
0 


where 


fe) = | [so cos wt ul +| [oo sin wt ue’ 


\ 


T T 
tan 0(w) = [ (t)sin wt a+ (t) cos wt dt. 
0 0 


Then 
Lh i 
[ecra=_ | se) Me 
0 0 
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whereby the time integral is transformed into an integral with 
respect to frequency. 

A corollary of this theorem is as follows: 

If two functions ¢1(t), ¢2(t) supposed to exist only in the epoch 
0<t<T, are formulated by the Fourier integrals 


gilt) = [ne -cos [wt —61(w) |dw 
J0 
d() =F | falo)|-coslut—0x(u)lde, 


e (0) 
then 


ea Pes , 
[, e@eutat=_ | fle) -[falw)|- cos @s—0s)de, 
0 0 


The application of these theorems to circuit theory proceeds 
as follows :— 

If an electromotive force E(t), supposed to exist only in the 
epoch 0<t<T, is applied to a network of complex impedance 
Z (iw) = Z(iw)\e®, we know from the preceding discussion of the 
Fourier integral, that the electromotive force H(t) and current 
I(t) are expressible as the Fourier integrals 


E(t) = {| wifi one (are 0 Gada (356) 
0 

1” lf) Ms 

I(t) =| [Z (ies) COs (wt — A(w) B(w))dw. 
It follows at once from age theorem that 
as ” \f@)? 
Pat== peed 
Hl , Gaye w. (357) 


Now let J, be the current absorbed in branch n; let 2(iw) = 
lJe(iw) |e be the impedance of that branch and let Z,,(¢) be the 
potential drop across that branch. It follows at once from the 
corollary to Rayleigh’s mae that 


We I. B.()i,(pdiae | Ol Gay cosetaeecas 
0 ™ Jo |Z(tw)° 

Formulas (357) and (358) formulate the mean square current 

and mean power absorbed by the branch of the network under 

consideration, and enable us to calculate these quantities, even 

in the case of complicated networks, with a minimum of labor. 

Formula (357) is particularly well adapted to computation 
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because the integrand is everywhere positive, permitting, in 
most problems, of easy numerical integration, whereas the 
analytical solution may be complicated. 

Selective Circuits and Interference.—An important applica- 
tion of the preceding formulas is to the theory of selective circuits 
with particular reference to the protection they afford against 
irregular and transient disturbances such as “‘noise”’ in telephone 
systems and “‘static interference” in radio systems. In these 
cases the disturbance or interference is essentially random in 
character; its instantaneous wave form is unknown, and it can 
be specified only by its “frequency distribution.’’ Under these 
circumstances the problem is a statistical one and the only 
information deducible is the mean power and mean square current 
absorbed by the circuit from the disturbance. This is calculable 
from the preceding formulas when the statistical or average 
frequency distribution of the disturbance is determined or 
specified. A brief sketch of the analytical theory follows. 

Let a disturbance f(t), for convenience supposed to exist only 
in the epoch 0<t<7, be impressed on a network of impedance 
Z(iw). If the disturbance is specified for all values of ¢ in this 
epoch, it is theoretically possible to calculate F(w) and 6(w) and 
to express f(t) as the Fourier integral 


f= = il Gees LID (359) 


In any case even if insufficient information is supplied to calcu- 
late F(w) and 6(w) we know that f(t) is expressible as an integral 
of this form. From formula (358) we have 


| fal n| Pa 
i ae w, 
ia art h 


Now let the epoch 7 become indefinitely large; if the disturbance 
is random in character, the left hand side of this equation will 
approach a limit, equal to the mean square value of f(t), and we 
can write 

1 


fu=2 | R(w)dw (360) 
T Jo 
where R(w) is the limit of 7 (w) as T’ becomes very large. 


R(w) will be termed the statistical frequency spectrum of the 
disturbance. 
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Corresponding to the above the mean square current absorbed 
by the network is given by (See equ. (357)) 


iy (eA) 
rJo 2 (tw)? 


Formula (361) suggests a method for the determination of the 
frequency spectrum of the disturbance as follows: Let us suppose 
that the network is a selective circuit which absorbs current only 
in a narrow frequency band corresponding to w;<w<we Then 
if w2—w is made sufficiently small, formula (361) becomes 


Ph es 
I= 


dw. (361) 


il ae 
phe ee 
= Ron) Z(iw)|? (362) 
where w,, 1s defined by the equation 
= w.da Gs 
n= ee sos 363 
om Jy [2Gw)? i Z (iw)? Bee 


It follows from the preceding that measurement of the mean 
square current and calculation or measurement of the impedance 
characteristics of the selective circuit supplies the information 
necessary to compute the statistical frequency spectrum of the 
disturbance. 

Returning to formula (361) let us suppose that the selective 
circuit is designed to receive a band of frequencies w1<w<w, 
and excludes all other frequencies. The mean square current 
is then given by 


ees “#2 Rw) re 
= r J. Zea (364) 
Sauipe 365) 
a (@ iby IZ (tw) 2 (3 ao 


where w1<w'<wy». 

If therefore the selective circuit is required to receive a band 
of frequencies corresponding to w;<w<ws, formulas (364) and 
(365) give the irreducible minimum of interference. That is, the 
selective circuit must necessarily receive the energy of the dis- 
turbance which lies in the frequency range it is designed to receive. 
A definite limit is therefore set to the protection against inter- 
ference which can be obtained by means of selective circuits. 

This limitation is inherent in the nature of the interference 
and the frequency requirements of the signalling system. 
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For further and more detailed application of formulas (357) 
and (358) the reader is referred to the following papers: 

Transient Oscillations in Electric Wave Filters, Bell System 
Technical Journal, July, 1923. 

Selective Circuits and Static Interference, Trans. A. I. E. E., 
1924. 

An Application of the Periodogram to Wireless (Burch and 
Bloemsma), Phil. Mag., Feb., 1925. 

The Theory of the Schrotteffekt (Fry), Journal Franklin 
Institute, Feb., 1925. 

The foregoing must conclude our very brief account of the 
Fourier Integral and its applications in Electric Circuit theory; 
an adequate treatment of this subject would require a treatise 
in itself, and is beyond the scope of the present work. All that 
has been attempted is to give a very brief introduction to its 
significance in physical problems and a few of its outstanding 
applications in circuit theory. The reader who is interested in 
pursuing this subject further is referred to a paper by T. C. Fry 
on “The Solution of Circuit Problems” (Phys. Rev., Aug., 
1919), which gives a rigorous discussion of the solution of the 
Fourier Integral by contour integration, together with some 
general forms of solution of the circuit problems. 
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